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' We carry out an extensive investigation of conservation laws and potential symmetries for the 

. class of linear (l + l)-dimensional second-order parabolic equations. The group classification 

of this class is revised by employing admissible transformations, the notion of normalized 



(N 



^ , classes of differential equations and the adjoint variational principle. All possible potential 

conservation laws are described completely. They are in fact exhausted by local conservation 
laws. For any equation from the above class the characteristic space of local conservation laws 
' is isomorphic to the solution set of the adjoint equation. Effective criteria for the existence 

of potential symmetries are proposed. Their proofs involve a rather intricate interplay 
(— I ' between different representations of potential systems, the notion of a potential equation 

^ |. associated with a tuple of characteristics, prolongation of the equivalence group to the whole 

potential frame and application of multiple dual Darboux transformations. Based on the 
tools developed, a preliminary analysis of generalized potential symmetries is carried out and 
then applied to substantiate our construction of potential systems. The simplest potential 
symmetries of the linear heat equation, which are associated with single conservation laws, 
are classified with respect to its point symmetry group. Equations possessing infinite series 
^■f^ , of potential symmetry algebras are studied in detail. 

> 

^ ! 1 Introduction 

. In the present paper we classify local and potential conservation laws and potential symmetries 

' of linear (1 -|- l)-dimensional homogeneous second-order parabolic equations of the general form 

O ■ 

ut = A{t, x)uxx + B{t, x)ux + C{t, x)u, (1) 

I where A = A{t, x), B = B{t, x) and C = C{t, x) are arbitrary smooth functions, A ^ 0. 

^ • This class contains a number of physically important subclasses that are widely investigated 

and may be applied in many situations. Probably the most famous examples are the Kolmogorov 
equations (C = 0) and adjoint to them the Fokker-PIanck equations {A^x — -|- C = 0) which 
are often considered as main equations of continuous Markov processes. Fokker-Planck equa- 
tions appeared first in [27] where the Brownian motion in the radiation field was studied, and 
in [52j where one of the first systematic attempts of constructing a complete theory of fluctua- 
tions was made. They are also derived from the Boltzmann equation in the limit of large impact 
parameters [H]. The rigorous mathematical substantiation of the Fokker-Planck equation in 
the framework of probability theory was given by Kolmogorov 02]. Now the Fokker-Planck 
equations form a basis for analytical methods in the investigation of continuous Markov pro- 
cesses. Fokker-Planck equations with different coefficients describe the evolution of one-particle 
distribution functions of a dilute gas with long-range collisions, Brownian motion without drift, 
problems of diffusion in colloids, population genetics, financial markets, quantum chaos, etc. 

[231 [301 ED Eg. 

Local conservation laws of linear parabolic equations are, indeed, well understood. More 
precisely, it is well known that the space of characteristics of the local conservation laws of 
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a linear system of partial differential equations comprises those functions of the independent 
variables which solve the adjoint system. However, to the best of our knowledge, the statement 
that all characteristics are equivalent to such functions was proved, in an explicit way, only for 
the linear heat equation [JT] so far. Below this statement is extended to the entire class ([1]). 
Moreover, it is proved that linear second-order parabolic equations have no purely potential 
conservation laws. In other words, for these equations potential conservation laws of any level 
are exhausted by local conservation laws. This generalizes the analogous statement from [60] 
on the linear heat equation and finally solves the problem on potential conservation laws in 
class ([T|). 

Let us emphasize that in fact it is not a common situation when all the characteristics 
of the local and, especially, potential conservation laws of a linear system can be taken as 
functions of the independent variables alone, i.e., the associated conserved vectors are linear 
in the unknown functions and their derivatives. Even the usual wave equation and third-order 
parabolic equations have characteristics essentially depending on derivatives. An example of a 
third-order parabolic equation with such characteristics is given in Note [13] of the present paper. 
Moreover, a self-adjoint linear system possessing a quadratic conservation law has an infinite 
series of such conservation laws. See, e.g., the section on symmetric linear systems in Chapter 5 
and the concluding remarks of the corresponding chapter in ^B] . 

In contrast to conservation laws, potential symmetries have not been sufficiently investigated 
even for simplest equations from class ([T]). Thus, e.g., for the linear heat equation potential 
symmetries were studied only in the case of the single characteristics 1 |T2l [TO] [61] and x [38] . 
In [38l[63llS7] potential symmetries of the Fokker-Planck equation ut = Uxx + {xu)x, associated 
with the characteristic 1, were found. First-order conservation laws of the Fokker-Planck equa- 
tions of the form ut = Uxx + {B(x)u)x and potential symmetries of such equations, associated 
with the characteristic 1, were investigated in [63]. Note that the idea of symmetry extension via 
involving potentials and pseudopotentials in the transformations as new (nonlocal) dependent 
variables was already presented in the monograph by Edelen [22j . The concept of potential 
symmetry was explicitly formulated first by Bluman et al and was subsequently applied in in- 
vestigations of important classes of partial differential equations [12^ I13j. The related notion of 
quasilocal symmetry was proposed in [Tl[2]. A systematic procedure of constructing quasilocal 
symmetries of (l+l)-dimensional evolution equations was described in [SJ [81] and its applica- 
tion was illustrated by nontrivial examples. This procedure is based on the exhaustive group 
classification of such equations. To the best of our knowledge, the problem of finding criteria 
for the existence of potential symmetries for classes of differential equations was first posed by 
Pucci and Saccomandi [63] , 

The problem of a complete description of potential symmetries is very difficult to solve not 
only for classes of equations but even for single equations. In particular, it includes studying 
symmetry properties of infinite series of potential systems associated with tuples of an arbitrary 
number of linearly independent characteristics. We propose effective criteria for the existence 
of potential symmetries of equations from class ([l|) . They are subsequently applied to the clas- 
sification of simplest potential symmetries of the linear heat equation and the separation of 
subclasses possessing infinite series of purely potential symmetry algebras. The framework of 
potential symmetries in class ([I|) appears to be closely related to the theory Darboux transfor- 
mations [36] in the same class. 

Our paper is organized as follows: Using as a guideline the notion of normalization of classes 
of differential equations, in Section[2]we review and extend the classical results on Lie symmetries 
and equivalence transformations of class ([T|), introduced by Lie [l5] and Ovsiannikov [51]. In 
particular, difficulties arising under group classification of the Kolmogorov and Fokker-Planck 
equations are satisfactorily explained in terms of normalized classes. The necessary definitions 
and statements on conservation laws and their characteristics, the equivalence of conservation 
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laws with respect to transformations groups and sets of admissible transformations and potential 
systems are collected in Section [3] for convenient reference. Special attention is paid to the 
equivalence of conservation laws with respect to transformation groups (Subsection IS.Sp as this 
notion is essential in studying the potential conservation laws of equations from class ([1]). 

The local conservation laws of these equations are exhaustively described in Section [5] based 
on the direct method. Since for any equation under consideration the characteristic space of 
local conservation laws is isomorphic to the solution set of the adjoint equation, in Section [5] 
the adjoint variational principle is specified for class ([I]) and extended to the corresponding set 
of admissible transformations. A number of auxiliary statements on admissible transformations 
of second-order evolution systems is proved. The adjoint variational principle is then applied to 
the group classification of the Fokker-Planck equations. 

The main result on potential conservation laws in class ([T|) is presented in Section [6j Namely, 
it is proved that the local conserved vectors of potential systems are equivalent to local conserved 
vectors of the corresponding equations. This also provides a complete description of potential 
systems, allowing us to initiate the investigation of potential symmetries. 

The simplest potential symmetries considered in Section [7] form a subject whose investigation 
generates a number of ideas on a special technique of working with class ([1]). The attribute 
'simplest' refers to the fact that these symmetries are associated with single characteristics (as 
opposed to simplicity of calculation). Thus, in Section [8] the simplest potential symmetries of 
the linear heat equation are classified with respect to its point symmetry group. There are only 
two inequivalent characteristics a = 1 and a = x giving simplest purely potential symmetries 
of the linear heat equation. This is the only example in the literature with an exhaustive 
investigation of at least simplest potential symmetries of equations from class ([1]). The obtained 
classification directly leads to the complete description, e.g., of the simplest second- level potential 
symmetries of the Burgers equation and the simplest potential symmetries of the equations 
which are equivalent to the linear heat equation with respect to point transformations. The 
classification of simplest potential symmetries of the Fokker-Planck equation ut = Uxx + ixu)x 
is presented for illustration. 

Section [9] is devoted to the construction of the potential frame over class ([I]) and studying 
its (non-symmetry) properties. Different kinds of potentials, potential systems and potential 
equations (p-order and p-level, usual and modified ones) associated with characteristic tuples 
are defined. Explicit expressions for all these object are found. The multiple dual Darboux 
transformation provides a clear connection between components of the potential frame. This is 
why enhanced statements on Darboux transformations in class are also presented. Probably 
the most important components of the potential frame are the so-called modified potentials 
and modified potential equations. In contrast to other components, they are invariant under 
nonsingular linear combining of characteristics in the associated tuples. Moreover, as proved 
in Section \W\ Lie symmetry analysis of potential systems is reduced to group classification of 
modified potential equations with respect to the equivalence group of class ([T]) prolonged to the 
whole potential frame. Another result of this section is the possibility and precise realization of 
the prolongation. A statement on generalized potential symmetries of equations from class ([1]) 
is also proved. This serves to substantiate the use of the canonical form of the conserved vectors 
in the construction of the potential frame. 

After analyzing Lie invariance of potential systems, we formulate, in different terms, cri- 
teria on the existence of general potential symmetries. Their effectiveness is demonstrated in 
Section [11] via the construction of wide subclasses of class ([1]) whose equations admit infinite 
series of potential symmetry algebras of arbitrarily large order. The multiple auto-Darboux 
transformation is used in this construction as a powerful auxiliary tool. 

In the final section the results of the paper are summarized and some open problems on 
potential symmetries of equations from class ([T]) are formulated and discussed. 
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We will refer to formula ([T]) and similar ones describing classes of systems of differential 
equations in a twofold manner, namely either as to a whole class (then the arbitrary elements 
are assumed to run through all possible values) or as a single equation from this class (then the 
arbitrary elements are assumed to take fixed values). 

By default, the indices i, j and k run from 1 to n, the indices a and b run from 1 to m and 
the indices s, a and ? run at most from 1 to p. Additional or other constraints on indices are 
indicated explicitly. The summation convention over repeated indices is used unless otherwise 
stated or it is obvious from the context that indices are fixed. 

2 Group classification 

The complete group classification of equations ([T|) was performed by Sophus Lie ^45j as a part 
of the more general group classification of linear second-order partial differential equations in 
two independent variables. A modern treatment of the subject is given in [51]. There exist 
also a number of papers rediscovering results of Lie and Ovsiannikov [45^1511 partially (see, e.g., 

[ini [la ia Eg Eg ED ESI [731 EH ) . 

Since our investigation on potential symmetries of equations from class ([T]) is essentially 
based on the above results, we review them for the reader's convenience. Another reason of 
their consideration is to justify the choice of a suitable subclass whose investigation allows 
the description of potential symmetries and potential conservation laws in the whole class ([1]). 
Moreover, we extend these results to the framework of admissible transformations in classes of 
differential equations. The normalization properties of the class of linear parabolic equations 
and its subclasses with respect to point transformations are also studied. 

Roughly speaking, an admissible transformation in a class of systems of differential equations 
is a point transformation connecting at least two systems from this class (in the sense that one 
system is transformed into the other by the transformation). The class is called normalized if 
any admissible transformation in this class belongs to its equivalence group and is called strongly 
normalized if additionally the equivalence group is generated by transformations from the point 
symmetry groups of systems from the class. The set of admissible transformations of a semi- 
normalized class is generated by the transformations from the equivalence group of the whole 
class and the transformations from the point symmetry groups of initial or transformed systems. 
Strong semi-normalization is defined in the same way as strong normalization. Any normalized 
class is semi-normalized. Two systems from a semi-normalized class are transformed into one 
another by a point transformation iff they are equivalent with respect to the equivalence group 
of this class. See [Ml ESI \57\ [58] for precise definitions and statements. 

To begin with, consider the class of inhomogeneous equations corresponding to ([T|) of the 
general form 

ut = A{t, x)uxx + B{t, x)ux + C(t, x)u + D{t, x), (2) 

where A{t,x), B(t,x), C{t,x) and D{t,x) are arbitrary smooth functions, A{t,x) 7^ 0. A moti- 
vation for the provisional extension of the class under consideration is that class ^ has a nicer 
normalization property than the initial class. 

Any point transformation T in the space of variables {t,x,u) has the form 

i=T\t,x,u), X = T''{t,x,u), n = r"(t,x,M), 

where the Jacobian \d{T^ ,T^)/d{t, x,u)\ does not vanish. 

Lemma 1. A point transformation T connects two equations from class ^ iff Tj: = = 0, 
T^ = T{t), T'^ = X{t,x), T"" = U\t,x)u + U'^{t,x), (3) 
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where T , X, and are arbitrary smooth functions of their arguments such that TtX^U^ 7^ 0. 
The arbitrary elements are transformed by the formulas 



Here L = dt — Adxx — Bdx — C is the second- order linear differential operator associated with 
the initial (non-tilde) equation. 

Proof. The proof is based on the direct method. We recalculate the old derivatives in the 
new variables, substitute the obtained expressions in the initial equation and then split it on 
the manifold of the transformed equation. As a result, we derive the determining equations 
for the components of T. The calculations can be simplified if we take into account known 
restrictions for transformations between equations from wider classes. Thus, equations of the 
form ([2]) are evolutionary. Any transformation between two evolution equations satisfies the 
conditions Tj: = = [41j. Then, a transformation between quasi-linear equations satisfies 
the conditions = and = (see, e.g., [591 [62]). The remaining determining equations 
implies the formulas for transforming the arbitrary elements. □ 

Corollary 1. Class ([2|) is strongly normalized. The equivalence group G'^^ of class ([2|) is 
formed by the transformations determined in the space of variables and arbitrary elements by 
formulas ([3]), ([1]) and ([5]), where T, X, and are arbitrary smooth functions of their 
arguments such that TtX^U^ 7^ 0. 

Using transformations from G^j^, we can gauge arbitrary elements of class ([2]). Thus, applying 
the equivalence transformation with T = t, X = x, = 1 and a solution of the equation 
LU^ = —D, we obtain the standard conversion of the inhomogeneous equation Lu = D to the 
homogeneous one Lu = 0. As a result, class ([2]) is mapped to class ([T]). Unfortunately, the 
normalization property is broken under this mapping. 

Corollary 2. A point transformation T connects two equations from class ([TJ iff its components 
are of the form ([3|), where T, X and are arbitrary smooth functions of their arguments such 
that TtXxU^ 7^ and additionally U^/U^ is a solution of the initial equation. The arbitrary 
elements are transformed by formulas (Uj) . 

Corollary 3. Class ([T|) is strongly semi-normalized. The equivalence group G~ of class ([1]) is 
formed by the transformations determined in the space of variables and arbitrary elements by 
formulas ([3]), where T, X and are arbitrary smooth functions of their arguments such 
that TtXxU^ 7^ and = additionally. 

Note 1. A similar relation holds between normalization properties of general classes of inhomo- 
geneous and the corresponding homogeneous linear systems of differential equations. Namely, 
consider a class /2hom of homogeneous linear systems of / differential equations of the form 
LffU = 0, for m unknown functions u = (n^, . . . , u^) of n independent variables x = (xi, . . . , Xn). 
Here Lg = (Lg"), fi = 1, . . . ,1, is a matrix differential operator parameterized with 6 running 
through a parameter set. Let the corresponding class Ci^h of inhomogeneous systems be nor- 
malized and its equivalence group G~jj consist of transformations projectable in x and affine 
in u, i.e. the transformations of {x,u) have the form x = X[x) and u = U"'^{x)u^ + U°'^{x) 
for any T G G~j^. Then the class £hom is semi-normalized. The equivalence group G^^^ of 
^hom is isomorphic to the subgroup of G^j^ formed by the transformations with (C/"*C/"°) run- 
ning through the intersection of the solution sets of the systems from Chom, where {U"'^) is the 
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inverse matrix of {U""^). Often this implies that U""^ = 0. The additional condition for the 
admissible transformations in £hom is that {tj°'^U^^) is a solution of the initial system (with 
fixed values of the arbitrary elements). Therefore, normalization is broken under restricting to 
>Chom due to the presence of the linear superposition principle. This justifies the consideration 
of inhomogeneous linear systems in the framework of admissible transformations. 

Another possibility is to gauge the arbitrary element A in class ([2]) to 1 with a transformation 
of form ([3]), where Tt = sign A, Xx = \A\~^/'^^ C/^ = 1 and = 0. The admissible transfor- 
mations in the subclass of ^ with A = \ are those transformations ^ which preserve the 
condition A = 1, i.e., which additionally satisfy the condition = (T^)'^. 

Corollary 4. A point transformation T connects two equations from class ([2|) with A = A = 1 iff 

T' = T{t), T'^ = X = ±y^T^x + C{t), T"" = U\t,x)u + U^{t,x), (6) 

where T, ^, and are arbitrary smooth functions of their arguments such that Tt > and 
7^ 0. The transformations of this form, prolonged to the arbitrary elements B, C and D by 
formulas ^ and ([5]) constitute the equivalence group of the subclass of ([2]) with A = 1. This 
subclass is strongly normalized. 

Analogously to the case of the entire class ([2]), we can convert the inhomogeneous equations 
with ^ = 1 to the homogeneous ones. As a result, the subclass of ([2]) with A = 1 is mapped to 
the subclass of ([T]) satisfying the same condition. The normalization property is again broken 
under this mapping. A point transformation T connects two equations from class ([T]) with 
A = A = 1 iS it has the form adduced in Corollary [J] and additionally U^/U^ is a solution of 
the initial equation. The subclass of ([T]) with A = 1 is strongly semi-normalized. Its equivalence 
group is formed by the transformations ([6]) with = 0, prolonged to the arbitrary elements B 
and C by formulas 

The arbitrary elements A and B can be simultaneously gauged to 1 and 0. The subclass 
of ([2]) with (A, B) = (1, 0) is a restriction of the one with A = 1 and is investigated in a similar 
way. 

Corollary 5. A point transformation T connects two equations from class ([2]) with A = A = 1 
and B = B = iff it has the form where additionally 

U' = 9{t) exp ^ -\^x] 

and T, 9 and are arbitrary smooth functions of their arguments such that Tt > and 9^0. 
The transformations of this form, prolonged to arbitrary elements C and D by formulas (H]) 
and dl]) constitute the equivalence group of the subclass of ([2]) with A = 1 and B = 0. This 
subclass is strongly normalized. 

The inhomogeneous equations with A = 1 and B = are mapped to the homogeneous ones 
in the standard way. Hence, any equation from class ([2]) or class ([TJ can be reduced by a 
transformation from the corresponding equivalence group to an equation of the general form 

ut - Uxx + y{t,x)u = 0. (7) 

The normalization property is broken for the class ([7]). A point transformation T connects 
two equations from class d?]) iff it has the form adduced in Corollary [5] and additionally U^/U^ 
is a solution of the initial equation. At the same time, class ([7]) is strongly semi-normalized. 
Its equivalence group is formed by the transformations with = 0, prolonged to the 
arbitrary elements V = —C by formulas (j4]). Therefore, the functions parameterizing Gi depend 
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only on t. The narrower equivalence group under preserving certain normalization properties 
suggests class ([7]) as the most convenient one for group classification. Moreover, solving the group 
classification problem for any of the above classes is reduced to solving the group classification 
problem for class d?]). This is why we formulate the main result on admissible transformations 
in class ([7]) theorem. 

Theorem 1. Class ([7]) is strongly semi-normalized. Any transformation from the equivalence 
group Gi of class (j7]) has the form 



where a = (T{t), C, = C,{t) and 9 = 9{t) are arbitrary smooth functions, a9 ^ 0. 

For our further considerations we need to introduce terminology connected with the symmetry 
structure of linear equations. In view of the linear superposition principle, the point symmetry 
group G{C) and the maximal Lie invariance algebra q{C) of any (homogeneous) linear differen- 
tial equation (or any system of such equations) C have certain properties [5l]- Namely, G{C) 
contains the translations u ^ u + ef of the unknown function u by an arbitrary solution f of C 
and the scale transformations u ibe^-u superimposed with reflection of u. These sets of trans- 
formations form subgroups of G{C) both separately and simultaneously. The whole subgroup 
of symmetry transformations associated with the linear superposition principle of C is called 
the trivial symmetry group of the linear differential equation C and will be denoted by G^"^{C). 
The algebra q{C) contains the corresponding operators fdu and udu which form the Lie algebra 
g*"^(£) called the trivial invariance algebra of the linear differential equation C Hereafter the 
function / runs through the solution set of C 

The ideal Q°°{J~.) = {fdu) of g^'^"{C) is cafled the (trivial) inflnite-dimensional part of q{C). 
The corresponding normal subgroup of G{C) is similarly denoted by G^{C). This notation is 
justified in the following way: Usually [51] the maximal Lie invariance algebra %{C) of a linear 
differential equation C can be represented in the form 



Here q°°{C) is an (infinite-dimensional) Abelian ideal of the algebra q{C) and 0'^^^(£) is its finite- 
dimensional subalgebra spanned by the Lie invariance operators of £, which are projectable to 
the independent variables and whose coefficients of 5„ depend linearly on u. In particular, the 
above representation is true for the equations from class ([1]). Then a similar representation 
also holds for the group G{C). Namely, G{C) = G''''\C) x G°°{C), where G^{C) is a normal 
subgroup of G{C) and G^^^{C) its subgroup. 

The operators from q'^^^{C) will be called essential symmetry operators of the linear differential 
equation C since mainly they and the corresponding finite transformations are useful for group 
analysis of C Such operators can be found in a particular way via the commutation relation with 
the differential operator associated with the equation C [28] (see also Section [5]) that provides 
possibilities for various generalizations of the notion of symmetry operators. They are employed 
for finding finite symmetry transformations of £, for constructing exact solutions via the Lie 
reduction procedure to differential equations with fewer independent variables and for the direct 
generation of new exact solutions by acting on known ones [51]. The algebra q'^^^{C) will be 
called the essential Lie invariance algebra of the linear differential equation C 

The intersection q"^^ {C)(^q^'^" {C) = {udu) is contained in the center of q'^^^{C). The dimension 
of Q^^^ {C) / {udu) is called the dimension of extension of the maximal Lie invariance algebra or 
the number of independent nontrivial Lie symmetry operators. 




) 



(8) 



fl(/:) = 0-«(£)5 0-(£). 
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In terms of these notations we can formulate Corollary [3] and similar statements more pre- 
cisely. Any point transformation between two equations from class ([TJ is the composition of a 
trivial symmetry transformation from G°° of the initial equation and a transformation from G"" . 

The results on the group classification of class ([Tj) can be formulated in the form of the 
following theorem [451 EI] • 

Theorem 2. The kernel Lie algebra of class ([7]) is {udu)- Any equation from class ([7]) is 
invariant with respect to the operators fdu, where the parameter- function f = f{t,x) runs 
through the solution set of this equation. All possible Gi -inequivalent cases of extension of the 
maximal Lie invariance algebra are exhausted by the following ones (the values of V are given 
together with the corresponding maximal Lie invariance algebras): 

1. V = V{x): {dt, udu, fdu); 

2. V = /I / 0: {dt, D, n, udu, fdu); 

3. v = o: {dt, d^, a, D, n, udu, fdu). 

Here D = 2tdt + xd^, U = At^dt + Uxd^ - (x^ + 2t)udu, G = 2td^ - xudu- 

Note 2. It is assumed in case 1 of Theorem [2] that the value V = V{x) is -inequivalent to 
the value /ux"^, where ^ € M. 

Note 3. Theorem [2] can be reformulated for the entire classes ([I]) and ([2]) if G^-equivalence is 
replaced by G~- and G~j^-equivalences correspondingly. A similar reformulation is possible also 
for subclasses with A = 1. Let us emphasize that the group classification in a semi-normalized 
class with respect to its equivalence group is identical to the classification up to all admissible 
point transformations. 

Corollary 6. For any equation C from class ([T|) dim0°'^'^(£) G {1,2,4,6}, i.e., the number of 
independent nontrivial symmetries belongs to {0, 1,3,5}. // dimg^****(/^) = 6 then the equation C 
is G"" -equivalent to the linear heat equation ut = Uxx- 

Note 4. The presented way of gauging the arbitrary elements is optimal for group classifica- 
tion. The hierarchy of normalized classes of inhomogeneous equations and the corresponding 
semi-normalized classes of homogeneous equations are constructed. Due to its properties, the 
subclass ([7]) is convenient for solving the group classification problem. The obtained results can 
be extended in an obvious way to all classes from the hierarchy. Different choices of gauges for 
the arbitrary elements (e.g., reduction to the 'Kolmogorov' or 'Fokker-Planck' form) may lead 
to a considerable complication of the problem. 

Consider the group classification problem for the 'Kolmogorov' form (C = 0) of equations 
from class ([1]). (It follows from results of Section [J] that the symmetry analysis of the 'Fokker- 
Planck' form, being adjoint to the 'Kolmogorov' form, is reduced to an investigation of the 
'Kolmogorov' form.) Note that the signs of A and B are inessential under symmetry investigation 
due to the presence of equivalence transformations alternating the signs. 

The gauge G = totally breaks the normalization properties. Indeed, a point transforma- 
tion T connects two equations from the class ([1]) with G = G = iff its components are of the 
form ([3]), where T and X are arbitrary smooth functions of their arguments such that TtXx ^ 0, 
C/^ 7^ and additionally and U^/U^ are solutions of the initial equation. The arbitrary 

elements A and B are transformed by formulas The equivalence group of the subclass 
of ([1]) with G = consists of only those transformations of the form 1^ with U^,U^ = const. 
Therefore, this subclass is not semi-normalized. There exist equations in it, transformed into 
one another by a point transformation, which are inequivalent with respect to the equivalence 
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group. The structure of admissible transformations which are not generated by transformations 
from the equivalence group is quite complicated. That is why it seems too difficult to present a 
classification for the subclass with respect to its equivalence group. The additional gauge A = 1 
does not improve the situation. A classification up to its set of admissible transformations is 
derived from Theorem [2] by mapping the listed equations to the 'Kolmogorov' form. 

Corollary 7. The kernel Lie algebra of the subclass of ([1]) with C = is {udu,du)- Any 
equation from this subclass is invariant with respect to the operators fdu, where the parameter- 
function f = f{t,x) runs through the solution set of this equation. By a point transformation it 
is reduced to an equation with A = 1 from the same subclass. All possible cases of extension of 
the maximal Lie invariance algebras in this subclass are exhausted, up to point transformations, 
by the following ones (in all the cases A = 1; the values of B are given together with the 
corresponding maximal Lie invariance algebras): 

1. B = B{x): {dt, udu, fdu); 

2. B = ux~^, u ^ 1, u ^ 2: {dt, D,Il — 2utudu, udu, fdu); 

3. B = — 2xtan(xln|x|)), {dt, 2D — xBudu, H — 2txBudu, udu, fdu); 

4. B = 0: {dt,d^,G, D,U,udu, fdu). 

Here D = 2tdt + xd^, H = At^dt + Uxd^, - (x^ + 2t)udu, G = 2td^ - xudu- 

Corollary [7] shows that even up to all admissible transformations, the group classification of 
equations in the 'Kolmogorov' form is more complicated than the group classification in class ([7]). 
Parameters of equations are explicitly included in expressions of symmetry operators. Case 2 
of Theorem [2] is split into two cases of Corollary [7] (case 2 with v = 1 + ^Jl + 4/i if 4^ ^ — 1 
and case 3 with x = y^— 1/4 — /i if 4/i < —1). These cases can be united only over the complex 
number. The equations with B = vx~'^ and B' = u'x~'^ {A = A' = 1) are equivalent with respect 
to point transformations iff z^ + i^' = 2. The corresponding transformation is u' = x'^~^u, t and x 
remain unchanged. In particular, the equation with B = 2x~^ is reduced by the transformation 
u' = xu to the linear heat equation {B = 0). That is why the parameter ly should be constrained 
in case 2. The form of the arbitrary element B is not simple in case 3 and cannot be simplified 
in the real case. Therefore, it is preferable to carry out the symmetry analysis of equations in 
the form ^ and then to derive results for the 'Kolmogorov' form. 

The group classification of the subclass with C = with respect to its equivalence group 
can be obtained from the classification presented in Corollary [7| by extending the classification 
cases by essential admissible transformations which are not generated by the equivalence group. 
Such transformations have the form i = t, x = x, u = U^(t,x)u, where ^ and is an 

arbitrary solution of the initial equations. 

3 Theoretical background on conservation laws 

To begin with, we present the necessary theoretical background on conservation laws and po- 
tential systems, basically following [TBI B9l [601 178]. 

3.1 Definition of local conservation laws 

Let C he a system L{x, U(^p^) = of / differential equations = 0, . . . , = for m unknown 
functions u = {u^, . . . , u^) of n independent variables x = {xi, . . . , Here U(p) denotes the set 
of all the derivatives of the functions u with respect to x of order not greater than p, including u 
as the derivative of order zero. Let denote the set of all algebraically independent differential 
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consequences that have, as differential equations, orders not greater than k. We identify 
with the manifold determined by in the jet space 

Definition 1 . A conserved vector of the system £ is an n-tuple F = {F^ (x , Uj-^^ ) , . . . , (x , Uj-^) ) ) 
for which the divergence DivF := DiF^ vanishes for all solutions of C, i.e., 

DivF|^=0. (9) 

In Definition [1] and below Di = D^, . denotes the operator of total differentiation with respect 
to the variable Xj, i.e., Di = dx^ + u'^^du'^, where u% and j stand for the variables in jet 
space which correspond to the derivatives d^"^u°'/dxi^ . . . dx'^" and du'^/dxi, a = (ai, . . . , a„), 
aj G N U {0}, \a\: = ai + • • • + an- We use the summation convention for repeated indices and 
assume any function as its zero-order derivative. The indices i, j and k run from 1 to n, the 
index a runs from 1 to m. The notation means that the values of V are considered only on 
solutions of the system C. 

Heuristically, a conservation law of the system C is an expression Div F vanishing on the 
solutions of C The more rigorous definition of conservation laws given below is based on the 
factorization of the space of conserved vectors with respect to the subspace of trivial conserved 
vectors. Note that there is also a formalized definition of conservation laws of £ as (n — 1)- 
dimensional cohomology classes in the so-called horizontal de Rham complex on the infinite 
prolongation of the system C \16\ \76\ [77] . The formalized definition is appropriate for certain 
theoretical considerations and reduces to the usual one after local coordinates are fixed. 

Definition 2. A conserved vector F is called trivial if F^ = F'' + F^, where F^ and -F* are, like 
F', smooth functions of x and derivatives of u (i.e. differential functions), F^ vanishes on the 
solutions of C and the n-tuple F = {F^, . . . , F") is a null divergence (i.e., its divergence vanishes 
identically) . 

The triviality concerning conserved vectors vanishing on solutions of the system can easily 
be eliminated by restricting to the manifold of the system, taking into account all its differential 
consequences. A (local) characterization of all null divergences is given by the following lemma 
(see e.g. [39]). 

Lemma 2. The n-tuple F = (F\...,F"), n ^ 2, is a null divergence (Div F = 0) iff there 
exist smooth functions v'-' of x and derivatives of u, such that v^^ = —v^^ and F^ = Djv^K 

The functions v^^ are called potentials corresponding to the null divergence F. If n = 1 any 
null divergence is constant. 

Definition 3. Two conserved vectors F and F' are called equivalent if the vector- function F' — F 
is a trivial conserved vector. 

If £ is a system of ordinary differential equations (n = 1) then the conserved quantities (first 
integrals) F and F' are equivalent by Definition [3] if their difference is constant on the solutions 
of £. 

In the case of two independent variables we re-denote {F^, F"^) {F, G) and (xi, X2) {t, x). 
Two conserved vectors {F, G) and {F', G') are equivalent if there exist functions F, G and H 
of t, X and derivatives of u such that F and G vanish on £(fc) for some k and 

F' = F + F + D^H, G' = G + G- DtH. 

The above definitions of triviality and equivalence of conserved vectors are natural in view 
of the usual "empiric" definition of conservation laws of a system of differential equations as 
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divergences of its conserved vectors, i.e., divergence expressions which vanish for all solutions of 
this system. For example, equivalent conserved vectors correspond to the same conservation law. 
It allows us to formulate the definition of conservation law in a rigorous style (see, e.g., [16ti78j). 
Namely, for any system L of differential equations the set CV(£) of conserved vectors of its 
conservation laws is a linear space, and the subset CVo(>C) of trivial conserved vectors is a 
linear subspace in CV(£). The factor space CL(£) = CV(£)/ CVo(/2) coincides with the set of 
equivalence classes of CV(£) with respect to the equivalence relation introduced in Definition [3l 

Definition 4. The elements of CL(£) are called conservation laws of the system £, and the 
factor space CL(£) is called the space of conservation laws of £. 

This is why we understand the description of the set of conservation laws as finding CL(£), 
which in turn is equivalent to constructing either a basis if dimCL(£) < cxd or a system of 
generators in the infinite dimensional case. The elements of CV(£) which belong to the same 
equivalence class giving a conservation law are all considered as conserved vectors of this 
conservation law, and we will additionally identify elements from CL(£) with their represen- 
tatives in CV(£). For F £ CV(£) and £ CL(£) the notation F e will mean that F is 
a conserved vector corresponding to the conservation law J^. In contrast to the order r^? of a 
conserved vector F as the maximal order of derivatives explicitly appearing in F, the order of 
the conservation law T is defined as minjri? | F G T^. By linear dependence of conservation 
laws we mean linear dependence of them as elements of CL(i2). Therefore, in the framework of 
the "representative" approach conservation laws of a system £ are considered linearly dependent 
if there exists linear combination of their representatives which is a trivial conserved vector. 

3.2 Characteristics of conservation laws 

Let the system £ be totally nondegenerate [49]. Then an application of the Hadamard lemma 
to the definition of conserved vector and integration by parts imply that the divergence of any 
conserved vector of £ can always be represented, up to the equivalence relation of conserved 
vectors, as a linear combination of the left hand sides of the independent equations from £ with 
coefficient functions on a suitable jet space 

DivF = A''L''. (10) 

Here the order k is determined by £ and the allowable order of conservation laws, [i = 1,1. 

Definition 5. Formula (jlOp and the Z-tuple A = (A^,...,A') are respectively called the char- 
acteristic form and the characteristic of the conservation law associated with the conserved 
vector F. 

The characteristic A is trivial if it vanishes for all solutions of £. Since £ is nondegenerate, 
the characteristics A and A satisfy (jlOp for the same F and, therefore, are called equivalent iff 
A — A is a trivial characteristic. Similarly to conserved vectors, the set Ch(£) of characteristics 
corresponding to conservation laws of the system £ is a linear space, and the subset Cho(£) of 
trivial characteristics is a linear subspace in Ch(£). The factor space Chf(£) = Ch(£)/ Cho(£) 
coincides with the set of equivalence classes of Ch(£) with respect to the above characteristic 
equivalence relation. 

The following result [l9] forms the cornerstone for the methods of studying conservation laws, 
which are based on formula (llOp . including the Noether theorem and the direct method in the 
version by Anco and Bluman [U [5] . 

Theorem 3 ([49j). Let C be a normal, totally nondegenerate system of differential equations. 
Then the representation of conservation laws of £ in the characteristic form ([TO]) generates a 
one-to-one linear mapping between CL(£) and Chf(£). 
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Using properties of total divergences, we can exclude the conserved vector F from (jlOp and 
obtain a condition for the characteristic A only. Namely, a differential function / is a total 
divergence, i.e. / = Divi*" for some n-tuple F of differential functions iff E(/) = 0. Here, the 
Euler operator E = (E^, • • • , E™") is the m-tuple of differential operators 



E" = (-Z))"5„a, a = l,m, 

where = (— Di)"^ . . . (— Dm)"™; a = (aij-'-jOn) runs through the multi-index set 

(ajENU {0}). Therefore, the action of the Euler operator on (jlOp results in the equation 

E(A^L'^) = DI(L) + DHA) = 0, (11) 

which is a necessary and sufficient condition on characteristics of conservation laws for the sys- 
tem C. The matrix differential operators and are the adjoints of the Frechet derivatives D;^ 
and D^, i.e.. 

Since D^(^) = automatically holds on solutions of >C, equation (fTTI) implies a necessary condi- 
tion for A to belong to Ch(£): 

D2(A)|^=0. (12) 

Condition (|12p can be considered as adjoint to the criterion \^i^{ri)\^= for infinitesimal invari- 
ance of £ with respect to an evolutionary vector field having the characteristic r] = [rj^, . . . , r/™). 
This is why solutions of ()12p are sometimes called cosymmetries or adjoint symmetries. 



3.3 Equivalence of conservation laws with respect to transformation groups 

We can substantially simplify and systematize the classification of conservation laws by addition- 
ally taking into account symmetry transformations of a system or equivalence transformations 
of a whole class of systems. This problem is similar to that of group classification of differential 
equations. The following statement on transformations of equations in the conserved form is 
true (see, e.g., [60]). 

Proposition 1. Any point transformation T maps a class of equations in the conserved form 
into itself. More exactly, the transformation T: x = T^{x,u), u = T^{x,u) prolonged to the 
jet space transforms the equation DiF^ = to the equation DiF^ = 0. The transformed 

conserved vector F = {x,U(jY F) is determined by the formula 

D ■ X- ■ ~ 1 

F\x,u,^))= F^{x,uir)), i-e. F{x,urr)) = 1^ ^ {D^x)F{x,mr)) (13) 

^ ' \F>xX\ ' ^ ' \F>xX\ ^ ' 

in matrix notation. Here \F)xx\ is the determinant of the matrix D^x = (DxjXi). 

Proof. We give two equivalent versions of the proof. The first is direct and based on the usual 
definition of conservation laws. The second is closer to the formal definition and naturally 
involves the technique of differential forms. 

We prolong the transformation T to the jet space in the standard way [l9] , i.e. , we recal- 

culate all derivatives up to order r+1 in the new ('tilde') variables: U(^r+i) = P^\r+i)'^^i^^ '"(r+i))- 
Since the equation DiF'^ = is linear in F, the transformation for the tuple F is found in the 
form F^ = G^^ FK The smooth functions G^^ of x and the derivatives of u should be selected by 
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the condition -Dj-F* = ADiF^, where A also is a smooth function of x and the derivatives of u. 
Let G = (G*-') be the inverse matrix to the matrix G = (G*-'), A = 1/A. Then 

DjF^ = Dj{&^F^) = {Dj&^)F^ + {DjXi)&^biF^ = ADiF' 

for any F iff {DjXi)&^ = A(5jfc, Dj&^ = 0. Here 6ik is the Kronecker delta. The first set of 
equations on &^ implies that G = A{Dxx)^^ = k{Dxx). Substituting these expressions for &^ 
into the second set of equations we get: 

DjG^^ = DjikbkXj) = Dj{k)bkXj + A{DjXk')bk'bkXj = 0. 

(Note that Dj = {DjXk')Dk'-) After dividing the result by A and convolving it with DiXk, we 
obtain 

-{DiXk){DjXk')bk'bkXj = -{DjXk'){DiXk)bkbk'Xj = -{DjXk')Dibk'Xj 



A 

= -ti{{Dxx)D,{Dx3:)-^) = tv{{D,Dxx){Dxx) 



\D,rX\ 



where we have used the commutation property of the total derivative operators and the well- 
known equalities for matrix derivatives {A~^)' = —A^^A'A^^ and tr(yl'^~-'^) = Here 
the prime denotes the derivative with respect to a parameter, ti A and |^| are the trace and 
the determinant of a square matrix A, respectively. The above equations imply that A = 
l-D^^iil up to an arbitrary nonzero constant multiplier (which is inessential) and, therefore, G = 
\Dxx\{Dxx)~^ , i.e., G = \Dxx\^^ DxX. 

The second version of the proof is much simpler and indeed justifies the first one. We associate 
any tuple F with the differential form lof = {—iy~^F^ dxiA- ■ • A&q, A- • ■f\dxn in the 'horizontal' 
de Rahm complex |16] called also D-complex [49j over the space of the independent variable x 
and the dependent variable u. Hereafter the notation c&jj means that the term dxi is absent 
in the corresponding external product. In the D-complex the differential of the usual de Rahm 
complex is replaced by the total differential D. Due to the invariance of differential forms under 
transformations of variables, 

ujp = {-ly^^F^ dxi A • • • A A • • • A dxn 

= {-ly-'^F'M'^ dxi A • • • A A • • • A dx„ = LOp, 

where M^^ is the minor of the element DjXi in the matrix D^x. Therefore, F^ = (— 1 M*-'^*, 
i.e., F = \Dxx\~^(Dxx)F. Applying the total differential, we also have Du;^ = Dljf, where 

Dlof = {DjF^)xi A • • • A dxn, 

Dcvp = {biF'') dxi A • • • A dxn = (A-^*)|-Dxx| dxi A • • • A dxn, 

i.e., biF^ = \Dxx\-^DjF^. □ 

Note 5. In the case of one dependent variable (m = 1) T can be a contact transformation: 
X = T^(x,M(i)), ti^i) = pr(-;^)T"(x, M(i)). The proof is entirely analogous. Similar remarks apply 
to the statements below. 

Note 6. After [60], formula (jl3p and related results were also presented in [111 I15j . The above 
proofs of Proposition [1] essentially differ from the analogous proof in [15] in that the formula for 
transformations of conserved vectors is derived during the proof while in [15] it was proved a 
posteriori. 
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Definition 6. Let G be a symmetry group of the system £. Two conservation laws with the 
conserved vectors F and F' are cahed G-equivalent if there exists a transformation T £ G such 
that the conserved vectors F = and F' are equivalent in the sense of Definition [3l 

Any transformation T € G induces a linear one-to-one mapping Zj, in CV(£), transforms 
trivial conserved vectors only to trivial ones (i.e. CVo(vC) is invariant with respect to Zj.) and 
therefore induces a linear one-to-one mapping 7f in CL(£). It is obvious that 7f preserves linear 
(in)dependence of elements in CL(£) and maps a basis (a set of generators) of CL(£) in a basis 
(a set of generators) of the same space. In this way we can consider the G-equivalence relation 
of conservation laws as well-defined on CL(£) and use it to classify conservation laws. 

Proposition 2. Any point transformation T between systems C and C induces a linear one- 
to-one mapping % from CV(£) onto CV(£), which maps CYo{C) into CVo(-C) and generates a 
linear one-to-one mapping 7f from CL(£) onto CL(£). 

Corollary 8. Any point transformation T between systems C and C induces a linear one-to-one 
mapping from Chf(>C) onto Chf(£). 

It is possible to obtain an explicit formula for the correspondence between characteristics of C 
and C This formula obviously depends on representation of £ and C as systems of differential 
equations. 

Proposition 3. Let T be a point transformation of a system C to a system C and L'^ = A^'^L'^ , 
where A^^ = A^^"D°, A'^'^" are differential functions, and a = (ai,...,a„) runs through the 
multi-indices (ai G N U {0}J, = 1,1. (The number of a's for which A'^'^" ^ is finite.) 
Then the transformation T induces the linear one-to-one mapping from Ch(£) onto Ch(£), the 
inverse of which is defined by the formula 

X^' = A^''*{\D^x\~X''). 

Here A''"* = (--D)" • A'^'^" is the adjoint to the operator A"!^. 

Proof. By the definition of characteristics, for any A G Ch(£) there exists F £ CV(£) such that 
X^^L^^ = D^F^. We take the preimage F E CV(£) of F with respect to the mapping induced 
by T. Then 

DjF^ = \D^x\ DiF' = \D^i\ A^A'^^L^ = A''''*{\D^x\ \^)L'' + DiF' = X'^L'' + DiF\ 

where A'^ = A'^'^*{\Dxx\X'^) and each vanishes on the solutions of C, i.e., the tuple F = 
{F^, . . . ,F^) is a trivial conserved vector of C. It means that the tuple A = (A-*^, . . . , A') is a 
characteristic of the system £, associated with the conservation law containing the preimage 
of the conserved vector F. Since the matrix-operator A = (A^'^) has an inverse which is of a 
similar form, the induced mapping of characteristics is one-to-one. The linearity of this mapping 
is obvious. □ 

Note 7. A^^'^ = for \a\ > in a number of cases, e.g., if C and C are single partial differential 
equations (/ = !). Then the operators A^*^ are simply differential functions (more precisely, the 
operators of multiplication by differential functions) and, therefore, A^^^* = A"^. 

Consider the class C\s of systems Ce'- L{x,U(^p-j,9{x,U(^p-^)) = parameterized with the 
parameter-functions 6 = 9{x,U(^p)). Here L is a tuple of fixed functions of x, n(p) and 6. 9 de- 
notes a tuple of arbitrary (parametric) functions 9{x, U(^p-^) = {9^{x, U(^p^), . . . , 9^{x, U(p))) running 
through the solution set S of the system S{x,Ui^p-^,9(^q-^{x,U(^p^)) = 0. This system consists of dif- 
ferential equations for 9, where x and ti(p) play the role of independent variables and ^(g) stands 
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for the set of all the partial derivatives of 9 of order not greater than q. Sometimes the set 
S is additionally constrained by the non-vanishing condition S"(x, U(p), ^ with 

another tuple S' of differential functions. (See also [56j for other nuances in the rigorous def- 
inition of classes of differential equations.) In what follows we call the functions 6 arbitrary 
elements. Denote the point transformation group preserving the form of the systems from C\s 
as G~ = G~(L,5). 

Let P = P{L,S) denote the set of all pairs each of which consists of a system Cg from C\s 
and a conservation law of this system. In view of Proposition [21 the action of transformations 
from C on C\s and {CV(£5i) \ 9 S} together with the pure equivalence relation of conserved 
vectors naturally generates an equivalence relation on P. 

Definition 7. Let 9,9' £ S, T £ CL{Cg), T' G CL(£e'), F ^ T , F' £ P. The pairs {Ce,T) 
and {Cqi^J-') are called -equivalent if there exists a transformation T £ G~ which transforms 
the system Cg to the system Cq' and such that the conserved vectors F = and F' are 
equivalent in the sense of Definition [31 

The classification of conservation laws with respect to C will be understood as classification 
in P with respect to the above equivalence relation. This problem can be investigated in a way 
similar to group classification in classes of systems of differential equations, especially if it is 
formulated in terms of characteristics. Namely, we construct firstly the conservation laws that 
are defined for all values of the arbitrary elements. (The corresponding conserved vectors may 
depend on the arbitrary elements.) Then we classify, with respect to the equivalence group, the 
arbitrary elements for each of which the system admits additional conservation laws. 

In an analogous way we also can introduce equivalence relations on P which are generated 
either by generalizations of usual equivalence groups or by all admissible point or contact trans- 
formations (also called form-preserving in [H]) in pairs of equations from C\s- 

Note 8. It can easily be shown that all the above equivalences are indeed equivalence relations 
(i.e., are reflexive, symmetric and transitive). 

3.4 Action of symmetry operators on conservation laws 

If the system C admits a one-parameter group of transformations then the infinitesimal generator 
Q = ^"^di + rj'^dua of this group can be used for the construction of new conservation laws from 
known ones. Namely, differentiating equation (I13p with respect to the parameter e and inserting 
the value e = 0, we obtain a new conserved vector expressed via the coefficients of the operator Q 
and the known conserved vector. 

Proposition 4. If Q = C^i + rf'^dua is a Lie symmetry operator of the system C and F G CV(£) 
then the differential functions 

r = -Q^,)F^ + iDje)F= - {Dje)F' (14) 

also are components of a conserved vector of C. Here Q(^.) denotes the r-th prolongation 14 9[ \51 ^ 
of the operator Q. 

Note 9. In contrast to formula (jl3p . formula (jl4p is well-known and extends directly to gener- 
alized symmetries. See, for example, [321 I^Ol [39] for generalized symmetries in the evolutionary 
form (^* = 0) and [17] for the general case. Below we show that in fact it is enough to restrict 
oneself to the version of the formula for the evolutionary form of symmetries. It was used in [JO] 
to introduce a notion of basis of conservation laws as a set which generates a whole space of 
conservation laws by the action of generalized symmetry operators and the operation of linear 
combination. Here we give formula ()14p only through its connection with formula (1131) . 
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There is a well-defined equivalence relation on the space GS(£) of generalized symmetries 
of a system £ of differential equations [50]. Namely, generalized symmetry operators Q and 
Q' of the system C are called equivalent if the difference of their evolutionary forms vanishes 
on the solutions of C. The corresponding factor-space will be denoted by GSf(£). The equiva- 
lence relation on GS(£) agrees with the equivalence relation on CV(£) in view of the following 
statement. 

Proposition 5. The action of equivalent generalized symmetry operators on equivalent con- 
served vector generates equivalent conserved vectors. 

Proof. We show at first that a generalized symmetry operator Q and its evolutionary form Q 
generate equivalent conserved vectors, acting on the same conserved vector. Indeed, since 

Q{oo) = Q{oo) + f A then 

= -Q^r)F' + DjiCF^ - ^F') - CDjF^ = + + F\ 

where are components of the conserved vector F = Q[F] defined by (fT4|) . and Q(oo) and (5(oo) 
denote the formal infinite prolongations of the operators Q and Q, respectively. The differential 
functions = Dj{^'^F^ — F^) are the components of a null divergence and = —S^^DjF^ 
vanish on the solutions of C. Therefore, the conserved vectors F and Q[F] = —Q(^j,-^F are 
equivalent. 

If the difference of generalized symmetry operators in the evolutionary form vanishes on the 
solutions of C, the difference of their actions on a conserved vector obviously have the same 
property. 

The action of any generalized symmetry operator in evolutionary form on a trivial conserved 
vector results in a trivial conserved vector. Indeed, consider an operator Q £ GS(£) in evo- 
lutionary form of order r and a trivial conserved vector F of the system C, i.e., F = F + F, 
where F\j,= and DivF = 0. In view of the Hadamard lemma and the condition F\^= 0, each 
component of F is presented in the form = X^^^U^. Here r is the order of F, X^^ and L^^, 
fi = 1, . . . ,lr, are smooth functions in the jet space J^'''\ the tuple {L^^, . . . , U'^^) determines the 
manifold £(^) in Then Q(^,)F* = L'^^'Q(^,.)X'^' + A*'^Q(^,)L''^ = on where r' = r + r, 

since Q G GS(£). This means that Q[-^]|£= 0. The components of the null divergence F are 
represented in the form F^ = DjV^^ , where v^^ are smooth functions of x and derivatives of u 
such that v^^ = —v^^. The equality 

Q[F] = -Qioo)Djv'^ = -DjQioo)v'' = -D.v'i 

implies that Q[F] also is a null divergence since t"*-^ = Q(oo)V^'' are differential functions and 
= Therefore, Q[F] = Q[F] + Q[F] is a trivial conserved vector as the sum of the tuple 

Q[F] vanishing on solutions of C and the null divergence Oi-^]- 

Let Q,Q' e GS(£), F,F' e CV(£), Q - Q' and F ~ F'. Then Q[F] ~ Q'[F'] since 

Q'[F'] - Q[F] = iQ' - Q')[F'] - (Q - Q)[F] + (Q' - Q)[F'] + Q[F' - F] 

is a trivial conserved vector of the system C in view of the above considerations. □ 

Corollary 9. For any system C formula (fT^ gives a well-defined action of elements from 
GSf(£) on conservation laws of C 

That is why formula (|14p is usually presented only for generalized symmetry operators in the 
evolutionary form = 0) and only inequivalent generalized symmetry operators should be used 
to generate new conservation laws from known ones. Note additionally that the application of 
formula (|14p does not guarantee the construction of nontrivial conserved vectors from nontrivial 
ones [32l|40l[49]. 
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3.5 Potential systems 



If the local conservation laws of a system C of differential equations are known, we can apply 
Lemma [2] to conservation laws constructed on the set of solutions oi 11, = C?. In this way we 
introduce potentials as additional dependent variables. Then we attach the equations connecting 
the potentials with the components of the corresponding conserved vectors Xo C? . If n > 2 the 
attached equations of this kind form an underdetermined system with respect to the potentials. 
Therefore, we can also add gauge conditions on the potentials to . In fact, such additional 
conditions are absolutely necessary in the case n > 2. It was proved in Theorem 2.7 of [3j for 
a quite general situation that every local symmetry of a potential system with potentials which 
are not additionally constrained is projectable to a local symmetry of the initial system, i.e., 
such a potential system gives no nontrivial potential symmetries. Moreover, each conservation 
law of such a system is invariant with respect to gauge transformations of the potentials [6]. 

We have to use linearly independent conservation laws since otherwise the introduced poten- 
tials will be dependent in the following sense: there exists a linear combination of the potential 
tuples, which is, for some r' G N, a tuple of functions of x and it(r') only. 

Then we exclude the superfluous equations (i.e., the equations that are dependent on equa- 
tions from and the attached equations simultaneously) from the extended (potential) sys- 
tem C^, which will be called a potential system of the first level. Any conservation law of is 
one of C^. We iterate the above procedure for to find its conservation laws which are linearly 
independent with those from the previous iteration and will be called potential conservation laws 
of the first level. 

We continue this process as long as possible (i.e., the iteration procedure has to be stopped if 
all the conservation laws of a potential system £,^~^^ of the {k + l)-st level are linearly dependent 
with the ones of C^) or inductively construct infinite chains of conservation laws. This pro- 
cedure may yield purely potential conservation laws of the initial system C, which are linearly 
independent with local conservation laws and depend explicitly on potential variables. The 
idea of this iteration procedure can be traced back to the well-known paper by Wahlquist and 
Estabrook [79] . 

Any conservation law from the previous step of the iteration procedure will be a conservation 
law for the next step. Conservation laws which are obtained on the next step and depend only 
on variables of the previous step are linearly dependent with conservation laws from the previous 
step. It is also obvious that the conservation laws used for the construction of a potential system 
of the next level are trivial on the manifold of this system. 

Since gauge conditions on potentials can be chosen in many different ways, an exhaustive 
realization of the above iteration procedure is improbable if n > 2. 

The best way to calculate conservation laws on each level is to apply the direct method 
of finding conservation laws. One can distinguish four versions of this method depending on 
what condition on conservation laws (([9]), (llOp . (jlll) or ()12p ) is taken as a basis for performing 
calculations [U O [T6l \59\ I80j . Each of these four versions of the direct method has its advan- 
tages and disadvantages in applications and concerning implementation in computer algebra 
programs ^80j . The version involving characteristics is close to the the symmetry group method 
by Noether, which is applicable only in the case of Euler-Lagrange equations and is effective if 
the generalized symmetry algebra of the system under consideration is already known. 

The case of two independent variables is distinguished by the possible (constant) indetermi- 
nacy after the introduction of potentials and also by the high effectiveness of the application of 
potential symmetries. This is why we consider some notions connected with conservation laws 
in this case separately. We denote the independent variables as t (the time variable) and x (the 
space variable). Any local conservation law has the form 




(15) 
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where Dt and are the operators of total differentiation with respect to t and x. The com- 
ponents F and G of the conserved vector [F, G) are called the density and the flux of the 
conservation law, respectively. The conservation law allows us to introduce the new dependent 
(potential) variable v by means of the equations 



determining f up to a constant summand. 

If £ is a single equation, it is a differential consequence of (jl6p in the case of a nonsingular 
characteristic and equations of form (jl6p combine into the complete potential system. As a rule, 
systems of this kind admit a number of nontrivial symmetries and so they are of great interest. 

Introducing a number of potentials for an iteration step in the case of two independent 
variables, we can use the notion of potential dependence, which is more general than the one 
based on linear dependence of conservation laws. 

Definition 8. The potentials v^, . . . , are called locally dependent on the set of solutions of 
the system C (or, briefly speaking, dependent) if there exist r' € N and a function H of the 
variables t, x, U(^^.i\, f-*^, . . . , v'^ such that H^s for some s and li{t^ x, w^, . . . , v^) = for 



any solution {u,v , . . . , v^) of the total system determining the set of potentials v"^ , . . . , v'^ (up 
to gauge transformations, i.e., up to adding negligible constants to the potentials). 

The proof of local dependence or independence of potentials for general classes of differential 
equations is difficult since it is closely connected with a precise description of the possible 
structure of conservation laws. 

Proposition [2] and equation (I16p imply the following statement |60j . 

Proposition 6. Any point transformation connecting two systems C and C of PDEs with two 
independent variables generates a one-to-one mapping between the sets of potential systems cor- 
responding to C and C This mapping is induced by trivial prolongation on the space of the 
introduced potential variables, i.e. one can assume that the potentials are not transformed. 

Corollary 10. The Lie symmetry group of a system C of differential equations generates an 
equivalence group on the set of potential systems corresponding to C 

Corollary 11. Let C\s be the set of all potential systems constructed for systems from the 
class C\s with their conservation laws. The action of transformations from G""{L,S) together 
with the equivalence relation of potentials naturally generates an equivalence relation on C\s. 

Note 10. Proposition[6]and its corollaries imply that the equivalence group for a class of systems 
or the symmetry group for a single system can be prolonged to the potential variables for any 
step of the direct iteration procedure. It is natural that the prolonged equivalence groups and 
symmetry groups are used to classify possible conservation laws, potential systems and potential 
symmetries in each iteration. 

Definition 9. Every Lie symmetry of a potential system is called a potential symmetry of 
the initial system. A potential symmetry operator is called nontrivial if it is not projectable 
to the space of independent and (original) dependent variables, i.e., if some of the coefficients 
corresponding to these variables explicitly depend on potentials. 

The notion of generalized (resp. nonclassical, resp. conditional, resp. approximate, etc.) po- 
tential symmetry is defined in an analogous way. 

Each tuple of independent conservation laws generates, in fact, an infinite series of potential 
systems associated with equivalent tuples of conserved vectors. These systems are connected via 
transformations of potentials of the form v'^ = v'^ -\- H'^, s = 1, . . . ,p, where p is the number of 



= F, vt = -G 



(16) 
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conservation laws in the tuple and the are functions of x and derivatives of u. Therefore, they 
are equivalent in the investigation of generalized symmetries of arbitrary orders. At the same 
time, the choice of representatives in sets of equivalent conserved vectors becomes significant if 
one considers generalized symmetries of a fixed order, e.g., in the case of Lie symmetries having 
order 0. For some classes of equations, a favorable choice can explicitly be given. 

4 Local conservation laws 

We look for (local) conservation laws of equations from class ([T]), applying the modification 
of the direct method which was proposed in [60]. Since equation ([T]) is two-dimensional, the 
constructed conservation laws will have the general form (jlSp . 

At first we prove a lemma on the order of local conservation laws for equations from class ([T]) . 

Lemma 3. Any local conservation law of any equation from class ([T]) is of first order and, 
moreover, it possesses a conserved vector with density depending at most on t, x, and u and flux 
depending at most on t, x, u and Ux- 

Proof. Consider a conservation law (jlSp of an equation of form ([1]). In view of equation ([T]) and 
its differential consequences, we can assume that F and G depend only on t, x and = d^u/dx^, 
k = 0, r', where r' ^ 2r. Suppose that r' > 1. We expand the total derivatives in (jl5p and 
take into account differential consequences of the form utk = F>^{Auxx + Bux + Cm), where 
utj^ = d'^^^u/dtdx^ , k = 0,r'. As a result, we obtain the following condition 

Ft + F^^D^iAuxx + Bux + Cu) + Gx + G^.u^+i = 0. (17) 

Let us decompose pT|) with respect to the highest derivatives Uk- Thus, the coefficients of Ur'+2 
and n^'+i give the equations F^^, = 0, Gu^, + AFu^,__^ = 0, which implies 

F = F, G = -AFu,.,_,Ur' + G, 

where F and G are functions t, x, u, ui, . . . , Ur'-i- Then, after selecting the terms containing 
u^,, we obtain that — = 0. It follows that F = F^u^'-i + where F^ and F° 

depend only on t, x, u, ui, . . . , Uj.i-2- 

Consider the conserved vector with density F = F — DxH and fiux G = G + DtH, where 
H = J F^dur'-2- This conserved vector is equivalent to the initial one, and 

F = F{t, X,U,Ui, . . . , Ur'-2), G = G{t, X,U,Ui, . . . , Ur'-l). 

Iterating the above procedure a suitable number of times, we obtain an equivalent conserved 
vector depending only on t, x, u and Ux, i.e., we can assume at once that r' ^ 1. Then the 
coefficients of Uxxx and Uxx in (HZ]) lead to the equations Fu^ = 0, + AFu = 0, implying 
F = F{t, X, u) and, moreover, G = —AF^Ux + G^, where G^ = G^{t, x,u). □ 

Note 11. A similar statement is true for an arbitrary (l-l-l)-dimensional evolution equation C 
of even order r = 2f, f S N. For example [32], for any conservation law of C we can assume up 
to equivalence of conserved vectors that F and G depend only on t, x and derivatives of u with 
respect to x, and that the maximal order of derivatives in F is less than f. 

The proof of Lemma [3] can easily be extended to other classes of (l-l-l)-dimensional evolution 
equations of even order and some systems connected with evolution equations |60] . 
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Theorem 4. For an arbitrary equation of the form ([T]) the space of all local conservation laws 
is generated by the conserved vectors 

(qu, —aAux + {{aA)x — aB)uj, (18) 

where the characteristic a = a{t, x) runs through the solution set of the adjoint equation 

at + {Aa)xx - {Ba)x + = (19) 

Proof. In view of Lemma [3] we look for (local) conservation laws of equations from class ([T]) in 
the form DtF{t, x, u) + DxG{t, x, u, Ux) = 0. First, we expand the total differentiation operators 
in the latter equality on the solution manifold of ([T]): 

Ft + Fu{Auxx + Bux + Cu) + Gx + GuUx + Gu^Uxx = 0, 

and split the obtained expression with respect to the unconstrained variable Uxx- Coefficients 
of the first power of Uxx give G = —AF^Ux + G^{t,x,u). Splitting the rest of the conservation 
law with respect to different powers of Ux yields 

Fun = 0, F^B - AxF^ - AFxu + = 0, Ft + CuF^ + = 0. 

Solving the obtained system up to the usual equivalence relation of conserved vectors, we obtain 
the conserved densities and fiuxes of the local conservation laws of ([T]) : 

F = au, G = —aAux + {{aA)x — aB)u, 

where a = a{t,x) is an arbitrary solution of equation (|19p . □ 

Note 12. It is well-known that the space of characteristics of the local conservation laws of a 
linear differential equation Lu = contains all functions of the independent variables which are 
solutions of the adjoint equation L*a = 0. In view of Theorem H] any local conservation law of a 
(1 + l)-dimensional linear second-order parabolic equation has such characteristics. Moreover, 
two different solutions of the adjoint equation give inequivalent characteristics corresponding to 
inequivalent conservation laws. Therefore, for any equation C from class ([T|) we can identify the 
solution space of the adjoint equation ()19p with the space Chf(£) of 'nontrivial' characteristics. 

Note 13. The above structure of the space Chf(£) is directly related to the restriction on the 
order of conservation laws for evolution equations of even orders. Third-order (l + l)-dimensional 
linear parabolic equations admit conservation laws with different characteristics. For example, 
for the equation ut = Uxxx the space of conservation laws of orders not greater than 3 is generated 
by the conservation laws with the conserved vectors {au, ), (m^,m^^ — 2uuxx) and 

{ux'^,Uxx — "^UxUxxx) and the characteristics a, 2u and 2ux, respectively. Here a = a{t,x) runs 
through the solution set of the same equation at = axxx- The conservation laws parameterized 
by a have the structure usual for linear equations but the two other conservation laws do not. 
Moreover, acting by the symmetry operators dx + Cudu and 3tdt + xdx + Gudu on solutions and 
conservation laws of the equation ut = Uxxx, we construct infinite series of conservation laws for 
it which are of arbitrarily large orders with conserved vectors quadratic in the derivatives of u. 

Let us emphasize once more that the function a = a{t, x) is a characteristic of the conservation 
law of ([I]) with the conserved vector ()18p and, therefore, an adjoint symmetry of equation ([1]). 
Equation (jl9p is the formally adjoint equation to ([1]). Any solution of (jl9p is an adjoint symmetry 
of and any adjoint symmetry of ([I]) depends only on t and x and is a solution of (fT9]l . 
Conservation laws of ([T]) are linearly independent iff the corresponding solutions of (jl9p are 
linearly independent. 
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Corollary 12. There is a one-to-one correspondence between local conservation laws and zero- 
order adjoint symmetries of equation ([T|). 

Proposition [H Corollary [3] and Theorem S] imply the following statement. 

Proposition 7. Any point transformation between equations from class ([T|) is canonically pro- 
longed to the characteristics of conservation laws of equations from this class by the formula 

(The constant k arises due to the linearity of the characteristic space. It is inessential and can 
be set equal to 1.) 

Note 14. Any equation ([1]) is reduced by an equivalence transformation to an equation of the 
form ([7]). In view of Proposition [2] this allows us to restrict ourselves to an investigation of 
conservation laws for the simpler reduced form d?]) of parabolic equations. The space of local 
conservation laws of an equation ([7]) is generated by the conserved vectors {au, —aux + Oxu), 
where the characteristic a = a{t,x) runs through the solution set of the associated adjoint 
equation at + axx — = 0. At the same time, it is by no means evident that for a further 
study of potential conservation laws and potential symmetries it suffices to consider only the 
reduced form. This point deserves additional investigation. 

The equivalence relation generated by the equivalence group on the set of pairs ' (equation from 
the class, its conservation law)' can be used for the normalization of parameters in different ways. 
We can maximally simplify the form of the equations under consideration as in Notell4[ Another 
way is to simplify equations and their conserved vectors or characteristics simultaneously. For 
example, let us act on equation ([T|) and its characteristic by the equivalence transformation 
from G~ with Tj = signA, Xx = \A\~^^'^, = q;|A|^/^. Then A = 1, a = 1 and, moreover, 
C = Bx since a should be a solution of the transformed adjoint equation {Aa)xx — {Ba)x + 
Ca = 0. As a result, the following proposition holds: 

Proposition 8. Any pair {C,J-), where C is an equation from class ^ and T G CL(£), 
is reduced by a point transformation from the equivalence group of class ^ to a pair {C,J-), 
where C is a Fokker-Planck equation ut = Uxx + {Bu)x and T is its conservation law with the 
characteristic a = 1, i.e., (n, —Ux — Bu) £ T . 



5 The adjoint variational principle 

The adjoint variational principle is well-known especially in physics. In the nonlinear case it is 
also called the composite variational principle [7] . It provides a way of constructing generalized 
Lagrangians for systems of differential or other equations which have no usual Lagrangians, e.g., 
for evolution equations. Thus, the adjoint variational principle for the linear heat equation was 
considered in the classical textbook by Morse and Feshbach [481 p. 313]. It was discussed in [66] 
for general linear operator equations. A comparison of the adjoint variational principle with 
other approaches to deriving generalized Lagrangians was presented in [751 pp. 341-342] along 
with a review of the literature on its application. A nonlinear version was proposed in p6j and 
used ibid, for the Navier-Stokes equations. The composite variational principle in its general 
formulation was given in [7j. See also )49i Exercise 5.27] and the discussion in section 2 of |17j . 

Briefly the composite variational principle can be described in the following way. Let L be 
a system 'U(p)) = of / differential equations = 0, . . . , L' = for m unknown functions 
u = {u^, . . . , n™) of n independent variables x = (xi, . . . , Xn)- (Here we employ the notations of 
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Section [3l Thus, the index // runs from 1 to I.) We introduce I auxihary dependent variables 
V = {v^, . . . ,v^) and construct the Lagrangian 

2, = vL = 

The Euler-Lagrange equations for the corresponding functional are 
= 0, Dl{v) = 0, 

where the matrix differential operator is the adjoint of the Frechet derivative (see subsec- 
tion [32])- The system D*i^{v) = will be called variationally adjoint to the system C. Only 
in the linear case the variationally adjoint equations coincide with the usual adjoint equa- 
tions L*v = [7]. Therefore, simultaneously extending the tuple of dependent variables with / 
auxiliary variables and the system with the variationally adjoint equations always results in the 
system formed by the Euler-Lagrange equations of a special functional. In physical terms this 
means that a dissipative system possessing a usual friction is considered simultaneously with 
its 'mirror reflection' possessing a negative friction and absorbing the energy lost by the initial 
system [ISl Section 3.2]. As a result, the total energy of the extended systems is conserved. 
The trick is quite artificial from the physical point of view but nevertheless it has a number of 
applications. It allows one to operate with dissipative systems as if they were conservative. 

Interest in applying the adjoint variational principle in the framework of symmetry analysis 
has arisen quite recently [3H [351 IMl EZ]. Thus, a nice result on prolongation of symmetries of 
initial systems to the auxiliary variables of the adjoint variational principle was proved in |34j . 
The prolongation is carried out in such a way that the prolonged operators are variational sym- 
metries of extended Lagrangians and, therefore, Lie invariance operators of extended systems. 
At the same time, this result can be obtained for linear equations without the usage of the La- 
grangian technique. It can be strengthened for special classes of differential equations including 
class ([H). 

Let £ be a homogeneous linear differential equation for one unknown function li, i.e., m = 1 
and let C be presented in the form Lu = 0. Here L is the associated linear differential operator: 
L = A'^{x)d^ , where v = (z^i, . . . , is a multiindex, f j G N U {0}, \v\ := vi + ■ ■ ■ + Vn ^ P, 
d := {dx^i ■ ■ ■ ,dx,J, d'^ := d'^\ . . .5^^, A'^ = A'^{x) are smooth functions of x, and summation 
over the multiindex ly is understood. The adjoint K* to a linear differential operator K (over 
the real numbers) can be defined by means of the following formal rules. The adjoint to a linear 
combination of operators is the linear combination of the adjoint operators with the same coeffi- 
cients. {K1K2)* = K2KI, {dxi)* = —dxi and the adjoint operator to the multiplication operator 
by a fixed function is the multiplication operator by the same function. Thus, L* = (—dYA'^. 
The equation £*: L*v = is adjoint to C in both the usual linear and the variational sense. 
(According to the adjoint variational principle, only one auxiliary variable v should be intro- 
duced for the equation £,.) Since (L*)* = L, the adjoint equation to C* also can be assumed to 
coincide with the initial equation, i.e., (£*)* = C. The united system C H £*: Lu = 0, L*v = 
is the system of Euler-Lagrange equations of the Lagrangian vLu as well as of the equivalent 
Lagrangian uL*v and any linear combination of these Lagrangians. The Lagrangians vLu and 
uL*v are equivalent since vLu — uL*v is a total divergence, i.e., vLu — uL*v = DivF for some 
n-tuple F of differential functions bilinear in u, v and their derivatives. This formula also implies 
that V S Ch(£) if L*v = 0. The tuple F for a fixed solution v of the equation L*v = gives a 
conserved vector of the conservation law corresponding to this characteristic. 

Suppose that Q = ^^{x)dxi + V^{x)udu is an essential Lie symmetry operator of C and 
Q = —^^{x)dxi + r]^{x) is the associated first-order differential operator acting on functions of x. 
Employing the Hadamard lemma, we write the infinitesimal criterion Q(^p-jLu\Lu=o = in the 
form Q(^p-^Lu = XLu [49]. Here (^(p) is the standard p-th order prolongation of Q and in the case 
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under consideration A is a smooth function only of x. The infinitesimal criterion is equivalent to 
the operator equality [Q, L] = XL, where \ = rj^ — X and [Q, L\ = QL — LQ is the commutator 
of the operators Q and L. 

Proposition 9. Q = f (9^, + r]^udu G 0^''(/:) iff = Cdx, + 0\dy G 0"'^"(£*), w;/iere 

Moreover, Q := f 5^;. + ri^udu + O^vd^ £ 0'="'^(£ n £*) and Q is a variational symmetry operator 
of the Lagrangian vLu. 

Proof. Conjugating the operator equality [Q,L] = XL, we obtain [—Q* + X,L*] = XL*. At the 
same time, —Q* + A = —Cdxi + 6^ is the first-order differential operator acting on functions 
of X which is associated with . Therefore, G q'^^^{C*). Since (L*)* = L, the converse 
statement is true as well. The operators Q and have the same x-part and the system £ n 
£* is uncoupled. Hence Q G q'^^^{C H £*). This also follows from the fact that C and C* 
are the Euler-Lagrange equations of the Lagrangian vLu and Q is a variational symmetry 
operator thereof in view of the infinitesimal criterion of variational symmetry [49j. Indeed, 
Q(^) {vLu) + (vLu) = (vLu) (^i + A + ^ J = 0. □ 

The algebra 0°°(£n£*) is formed by the operators fdu + gd^, where the parameter- functions 
/ = f{t, x) and g = g{t, x) run through the solution sets of the equations C and C* , respectively. 
Any such operator is a variational symmetry operator of the Lagrangian vLu. It is projectable 
to {t,x,u) and {t,x,v) and its projections belong to Q°°{C) and q°°{C*). It can be assumed 
that 0°°(£ n £*) = 0°°(/:) e q'^{C*). The operators from 0~(£) and 0°°(£*) have a trivial zero 
complement to operators from 0°°(>Cn>C*) according to [34J, i-e. the algebras 0°°(£) and q°°{C*) 
cannot be obtained from each other via the adjoint variational principle. A slightly different 
situation obtains for the operators udu and vd^. An arbitrary linear combination of them belongs 
to q{C n C*) but only the operators proportional to udu — vdv are variational symmetries. At 
the same time, among such linear combinations only the operators from {udu + vdu) are trivial 
Lie invariance operators of the system £ n £*. So 0'"''(£ n £*) C q^™{C) Q^"^{C). 

Let us return to linear (1 -|- l)-dimensional second-order parabolic equations. Let C be an 
equation from class ([I]). Then the adjoint equation C* is of the form (fT9]) . i.e., the system £n£* is 

ut = Au^^ + Bu^ + Cu, vt + {Av)xx - {Bv)cc + Cv = 0. (20) 

where A = A{t,x), B = B{t,x) and C = C{t,x) are arbitrary smooth functions, A ^ 0. For 
class ([T|) we extend the adjoint variational principle to admissible transformations. At first 
we present auxiliary statements on admissible transformations for wider classes of systems of 
evolutionary equations. Their proofs are based on the direct method. 

Consider the class of systems of second-order evolutionary equations of the general form 

Ut = F{t,x,u^,u^x), (21) 

where u = {u^, . . . , u™), F = {F^, . . . , F™) and \dF /duxx\ 7^ 0. Any point transformation T in 
the space of variables {t,x,u) has the form 

i=T\t,x,u), X = T''{t,x,u), u"" =W{t,x,u), 

where the Jacobian \d{T* ,Ll) / d{t,x,u)\ does not vanish, U = {U^, . . . ,U^). In what follows 
the indices a, h and c run from 1 to m. 
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Lemma 4. A point transformation T connects two systems from class (j2ip iff = T^a = 0, 
i.e., T* = T(t), where T is an arbitrary smooth function oft such that Tt ^ 0. The arbitrary 
elements are transformed by the formula 

TtDxT"" TtD^r^' 

where Dx = dx + Ux^ub + • • • is the operator of total differentiation with respect to x. Therefore, 
class (|21|) is normalized. The equivalence group of class (|21|) is formed by the transformations 
determined in the space of variables and arbitrary elements by the above formulas. 

Analogously, consider the subclass of class (|21|) formed by systems whose arbitrary elements 
are linear in Uxx, i-e., 

= S''''{t, X, u, Ux)ulx + H'^it, X, u, Ux), (22) 
where \S\ / 0, S = {S"^), H = {H\...,H"'). 

Lemma 5. A point transformation T connects two systems from class ([22|) iff = T^a = 
and T.^a = 0, i.e., T* = Tit) and = X{t,x), where T and X are arbitrary smooth functions 
of their arguments such that TfXx ^ 0. Moreover, the Jacobian \dU/du\ ^ 0. The arbitrary 
elements are transformed by the formulas 

5. f*5 (*)"'. 

Tt ou \ou J 

n = L(^H+Ut- ^DxU -^s(^Y ( DxUx + W^DxU^a - ^Dxd] \ . 
Ttydu Xx du \du J \ X X u Xx J J 

Therefore, class (j22p is normalized. The equivalence group of class (j22p is formed by the trans- 
formations determined in the space of variables and arbitrary elements by the above formulas. 

Note 15. If 5" = diag(5]^, . . . , S"*) and S = diag(S^, . . . , S™") are diagonal matrix-functions 
then in view of Lemma [5] there exists a permutation C7 G such that = X.^Tf^J:''^''^ and 
U";,{T,°'^"'^ — S'') = 0. (There is no summation in the last formula.) This implies in the case 
E» / for any a / 6 that U^^ = if a{a) / b. 

Corollary 13. The subclass of class (j22p . formed by the systems linear in the derivatives, 
i.e., defined by the constraints S'^c = 0, H^c^t = on the arbitrary elements, is normalized. 
The equivalence group of this subclass is a subgroup of the equivalence group of class (j22p and 
is formed by the transformations in which additionally the parameter-functions lA"" satisfy the 
conditions U^,^, = 0, i.e., W = U''^{t,x)u^ + ^^{t^x). 

Corollary 14. The subclass of class (j22p . formed by the linear systems i.e., defined by the 
constraints S^c = S'^c = 0, H'^c^b^ = H'^c^b = H^c^b = on arbitrary elements, is normalized. 
The equivalence group of this subclass coincides with the group described in Corollary \13i 

Combining Note [T] and Corollary [14] implies the following statement on properties of the 
corresponding class of homogeneous linear systems. 

Corollary 15. The subclass of class (I22p formed by the homogeneous linear systems (i.e., the 
constraints on arbitrary elements are S!^c = S^c = 0, -f^^c^t = H'^c^p = H^c^b = and UxH^c + 
u^H^c = H""), is semi-normalized. The equivalence group of this subclass is a subgroup of 
the group from Corollary \13l which is formed by the transformations with U""^ = 0. A point 
transformation T connects two systems from this class iff = T^a = 0, T^a = and h(\ c = 0, 
i.e., r* = T{t), = X{t,x) and W = U''^{t,x)u^ + U'^^H^x), where T,X, U"^ ""and U"^ 
are arbitrary smooth functions of their arguments such that TtXx\U°'^\ ^ 0, and additionally 
^jjabjjW-^ is a solution on the initial system. Here [U"'^) is the inverse matrix of {U""^). 
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Note 16. Representing systems from class (j2ip whose arbitrary elements are linear in Uxx in 
the form (j22p . we replace each F° by the tuple of arbitrary elements {S°'^,H°'). In the situation 
under consideration the different representations give equivalent results. Thus, the formulas for 
the transformation of S""^ and H"^ are obtained from the analogous formulas for F"" by splitting 
with respect to Uxx- The corresponding sets of admissible transformations as well as equivalence 
groups are isomorphic. We can work with the above subclasses of class (|22p in the same way. 
Namely, we can constrain the functions S""^ and successively replace i?" by the expressions 
H'^^t, X, u)ul + G<^(i, X, n), H''\t, + G''\t, x)v}' + C^it, x) and H"-\t, + C^^t, x)vf' . 
The transformation formulas for the new arbitrary elements are constructed by splitting the 
transformation formula for H with respect to Ux or [ux^u). The equivalence groups of the 
different representations of these classes are isomorphic. 

Let us continue with class ()20p . At first we study the corresponding class of inhomogeneous 
systems, writing them in the form 

ut = eiAuxx + B^Ux + C^u + D^, vt = e2Avxx + B'^Vx + C^v + D"^ . (23) 

Hereafter m = 2, ei = 1, £2 = -1, = B, = C, B'^ = B - 2Ax, = -C + Bx - Axx, 
L"" = dt — SaAxx — B°^dx — C". All the arbitrary elements are smooth functions of t and x. 

Corollary [14] implies that any point transformation between two arbitrary systems from 
class (|20p has the form 

i=T{t), x = X{t,x), 

u = U^^{t, x)u + U^^{t, x)v + U^°{t, x), (24) 
v = U^^{t,x)u + U^^{t,x)v + U^°{t,x), 



where TtXx\U''''\ / 0. In view of Note [l5] we additionally have that A = SX^Tf^A, 6 = ±1. 
Moreover, C/^^ ^ jj2i ^ q -^^ ^ ^ ^ ^jii = f/22 ^ q ^ ^ _i ^^y.^ (61,62) = (1,2) if 

6 = 1 and (61, 62) = (2, 1) otherwise. Using these notations the transformation of the arbitrary 
elements can be written as 



Tt '' a rp^ \ XX jjaba ^ j ' rp^ 

Tjaba 1 Tjaba f TTaO 

= L^" — ^, D" = D''- + L''" — 



(25) 



Jjaba ' Tt \ 

The specific connections between the coefficients {B^, C^) and {B^, C^) imply one more equation 

XxU^'^U'^^^ = K, (26) 

where k is an arbitrary nonzero constant. (Compare with the formula from Proposition [71) 

Analyzing the obtained result, we deduce the following: Any admissible transformation in 
class ()23p is either the prolongation of a transformation between the first (or second) equa- 
tions of the related systems according to condition (j26p or the composition of such a prolonged 
transformation with the transformation given by the simultaneous transposition of the depen- 
dent variables and equations in the resulting system {u = v, v = u, A = —A, B = —A, 
C = —C + Bx— Axx)- In the second case connections between arbitrary elements involve 
their derivatives. That is why we should consider the extended equivalence group of class psp . 
admitting the dependence of transformations of arbitrary elements on their (and only their) 
derivatives. 
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Proposition 10. Class (j23p is normalized with respect to the extended equivalence group 
which is formed by the transformations described by the formulas (|24p , (|25p and (j26p . The usual 
equivalence group G~ of class (|23p is ^/le subgroup of G^^^, defined by the condition 5 = 1. Its 
'essential' part formed by the transformation with U""^ = can be obtained as the prolongation 
of the 'essential ' part of the equivalence group of the class ^ (or the class of the same equations 
written in the adjoint form ) according to condition (j26p . Conversely, the equivalence groups of 
the class ([2]) and the class of the same equations written in the adjoint form are the projection 
of G~ to the corresponding sets of variables and arbitrary elements. The complement of C 
in G^^ is formed by the compositions of elements from G~ and the simultaneous transposition 
of the dependent variables and equations. 

Similarly to the above classes of linear evolution systems, in view of Note [1] and Proposition [TUl 
we obtain that the corresponding class ()20p of homogeneous systems is semi-normalized in the 
extended sense. The extended/usual equivalence group of ([20]) is the projection (neglecting 
transformations for and D^) of the subgroup of the extended/usual equivalence group of (123]) . 
defined by the constraints W^^ = 0. The usual equivalence group of (I20p can be constructed 
via the prolongation of the equivalence group of class ([T|) (or the class of the same equations 
written in the adjoint form) according to condition ()26p . Conversely, the equivalence groups 
of ([I]) and the class of the same equations written in the adjoint form are the projections of the 
usual equivalence group of (I20p to the corresponding sets of variables and arbitrary elements. 
The additional condition for the admissible transformations in class ()20p is L^°{U'^^ /If^^") = 0, 
where there is no summation with respect to a. 

In the case of class ([T]) Proposition [9] can be reformulated in a more precise way since any 
essential symmetry operator Q = rdt + £,dx + rj^udu of an equation from class ([T]) has X = Tt- 
The same statement is derived as an infinitesimal consequence of formula (j26p after taking 
into account that any symmetry transformation of any equation from class ([1]) is an admissible 
transformation in this class. 

Corollary 16. Let C be an equation from class ([T]). Then Q = rdt + ^dx + r]^udu G Q^^^{C) iff 
Qt = rdt + Cdx + O^vd^ £ g*'''}/:*), where rj^ + 6^ = -^x- Moreover, 

Q = Tdt + idx + r]\du + e^vd^ g q'''\c n C) 

and Q is a variational symmetry operator of the Lagrangian vL^u (or uL'^v). 

Using equivalence transformations, we can gauge the arbitrary elements A, B and G in 
different ways analogous to gauging them in class ([T]). 

The gauge A = 1 preserves all the normalization properties, relations between equivalence 
groups etc. in both the inhomogeneous and homogeneous cases. Under the gauge ^ = 1 we obtain 
the condition 6X'^ = Tt. Therefore, the constant 6 becomes coupled with the sign of Tt (i.e., 
6 = signTf) and X = 6'\Tt\^^'^x + C(Oi where 6' = ±1. Here is an arbitrary smooth function 
of t. The discrete extended equivalence transformation transposing the dependent variables and 
equations is replaced by its composition with the discrete transformation of alternating the sign 
of t, which also is an extended equivalence transformation. 

The further gauging of i? to gives the conditions 

E^1^JJ^^, ,e.. f/... expf-..^.^-£^-4,.V 

C/'^*" 2 Xx w ag^^ 25' |rt|V2 J 

Here 9^ and 6^ are arbitrary smooth nonvanishing functions of t. A new nuance is that under 
the gauge {A,B) = (1,0) the equivalence transformations involving the transposition of the 
dependent variables and equations reduce to usual ones since they contain no derivatives of 
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arbitrary elements. In view of this, the subclass of class (|23p with the gauge {A,B) = (1,0) 
is normalized in the usual sense. The corresponding subclass of homogeneous systems (|20p is 
semi-normalized in the usual sense. To avoid the transposition transformations in this case, 
statements on relations between equivalence groups should be formulated in other terms, e.g., 
for continuous equivalence groups. 

Note that point symmetry groups of some systems from class (|20p / (j23p contain transforma- 
tions involving the transposition of u and v. Such groups cannot be projected to the symmetry 
groups of single equations in contrast to the Lie symmetry groups. For example, any system 
(j20p with A = 1, B = and Ct = is invariant with respect to the transformation t = —t, 

X = X, U = V, V = u. 

The gauge C = determines the subclasses of the systems of Kolmogorov and Fokker- 
Planck equations associated with each other. The reduced form of such systems is given by 
the gauge {A,C) = (1,0). Both the gauges break the normalization properties. Nevertheless, 
the admissible transformations in these subclasses can be described. They should satisfy the 
conditions L^'^(l/U"'^'^) = 0. The 'essential' admissible transformations (for which = 0) are 
obtained as the prolongations of the 'essential' admissible transformations of the Kolmogorov 
(or Fokker-Planck) equations according to condition (j26p . Conversely, the 'essential' admissible 
transformations in the class of the Kolmogorov (or Fokker-Planck) equations are the projections 
of the 'essential' admissible transformations of the associated systems, which do not involve the 
transposition of the dependent variables and equations, to the corresponding sets of variables 
and arbitrary elements. Therefore, the sets of the 'essential' admissible transformations of the 
Kolmogorov and Fokker-Planck equations are similar. This implies the following statement in 
view of Corollary [7] and Proposition [9j 

Consider the class of Fokker-Planck equations having the form 

vt = {A{t, x)v^)^ - {B{t, x)v)^. (27) 
(For the form to be canonical, we alternate the sign of t.) 

Corollary 17. The kernel Lie algebra of class (|27p is {vd^)- Any equation from this class 
is invariant with respect to the operators fdy, where the parameter-function f = f{t,x) runs 
through the solution set of this equation. Any equation from class (j27p is reduced by a point 
transformation to an equation with A = 1 from the same class. Up to point transformations, all 
possible cases of extension of the maximal Lie invariance algebras in class ()27p are exhausted 
by the following ones (in all the cases A = 1; the values of B are given together with the 
corresponding maximal Lie invariance algebras): 

1. B = B{x): {dt, vd,, fd,); 

2. B = ux-'^, V ^ 1, u j^2: {dt, D, U + 2utvd^, vd^, fd^); 

3. S = x~i(l-2xtan(xln|j;|)), x/0: {dt, 2D - {xB -2)vdy, U + 2txBvdy, vd^, fd^); 

4. B = 0: {dt,d,,G,D,U,vd,,fd,). 

Here D = 2tdt + xd^, U = At^dt + Uxd-, - (x^ + 2t)vd^, G = 2td^ - xvd^. 

Note that the kernel Lie algebra of (j27p is different from the one of the subclass of ([1]) with 
C = since the prolongation of the operators to v depends on the arbitrary elements. 

Analogously to the 'Kolmogorov' form, the group classification of class (j27p with respect to 
its equivalence group can be derived from the classification presented in Corollary [17] through 
extending the classification cases by essential admissible transformations which are not generated 
by the equivalence group. 
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6 Potential conservation laws 



As proved in TheoremlU any (l+l)-dimensional linear second-order parabolic equation possesses 
an infinite series of local conservation laws. Studying potential conservation laws of equations 
from class ([T]) , in view of Note [14] and Proposition [6] we may restrict ourselves to equations of 
the reduced form ([7]), i.e., ut — Uxx + = 0. Fixing an arbitrary p € N and choosing p linearly 
independent solutions a = (a^, . . . , a^) of the associated adjoint equation 



we obtain p linearly independent conservation laws of equation ([7]) with the conserved vectors 
{F^,G^) = {a^u,a^u — a'^Ux)- (Hereafter the indices s, a and <j run from 1 to p. Let us also 
recall that summation over repeated indices is assumed.) 

The potentials v = {v^ , . . . ,v^) introduced with these conservation laws by the formulas 



are independent in the sense of Definition [8] according to the following lemma. 

Lemma 6. For any equation ([7]) the potentials are locally dependent on the equation manifold iff 
the corresponding conservation laws and, therefore, the corresponding solutions of equation ([28ll 
are linearly dependent. 

Proof. The 'if part of the claim being obvious, we only prove the converse statement, arguing by 
contradiction. Suppose that the potentials f^, . . . , f ^ introduced with the linearly independent 
solutions a^, . . . , of (j28p are locally dependent. In case p = 1 we have local triviality of 
as a potential, i.e., f ^ can be expressed in terms of local variables and, hence, the corresponding 
conservation law is trivial. Therefore it suffices to investigate the case where the number of 
independent conservation laws is greater than 1. 

Without loss of generality we may assume that there exist r G N and a fixed function P 
of t, X, V = (w^, . . .v^) and U(^r) such that = P{t,x,v,U{j.)) for any solution of the combined 
system determining the whole set of potentials f^, (up to gauge transformations, i.e., 

up to adding negligible constants to the potentials). In view of equation ([7]) and its differential 
consequences, we may suppose that P depends only on t, x, v and Uk = d^u/dx^, k = 0, r', 
where r' ^ 2r. Let us apply the operator Dx to the condition = P{t,x,v,u,ui, . . . ,Uri): 

= Px + Pjjs'Vx + Pu^Uk+i- (The index s' runs from 2 to p.) Taking into account the equations 
Vx = a^u, we split the differentiated condition with respect to step-by-step in reverse order, 
beginning with the highest derivative. As a result, we obtain = 0, Px = and = P^s'Oi'^ , 
i.e., the functions a^, . . . , are linearly dependent over the ring of smooth functions of t. In 
view of Corollary [18] (see below), this means that the functions a^, . . . , are linearly dependent 
in the usual sense, contradicting the assumed independence of the conservation laws. □ 

Let W {(f^ , . . . , (f'') denote the Wronskian of the functions ip^, with respect to the 

variable x, i.e. W{ip^, ...,</?') = det{ipl_-^)l.j^i. 

Lemma 7. The solutions (p^ = (p'^{t,x), y?' = ip)\t,x) of a {1 + 1)- dimensional linear 
evolution equation L(p = of arbitrary order are linearly dependent iff W{ip^, . . . , if'') = 0. 

Proof. Since the equation Lip = is linear and evolutionary, the operator L is the sum of dt 
and a linear differential operator with respect to x whose coefficients depend on t and x. If the 
functions p^, . . . , p>'' are linearly dependent then the equality W{p^, . . . ,p^) = is obvious. Let 
us prove the converse statement. 



(28) 



v„ = a u. 



vl = a^Ux — c^xU 



(29) 
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In the case I = 2 the condition W{(p^,ip'^) = imphes = C(/9^, where C is a smooth 
function of t. Acting on the latter equahty with the operator L, we obtain Ct^p^ = 0, i.e. 
C = const or 93^ = 0. In any case the functions ip^ and ip'^ are hnearly dependent. 

Suppose W{ip^, . . . , if'') = 0. Without loss of generality we can assume W{ip^, . . . , if''^^) / 0. 
(Otherwise we consider a smaller value of I.) Then ip'' = C^ip^, where are smooth functions 
of t and the superscript k runs from 1 to / — 1. Acting on the latter equality with the operator L 
results in the equation C^ip^ = which implies, in view of the condition W{ip^ , . . . , ^^~^) 7^ 0, 
= const, which concludes the proof. □ 

Corollary 18. Solutions of a [1 + 1) -dimensional linear evolution equation of an arbitrary order 
are linearly dependent in the usual sense iff they are linearly dependent over the ring of smooth 
functions oft. 

For equation ([7|) the complete set of first level potential conservation laws is indeed the union 
of conservation laws of systems (I29p corresponding to all possible values of p and p-tuples a. 



Theorem 5. Any local conserved vector of system (j29|) is equivalent on the manifold of sys- 
tem (j29|) to a local conserved vector of equation ^ . 

Corollary 19. For any linear (1 + 1)- dimensional second-order parabolic equation potential 
conserved vectors of any level are equivalent to local ones on the manifolds of the corresponding 
potential systems, and potentials of any level can locally be expressed via local variables t, x, u^^) 
(for some r ) and potentials of the first level only. 

Corollary 20. For any linear (1 + 1)- dimensional second-order parabolic equation any potential 
system of higher level is equivalent, with respect to point transformations nontrivially acting only 
on potential variables, to a potential system of the first level. 

In other words, the set of potential conservation laws of any linear (l + l)-dimensional second- 
order parabolic equation is exhausted by its local conservation laws. The set of locally indepen- 
dent potentials is exhausted by potentials of the first level. 

Following pO] where these statements were derived for the linear heat equation ui — Uxxt 
we present the proof of Theorem [5] in the form of a sequence of basic lemmas. In comparison 
with [60], only minor modifications of the proof are needed due to the reduction of equations 
from class ([T]) to the simpler from ([7]) by equivalence transformations. 

Lemma 8. Any local conservation law of system (j29p is equivalent to the one with the conserved 
vector {Ku,KxU — Kux) where the function K = K(t,x,v) is determined by the system 

Kt + Kxx -VK = Q, a'Kxvs - a^K^s = 0. (30) 



Proof. Consider a local conservation law of system (|29p in the most general form, where the 
conserved vector is a vector- function of t, x and derivatives of the functions u and v'^ from 
order zero up to some finite number. Taking into account system ()29p and its differential 
consequences, we can exclude the dependence of the conserved vector on any of the (nonzero 
order) derivatives of v" and the derivatives of u containing differentiation with respect to t. 
Similarly to Lemma[3]we can prove that the reduced conserved vector (F, G) does not depend on 
(nonzero order) derivatives of u and, moreover, F = F{t, x,v), G = —a^Fys (t, x, v)u+G^{t, x, v). 
The functions F and G^ satisfy the system 



a'a^F^s^. = 0, a'G% = 2a^F„. + a'Fxv^ , Ft + G° = 0. 



Let us pass on to the equivalent conserved vector {F, G), where F = F -\- D^H, G = G — DfH 
and H = H{t, x, v) is a solution of the equations = —F, Ht = G. (The variables v'^ are 
considered as parameters in the latter equations.) Then F = Ku, G = KxU— Kux- The function 
K = a^Hys depends on t, x and v and satisfies system (pOl) . □ 
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Lemma 9. Let the solutions a* = a^{t,x) and (3^ = P^{t,x) of equation (j28p satisfy the addi- 
tional condition — q*/3|. = 0. Then for any i, j G NU {0} 

a|/?| - a|/?| = 0. (31) 

Here and it what follows the subscripts i and j denote the i-th and j-th order derivatives with 
respect to x. 

Proof. We carry out an induction with respect to the value i + j. 

Equation ()3ip is trivial for i + j = 0, coincides with the additional condition for i + j = 1 
and is obtained from this condition by means of differentiation with respect to x if z + j = 2. 

Let us suppose that the Lemma's statement is true if i + j = m — 1 and i + j = m and prove 
it for i + j = m + 1. Acting on equation (I3ip where i + j = m — 1 with the operator dt + dxx 
and taking into account the conditions af + — Va'^ = and /3| + j3^^ — V (3'^ = 0, we obtain 
the equation a|_,_]^/3j^;^ — a^^j^-^fifj^i = 0. Therefore, the statement is true for i' + j' = m + 1, 
1 ^ i' ,j' ^ m (here i' = i + 1, j' = j + 1). It remains to prove the case i' = m + 1, j' = 
(or equivalently i' = 0, j' = m + 1). For these values of i' and j' the statement is obtained by 
subtracting the above equation a^/3f — af = from the results of differentiating equation ([3T]) 
where i = m, j = with respect to x. □ 

Lemma 10. // af/3j - a^f = /or ^ ? < j ^ p then W{a^, . . . ,aP, (3"") = for any a. 

Proof. Let M^^ denote the {p — l)-st order minor of W{a, (3"), which is obtained by deleting 
the (T-th and {p + l)-st columns corresponding to the functions a" and (3"^ and the i-th and 
j-th rows, where ^ i < j ^ p. (For convenience the rows of W{a, (3^) are enumerated from 
to p.) We multiply the equation a^/3| — a*/?| = by (— +'^+p+^Mj!^- and convolve with 
respect to the index pair In view of the Laplace theorem on determinant expansion we 

obtain W{a\a^.^a'^.,l3'^) = 0. Here the symbol means that the function a'^ is substituted 
instead of the function a"^, and summation over the index s is carried out. This latter equation 
immediately implies the lemma's statement since W {a\a^ .^a" , (3^) = for any fixed s ^ a. □ 

Lemma 11. The general solution of system (I30p can he written in the form 

K = a'H^s + (32) 

where H is an arbitrary smooth function of v, and (3^ = /3^{t,x) is an arbitrary solution of 
equation (i28l) . 



Proof. In view of Lemma [8] the functions a'^ and = K^s satisfy the conditions of Lemma [9] 
and, therefore, the ones of Lemma [TOl and the variables v are assumed as parameters. Since the 
a'^ are linearly independent, this implies K^cr = C'^^a'^, where are smooth functions of the 
variables v only. The expressions for the cross derivatives Ky<T^<. = C^/a* = Clla^ result in the 
equation C^<f = C^l which can easily be integrated: C^^ = P^a for some smooth function P'^ of 
the variables v. Substituting the expressions for C"^ in the equations on K and integrating, we 
obtain K = a^P^+P^, where 0^ = /3^{t, x) is a solution of equation ([28]) . The latter equality and 
the equation a^K^v" — a^K^o = together imply the equation (aja"^ — 0*^0!) (Pja — P^^) = 0. 
Analogously to Lemma [TOl we can state for any i,j £ NU {0} 

«a^^ - aJa^)(P,^. - P^.) = 0. (33) 
Let M°- denote the {p — 2)-nd order minor of W{a), obtained deleting the a'-th and ?'-th 

J ft 

columns corresponding to the functions 0" and and the i-th and j-th rows, where ^ i < 
j ^ p — \. (Enumeration of the rows from is used again for convenience.) We multiply the 
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equation ([33]) by (— 1)*+J+°"'+'^'Mj!^. and convolve with respect to the index pair The 
Laplace theorem on determinant expansion implies that 

^(«L^^a^: a^^a^')iPv^ - Pv^) = 0- (34) 

Here the symbol means that the functions a'^ and a*" are substituted instead of the func- 
tions a'^ and a'' respectively and we have summation over the indices a and q. For any fixed 
{(T, 7^ {o"',?'} we have W{q.\^„^^„i ^^^^^i) = 0. Since W{q.) 7^ in view of the linear inde- 
pendence of the functions a'^, equation (I34p results in P^^, — P^', = 0, i.e. P'^ = H^s for some 
smooth function H of v. □ 

In view of Lemma [TTl the conserved vector {Ku^K^u — Kux) from Lemma [8] has the form 
{f3^u+DxH, P^u — f3^Ux — DtH). Hence, it is equivalent to the conserved vector {(3^u, P^u — (3^Ux) 
which is also a conserved vector of equation ([7]). 

This completes the proof of Theorem [5j 



7 Simplest potential symmetries 

Since there do not exist any pure potential conservation laws of linear (1 -|- l)-dimensional 
second-order parabolic equations, the study of potential symmetries for class ([1]) is exhausted 
by the investigation of Lie symmetries of potential systems constructed with local conserva- 
tion laws only. Such systems are called potential systems of the first level. At first we classify 
so-called simplest potential symmetries |60j . which are the Lie symmetries of the simplest po- 
tential systems associated with a single local conservation law. As a result, inequivalent linear 
(1 -|- l)-dimensional second-order parabolic equations possessing nontrivial simplest potential 
symmetries are completely described. In the part on connections between different linear equa- 
tions our presentation is very close to the theory of Darboux transformations for linear evolution 
equations |46j . The potential systems obtained by employing an arbitrary finite number of local 
conservation laws are investigated in Sections [9l and [TOl A thorough study of the simplest poten- 
tial systems is necessary for understanding the general case since such systems are components 
of more general potential systems. 

Consider an arbitrary equation of the form ([T]). Introducing the potential v by the single 
conservation law of ([1]) with the characteristic a, we obtain the potential system 

Vx = au, vt = aAux — {{aA)x — aB)u. (35) 

(Let us recall that the characteristic a = a{t,x) is a solution of the adjoint equation (jl9p .) 
More precisely, for the conservation law we take the associated conserved vector of the canonical 
form (jlSp . The following statement justifies our choice of conserved vectors for the construction 
of potential systems. 

Lemma 12. In determining the simplest potential symmetries of equations from the class ([T]), 
it is sufficient to consider only conserved vectors of the form (jlSp . which are canonical repre- 
sentatives of the corresponding conservation laws. 

Proof. Only conserved vectors of the minimal order (with fluxes of order and flows of order 1) 
lead to nontrivial potential symmetries of equations from class ([I]). Indeed, consider a conserved 
vector of a higher order of an equation of the form ([T|). It necessarily possesses the representation 

{au + DxH, -aAux + {{aA)x - aB)u - DtH) , (36) 

where a = a{t,x) is a solution of the adjoint equation (jlSp and H is a differential function. In 
view of equation ([T]) the function H can be assumed to depend on t, x and derivatives of u with 
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respect to only x up to an order r, r > 0. Let Q = rdt + ^dx + r]du + Ody be a Lie symmetry 
operator of the corresponding potential system 

Vx = OLU + DxH, vt = aAux — {{aA)x — aB)u — DfH. 

The coefficients of Q are functions of x, u and v. Applying the infinitesimal invariance condi- 
tion [49^ [ST] to the first equation of the system and collecting coefficients of the unconstrained 
variable Vr+i = d^^^v/dx^'^^ , we obtain the equation t]^ — TyUt — S,vUx = 0. After splitting it 
with respect to derivatives of u, we have the condition rj^ = Ty = = 0, i.e., the operator Q is 
projectable to the space of {t,x,u). 

Each conserved vector of the minimal order (equal to 1) of equation ([1]) possesses the repre- 
sentation ()36p . where the function H depends only on t, x and u and, therefore, can be neglected 
due to the point transformation t = t, x = x,u = u and v = v — H , which has, up to similarity, 
no influence on Lie symmetries of the potential system. This finally gives the canonical form (jlSp 
of conserved vectors of equations from class ([TJ. □ 

The initial equation ([1]) for u is a differential consequence of system ()35p . Another differential 
consequence of (fSSl) is the equation 

Vt = Avxx + (^B-Ax- 2^^) Vx (37) 

on the potential dependent variable v, which is called the potential equation associated with the 
linear parabolic equation ([T]) and the characteristic a. There is a one-to-one correspondence 
between solutions of the potential system and the potential equation due to the projection 
(m, v) ^ V on the one hand and due to the formula u = Vx/a on the other. The correspondence 
between solutions of the initial equation and the potential system is one-to-one only up to a 
constant summand. 

Any linear system of the general form 

Vx = au, Vt = Pux + ^u, (38) 

where a, (3 and 7 are functions of t and x, af3 7^ 0, is the potential system of an equation from 
class (H]). Moreover, the corresponding initial and potential equations are uniquely determined as 
its differential consequences. The coefficients of the initial equation are defined by the formulas 

a a a 

and, therefore, a is a characteristic of the conservation law associated with the potential sys- 
tem psp since it satisfies the adjoint equation ()19p . The potential equation is represented in 
terms of (a, 7) as 

Vt = —Vxx + (— - P-^) Vx- (39) 

Lemma 13. Given arbitrary smooth functions a, (3 and 7 of the variables t and x, where 
a(3 7^ 0, the maximal Lie invariance algebras of system (j38p and equation (j39p are isomorphic. 
Namely, for any infinitesimal Lie symmetry operator Q = rdt + S,dx + Ody + rjdu of system (I38p 
its projection Q' = rdt + 5,dx + Ody to the variables (t, x, v) is an infinitesimal Lie symmetry 
operator of equation ([39]) . The coefficient rj is expressed via coefficients of the first prolongation 
of the operator Q' with respect to x. 
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Proof. A heuristic argument in favour of the lemma is the one-to-one correspondence between 
the solutions of the system and the equation. This argument is not sufficient since a point 
transformation on a set of solutions may induce a non-point transformation on an equivalent 
set. That is why the best way to proceed is by direct calculation. 

The infinitesimal invariance criterion [491 [51] implies the determining equations for the coef- 
ficients of a symmetry operator Q of system ([38]): 



Tx — Tu — Tv — 0, — — 0, Oil — 0^^ — 0, 

a \a/t \a/x a \a 

ix-r i — n-F— . (40) 

a a J a 

The equations Tu = = = ^ guarantee that the operator Q is projectable to the variables 
{t,x,v). The other determining equations excluding the last one form the complete system of 
determining equations for symmetry operators of equation ()39p . For any Lie symmetry operator 
Q' = rdt + S,9x + 09^ of equation ([39]) the operator Q = Q' + rjdu, where rj is defined by ([lO]) . is 
a Lie symmetry operator of system (I38p . Equation (I40p is re-written as 

f ^ nx ^"a; 

r]= { -6 - T^Vx - i^Vx 

where 6^ = 6x + O^Vx — (,xVx, i-e. t] is expressed via coefficients of the first prolongation of the 
operator Q' with respect to x in view of the equation u = Vx/ct- □ 

In view of Lemma [13] the investigation of Lie symmetries of any potential system associated 
with an equation from class ([1]) is reduced to the consideration the corresponding potential 
equation. The only possibility for pure potential symmetries to appear is connected with the 
coefficient rj. Namely, the condition r]^ ^ should be satisfied. Therefore, the Lie symmetry 
operator Q' of equation ([39]) induces a potential symmetry of the associated equation from 
class ([I]) iff (9^x / 0. 

It is convenient to use the other dependent variable w = v/a instead of v in the further consid- 
eration. The function w will be called the modified potential associated with the characteristic a. 
In terms of w the potential equation ([37p takes the form 

wt = Awxx + Bwx + Cw, (41) 

which will be called the modified potential equation, associated with characteristic a. Here 

at , , Q-xx . . N "a; o . / "a; 2 



A = A, B = B-Ax, C = - + A— + {B-Ax)—-2A 

a a a \ a 

Introducing the function tp = 1/a, we rewrite the equality for C in the form 

^pt - Mxx - {B - Ax)tpx 



C 



i.e. is a solution of the same equation (|4ip as w. (The value w = 1/a obviously is a solution 
of ()4ip since u = 1 is a solution of ([37p .) The potential system ([35p can also be written in terms 
of w and ^ instead of v and a: 

Wx - '^w = u, Wt - '^w = Aux + ( + B - Ax] u. (42) 
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The above representations of the potential equation and potential system are more suitable for 
the classification of Lie symmetries. In fact, the first equation of (j42p is the Darboux transfor- 
mation [46] of (j41|) to (dl). The Darboux transformation possesses the useful property of duality. 
We formulate this in a way slightly different from [46] . Denote the Darboux transformation 
constructed with the nonzero function ip by DT[?/;], i.e. 

Lemma 14. Let be a fixed nonzero solution of (j41|) and let the Darboux transformation 
T)T[w^] map ()4ip to equation ([T|). Then = is a solution of the equation ()19p adjoint to 

equation ([T|) and DT[a''] maps (jl9p to the equation adjoint to (j4ip . i.e., 

DT[u;"] ^ ^ ^ 

ut = Auxx + Bux + Cu < wt = Awxx + Bwx + Cw 

t 

at + {Aa)xx - {Ba)x + Ca = > at + {Aa)xx - {Ba)x + Ca = 0. 

Note 17. DT[a'^] will be called dual to the Darboux transformation DT[u)'']. Since the twice 
adjoint equation coincides with the initial one, the twice dual Darboux transformation is nothing 
but the initial Darboux transformation. Moreover, 'then' in the lemma can be replaced by 'if 
and only if. 

Lemma O can also be reformulated in terms of characteristics of conservation laws. Denote 
equations and (f^TT) by C and £ for convenience. 

Lemma 14'. If is a nonzero solution of C and DT[tt;°](£) = £ then a° = £ Chf(£) 

and DT[aO] : Chf(£) ^ Chf(£ ). 

Note 18. For any nonzero solution tp of (j4ip the Darboux transformation DT['(/;] is a linear 
mapping of the solution space of (j4ip to the solution space of ([1]) . The kernel of this mapping 
coincides with the linear span (ip). Its image is the whole solution space of ([1]). Indeed, for any 
solution n of ([T]) we can find a solution w of (I4ip . mapped to u, by integrating system ()42p with 
respect to w. System ()42p is compatible in view of equation ([1]). Therefore, DT[t/)] generates a 
one-to-one linear mapping between the solution space of (j41l) . factorized by {^p), and the solution 
space of (dl). 

In view of Proposition [7] any point equivalence transformation in class ([1]) is prolonged to 
characteristics of conservation laws of equations from this class. Then it can be prolonged to the 
first derivatives of potentials due to potential system (|35p , and under this prolongation Vx = Vx, 
vt = vt. Excluding trivial translations of the potential variable v with constant summands, we 
assume that v is not transformed under prolongation of the equivalence transformations. The 
extension to the functions w and ip is obvious. As a result, we obtain simultaneous equivalence 
transformations between the initial, potential and modified potential equations: 

i=T{t), x = X{t,x), u = U^{t,x)u, 

""xlTi' * = ^-^'^' i^ = XxU^i^. (43) 

where TtXxU^ / 0. 

Transformations (j43p preserve the determining equations derived in the proof of Lemma [13] 
as well as both the conditions O^x 7^ and O^x = 0, where 6 is the coefficient of in symme- 
try operators of potential systems or potential equations. This is why Lie and pure potential 
symmetries of any equation from class ([T]) are not mixed under transformations (143 p . and the 
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dimension of the factor-space of potential symmetry operators constructed with a single char- 
acteristic with respect to the subspace of Lie symmetry operators is not changed. Therefore, 
potential symmetries of equations from class ([1]) can be studied up to the equivalence relation 
generated by transformations from C. In particular, it is sufficient to consider only the re- 
duced form ([7]) of linear parabolic equations. The corresponding potential system and potential 
equation and their modifications also are simplified to 

Vx = au, vt = aux — oixU, (44) 
ax 

Vt - Vxx + 2 Vx = 0, 

a 

wt - Wxx + Pw = (45) 
Wx -w = u, Wt -w = Ux + —u. (46) 

where a = a{t, x) is a fixed solution of the reduced adjoint equation at + axx — Va = 0, 
ip = ip{t,x) is a fixed solution of ([451) . V = V{t,x) and 

p = V-2(—) =V + 2(— 



a ) X V -0 

The equivalence transformations dH]) in the reduced class d?]) form a subgroup in the equiv- 
alence group under special restrictions on the parameter- functions X and . They are 
prolonged by formulas (j43p to equivalence transformations of the whole reduced potential frame. 
Hence, the classification of potential symmetries of class ([7]) follows from the group classification 
of the same class in terms of (w;, P) instead of (u, V). 

Consider a Lie symmetry operator Q' = rdt + S,dx + Qdw of equation (|45]) . The coefficients 
of Q' are functions of t, x and w. The infinitesimal invar iance criterion |491 [5T] applied to 
equation (I45p implies the conditions 



1 11 

r = T(t), = -Ttx + a{t), C = Q^{t,x)w + C^{t,x), = --mx^ - -atx + ^{t) 

which do not involve the arbitrary element P, and the classifying equations 
Cl - Qlx + rPt + CPx + TtP = 0, - Cxx + = 0. 

Due to the first equation of system (I46p . the coefficient r] of du of the corresponding Lie symmetry 
operator Q of system ()46p can be expressed via coefficients of the first prolongation of the 
operator Q' with respect to x: 



^0 _ aO 



It is obvious that the operator Q' induces a pure potential symmetry of equation ([7]) iff ry^, / 0, 
i.e., 

r0t+CQx + ^Ttg^CL (47) 

where g = tpx/tp- This implies that trivial Lie invariance operators of the potential equations 
generate no pure potential symmetry operators of the initial equations. More precisely, the 
Darboux transformation DT [-;/'] maps the trivial part of the Lie invariance algebra of (|4ip / (|46p 
onto the trivial part of the Lie invariance algebra of ([I])/ ([7]). We therefore obtain the following 
result: 
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Proposition 11. Equation ([7]) admits a pure potential symmetry associated with the character- 
istic a only if the corresponding equation ()45p has a nontrivial Lie symmetry, i.e., P is equivalent 
to a stationary function. 

It may happen that nontrivial Lie invariance operators of the potential equations do not gen- 
erate any pure potential symmetry operators of the initial equations. We provide the symmetry 
criterion for Lie invariance operators and solutions of the potential equation leading to pure 
potential symmetries. Its formulation is independent of the equivalence transformations. 

Proposition 12. Suppose that Q' = rdt + S,dx + C^wd^ is a nontrivial Lie invariance operator 
of equation (|4ip . tp is a solution of (I41j) and DT['i/'] transforms equation (j4ip to ([T]). The 
operator Q' induces a pure potential symmetry operator of equation ([T|) iff the following obviously 
equivalent conditions are satisfied, where Q'[ip] = C^ip — Tipt — d^x-' 

1) Q' [ip] and ip are linearly independent; 

2) "0 is not an eigenfunction of Q' ; 

3) DT['0](Q'['0]) is a non-vanishing solution of equation ([T]); 

4) ^ is not invariant with respect to the operator Q' — Xwdw for some constant A. 

Proof. All conditions in the proposition are stable under equivalence transformations. It is 
enough to consider the reduced equations ([7]) and (fl6l) instead of (f^TI) and ([T]). Since Q'[ip] is 
the result of the action of the Lie invariance operator Q of equation (j46p on the solution ^ of the 
same equation then Q'lijj] also is a solution of (fl6ll and hence DT[V'](Q'[V']) is a solutions of ([7]). 
Let Q denote the prolongation of Q' to u. The expression for r/^, where r] is the coefficient of du 
in Q, is rewritten in terms of Q'[ip]: 

Therefore, r/^ = iff W^tp, Q'['4']) = 0, i.e. DT = 0. In view of Lemma[71 the condition 
W{^p, Q'[^p]) = is equivalent to the fact that the functions Q'[ip] and tp are linearly dependent, 
i.e. Q'lip] = Xip for some constant A or QaIV'] = 0) where Q'y^ = Q' — Xwdw For any A the 
operator belongs to the Lie invariance algebra of equation ([16|) . □ 

In view of Lemma [13] pure potential symmetries obviously exist if the dimension of the 
essential Lie invariance algebra of the potential equation is greater than the corresponding 
dimension for the initial equation. As shown by the statements below and the example considered 
in Section [51 this condition is not necessary. 

Theorem 6. A linear second-order parabolic equation admits simplest pure potential symmetries 
iff it is equivalent with respect to point equivalence transformations to an equation from class ([7]) 
in which 

where ip = 'ip{t,x) is a nonzero solution of the equation ipt — ip^x + P{x)'ip = and either P = 
/u = const, or {ipx/ip)t 7^ (this latter condition is equivalent toVt^O) if P is inequivalent 
to fix~'^ with respect to point equivalence transformations. In the case P = the initial equation 
possesses at least two potential symmetry operators which are linearly independent up to Lie 
symmetries. The associated characteristic equals l/ip. 

Proof. Let us briefly recall the main results of the above consideration which form the basis 
of our proof. The transformations from G~ generate a well-defined equivalence relation on the 
whole frame of simplest potential symmetries of equations from class ([T]). The investigation of 
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Lie symmetries of potential systems is reduced to the group classification of modified potential 
equations forming the same class as the initial equations and are connected with them via the 
Darboux transformation. To complete the proof, we have to verify condition (j47p for all the 
cases of the Lie~Ovsiannikov classification (see Theorem [2|) . 

Let P = Suppose that any of the operators dt, D and n does not satisfy condition ()47|) 

for some il). Then Qt = 0, xQx + = and AtxQx + ^tg = —2x, which implies the contradiction 
= a;. Hence for any solution of equation (l45ll at least one operator among dt, D and H satisfies 
condition and, therefore, induces a pure potential symmetry operator of ([7]). 

In the case /i = at least one operator among dx and G induces a pure potential symmetry 
operator of ([7|) since otherwise the equations Qx = and 2tQx = — 1 would imply the contradic- 
tion = — 1. Therefore, for P = we have at least two independent pure potential symmetry 
operator of 

There is only one independent nontrivial Lie symmetry operator dt for the general value 
P = P{x). If it does not induce a pure potential symmetry operator of d?]) then gt = and, 
therefore, Vt = Pt + Qxt = 0. Conversely, suppose that Pt = Vt = and gt simultaneously. 
Then Qxt = 0. If equation ([15|) has a solution satisfying the conditions {ipx/il>)t / and 
{ipx/'ip)tx = then it is equivalent to the linear heat equation (the case P = 0). □ 

Based on Theorem [6l we can formulate symmetry criteria on the existence of simplest poten- 
tial symmetries without involving equivalence transformations. 

Corollary 21. A linear second-order parabolic equation admits simplest pure potential symme- 
tries associated with its characteristic a only if the corresponding potential equation possesses 
nontrivial Lie symmetry operators. If the potential equation has more than one (three) inde- 
pendent nontrivial Lie symmetry operators then the initial equation admits at least one (two) 
independent simplest pure potential symmetry operators. 



8 Simplest potential symmetries of the linear heat equation 

Theorem [6] gives description of equations from class ([T]) , having nontrivial simplest potential 
symmetries. At the same time, there is another problem concerning simplest potential symme- 
tries: Given an equation from class ([T]), to describe all its characteristics leading to its nontrivial 
simplest potential symmetries. In what follows we consider this problem in detail for the linear 
heat equation 

ut = Uxx- (48) 

We recall that the maximal Lie invariance algebra of the linear heat equation is 

So = {dt, dx, 2tdt + xdx, 2tdx - xudu, 4t^(?i + Atxdx - (x^ -\- 2t)udu, udu, fdu). 

Here the function / = f{t,x) runs through the solution set of this equation. The operators 
from 00 generate the continuous symmetry transformations of (j48p . See, e.g., |49j for their 
explicit form. Equation (j48p also possesses a nontrivial group of discrete symmetry transfor- 
mations generated by two involutive transformations of alternate sign (t, x, u) (t, —x, u) and 
(t,x,u) {t,x,—u). The point symmetry group Go of ([l8j) is generated by both the continuous 
and discrete symmetries. The most general solution obtainable from a given solution u = u{t,x) 
by transformations from Go is of the form 

u= , =e 1+4-6* u{ £2,^^^ - ei] + f(t,x), 
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where ei, . . . , eg are arbitrary constants, £364 7^ and / = f{t, x) is an arbitrary solution of the 
hnear heat equation. The essential part Gq^^ of Go is formed by the transformations with / = 0. 

In view of Corollary [10] the transformations from Gq^^ prolonged by formulas (f43l) form 
equivalence groups on the sets of potential systems and potential equations associated with 
equation (j48]l and single characteristics. We identify these equivalence groups with Gq^'^. It is 
natural to investigate problems on potential symmetries of (j48p up to the equivalence relation 
generated by G^^. 

Lemma 15. The characteristic a of the linear heat equation gives an equation of the form ()45|) 
with Px = iff a = 1 mod G\ 



ess 

• 



Proof. Let = in the equation (j45p corresponding to the characteristic a. Integrating the 
equation —2{ax/a)x = P as an ordinary differential equation on a with the parameter t, we 
obtain the expression a = ^exp(Px^/4 + Q, where P, 9 and (" are functions of t and 6^0. 
Substituting of the expression for a to the the backward linear heat equation at + Oxx = 
and the subsequently splitting with respect to x result in the system of ordinary differential 
equations 

Pt = -P\ Ct = -PC, Ot = (c' + \p^ 0. 

In integrating this system, two cases are to be distinguished, corresponding to two cases for the 
characteristics: 

V = 0, a = 62 exp((^ix — 6ft) and V = —, a = 62 exp — — , 

t + oo A[t + do) 

where 5q, 61 and 62 are constants. In both these cases the characteristics are Gq*^*^ -equivalent to 
a = 1. 

The converse statement is obvious. □ 

Theorem 7. Lie symmetries of system (|44p . where a = a{t,x) is a (fixed) solution of the 
backward linear heat equation at + axx = 0, induce pure potential symmetries of the linear heat 
equation iff a £ {l,a;} mod G'q^. 

Proof. In view of Proposition [TTl a characteristic a leads to pure potential symmetries of ()48p 
only if the function P = —2{ax/a)x is G^" -equivalent to a function not depending on t. Let 
us apply a transformation T from G^, prolonged to the whole frame under consideration, such 
that Pt = 0. Equation is mapped by T to the equation ut — Uxx + Vu = 0, where Vxxx = 0. 
(We omit tildes over the transformed variables t and x for convenience.) The functions V and P 
are connected by the relation P — V = —2{ax/a)x = 2(^^/^)^, where a is the transformed 
characteristic and ip = 1/a. Considering t as a parameter, we integrate this relation as an 
ordinary differential equation on ■0: t/j = -ize'^^ ^ , where H is a. smooth function of x such 

that H' = P and P is a fourth-degree polynomial of x with coefficients depending on t such that 
Rxx = y ■ Substituting the obtained expression for ij: into the equation ■0* — i^xx + Pi^ = gives 
the following relation between H and R: H' — + RxH = Rt — Rxx + \Rx^ ■ The differential 
consequence 

RtxH = Rtt — Rtxx + RxRtx 

of this relation with respect to dt is essential for the further proof. There are two possible cases 
which should be considered separately. 

If Rtx 7^ 0, we express H from the differential consequence: 

Rtt RtTT ^ 

H=-^--^ + Rx, 

J^tx ^tx 
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and then differentiate once more with respect to t. Since = 0, we derive an equation only 
in R, which can be spht with respect to x since R is a polynomial of x. The resulting equations 
on the coefficients of R imply that Rtxx = and {Rttx/ Rtx)t = 0. Therefore, H is a third-order 
polynomial of x and 

P-V = H' -Rxx = ^ = const . 

Rtx 

This means that P is a function only of t in the old variables. Hence, a = 1 mod Gq^^ in view 
of Lemma [TSl 

The condition Rtx = implies Ru = and = {x)e~'^'' , where v is a constant and the 
function tp^ depends only on x. Since {ipx/ip)t = 0, in view of Theorem [6] the function P should 
be G J" -equivalent to the function fix~^ with a constant Otherwise there are no nontrivial 
potential symmetries associated with the characteristic a. Therefore, 

V = K2X^ + KiX + Kq, P = ^ + ll2X^ + HlX + /io, 

a = — = Xo\x\ f^'^ exp — — — x^ H — — x'^ H + Xix + i^t ] , 

ip \ 24 12 4 J 

where ^, //j, Kj, Ai, Aq = const, i = 1,2, 3. Then the equation at + dxx — Vo: = implies that 

K2=/i2, ni = ^ii, /i(/x-F2) = 0, ;uA = 0, {I - - iiq) = 2{kq - u - \^), 

Ki = A(ko - /io), 4k2 = (ko - /"o)^- 

const considered above. If 7^ then it follows 
= /ii = 0, 1/ = (ko /Uo)/2, i.e.. 



X. (The explicit form of the variable transfor- 

□ 

Potential symmetries associated with the characteristic a = 1 were studied by a number of 
authors [T2l[70l[6lJ. The corresponding potential is defined by the system = u, v} = Ux- 
Its maximal Lie invariance algebra is 

pi = {dt, dx, 2tdx - {xu + v^)du - xv^d^i, 2tdt + xdx - udu, 

Afdt + Atxdx - {{x^ + et)u + 2xv^)du - (x^ + 2i)v^d^^, udu + v^d^i, fxdu + fd^i). 

The potential equation Vt = v^.^ has the same form as the initial equation ()48p . That is why the 
algebras go and pi are isomorphic [61]. The basis operators of pi are obtained from the basis 
operators of go by re-denoting u ^ and then carrying out the first prolongation with respect 
to X in view of n = Any linear combination of operators from pi which contains the third 
or fifth basis operator is a pure potential symmetry operator of the linear heat equation. 

The case of the simplest nonconstant characteristic a = x was studied in |38j. The corre- 
sponding potential system v"^ = u, v1 = xux — u possesses the maximal Lie invariance algebra 

p2 = {dt, 2tdt + xdx - 2udu, ^t^dt + 4txdx - {{x^ + 6t)u + 2v^)du - (x^ - 2t)v'^dy2, 
udu + v'^d^2, x'^hxdu + hdui), 

where the function h = h{t, x) runs through the set of solutions of the potential equation 
Vt + 2x~^v'^ — v'^^ = 0. Any linear combination of operators from p2 which contains the third 
basis operator is a pure potential symmetry operator of the linear heat equation. 



The condition fi = leads to the case P — V = 
from the above equations that /i = —2, X = ki 

a = Aoxexp ( -(kq - /xo) X H tj 

After returning to the old variables, we have a 
mation depends on the value oi kq — fj.Q.) 
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Note that the pure potential symmetry operators from the algebra p2 differ from the ones from 
the algebra pi by both the explicit form and the nature of the potential variable associated 
with the characteristic a = x in contrast to associated with the characteristic a = 1. 

Corollary 22. All Cq^ -equivalent simplest pure potential symmetries of the linear heat equa- 
tion (|^8]) are exhausted by the operators from pi and p2 satisfying the condition r]y ^ 0, where rj 
is the coefficient of du and v is the corresponding potential (v^ or v^). 

Let us demonstrate how to extend the results obtained in this section to equations which are 
equivalent to the heat equation with respect to point transformations. Consider the Fokker- 
Planck equation 

ui = Hxx + {xu)x- (49) 
Its maximal Lie invariance algebra is 

00 = (<9t, e~*dx, e-'^^dt - e-'^^xd^ + e~'^^udu, e^d^ - e^xudu, 
e^^dt + e'^^xdx - e^^x^udu, udu, fdu), 

where the function / = f{t,x) runs through the solution set of the same equation. (We omit 
tildes over variables when it is understandable that they originate from the Fokker-Planck 
equation.) Equation (|^9|) is reduced to equation (|18|) by the point transformation 

T : t = -e^*, X = e*x, u = e~^u. 

According to Proposition [71 the transformation T acts identically on values of characteristics: 
a = a. That is why in view of Theorem [7] only two Gg'^'^-inequivalent characteristics = 1 and 
dP' = e*x of equation ()49p lead to nontrivial potential symmetries of this equation. Here Gq"^ is 
the essential part of the point symmetry group of equation (I49p . The corresponding potential 
symmetry algebras pi and p2 can be obtained in two ways. The first way is the direct calculation 
of the maximal Lie invariance algebras of the associated potential systems v\ = u, = Ux + xu 
and = e*xn, 5? = e^xux + e*(x^ — l)u. The second way is to map the potential symmetry 
algebras pi and p2 of the linear heat equation by the transformation inverse to T and trivially 
prolonged to the potentials and w^, respectively. Finally, the algebras pi and p2 have the 
form 

pi = ( dt, e~*dx, e~'^*dt - e~'^*xdx + e~'^^udu, e^dx - e*{xu + v^)du - e^xv^d^i, 
e^^dt + e^^xdx - e^\x^u + 2xv^ + 2u)du - e^\x^ + l)u^a^i, udu + v^d.,x, 
gxdu + gd^i ), 

h = {dt- udu, e-^'dt - e-^'xdx + e-^'udu, 

e^^dt + e^^xdx - e^\x^u + 26'^ + 2u)du - e^\x^ - l)v^9„2, udu + v'^du2, 
e^x~^hxdu + hdy2 ), 

where the functions g = g{t, x) and h = h{t, x) run through the solution set of the associated 
potential equations — — xv\ = and — + (2x~^ — x)i;| = 0, respectively. 
As a result, we obtain the following statement. 

Corollary 23. All -equivalent simplest pure potential symmetries of the Fokker-Planck 
equation (|49p are exhausted by the operators from pi and p2 satisfying the condition rji, ^ 0, 
where fj is the coefficient of du and v is the corresponding potential (v^ or v^). 

Corollary [23] essentially generalizes results of [381 IMl [67] on simplest potential symmetries of 
equation ([49]) . associated with the characteristic 1. 
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9 Preliminary analysis of general potential systems 



The investigation of general potential symmetries can be carried out in a way similar to that of 
simplest potential symmetries but all calculations are much more complicated. The idea is to 
find an analogue of the potential equations ()4ip in the case of general potential systems and to 
construct potential symmetries of the initial equations via prolongation of Lie symmetries of the 
potential equations. It is necessary also to prove statements on the behavior of symmetries of the 
potential frame under the action of equivalence transformations. This allows us to replace the 
study of class by the study of the simpler class ([7]). Since the problem is quite complicated, 
in this section we only find a convenient form of potential systems and construct potential 
equations each of which is equivalent, in a certain sense, to a whole potential systems. Along 
the way, we make intensive use of multiple Darboux transformations. These results create a 
basis for the symmetry analysis of the potential frame in the next section. 

Let us fix an equation ut = Auxx + Bux + Cu from class ([T|), an arbitrary p G N and 
p linearly independent solutions a^, . . . , of the adjoint equation (I19p . Then the conservation 
laws corresponding to these characteristics are linearly independent. For any s we introduce 
the potential using the conserved vector of the canonical form (jlSp . associated with a^. Our 
choice of the conserved vectors will be justified below by Corollary [30j As a result, we obtain 
the potential system 

= a'^u, vf = a'^Aux - {{a^A)x - a^B)u (50) 

corresponding to the characteristic tuple a = (a^, . . . , a^). Let us recall that the indices s, a 
and q run at most from 1 to p. Additional constraints on indices are indicated explicitly when 
needed. The summation convention over repeated indices is used unless otherwise stated or it 
is obvious from the context that indices are fixed. 

If the initial equation runs through class ([T]) and , . . . , runs through all linearly indepen- 
dent solutions of the adjoint equation then the associated potential systems form the potential 
frame of order p (and the first level) over the class ([T]). In general, by the (potential) order of an 
object we will mean the number of independent first-level potentials appearing in this object. 
Below we extend the potential frame with other objects. 

System (jSOp is homogeneous with respect to the index s. It is formed by p similar blocks. 
Each of them consists of a pair of equations in a potential and the initial unknown function u 
and possesses the structure of a simplest potential system. All the potentials v^, . . . , are 
on an equal footing. At first glance these features seem to be advantages of this representation 
of potential systems but after careful consideration a number of drawbacks become apparent. 
The total number of unknown functions in system (j50p equals p + 1 and the total number 
of equations is 2p. At the same time, the system is not "too" overdetermined since it has 
no nontrivial differential consequences. For any fixed s the corresponding pair of equations 
implies an equation only with respect to and a differential consequence equivalent to the 
initial equation. In fact, system (j50p contains only p + 1 independent equations but the num- 
ber of equations cannot be reduced to the minimal one in a symmetric way. It is not clear 
what a potential equation corresponding to the whole system (|50p should be. Another argu- 
ment in favour of modifying system (jSOp comes from group analysis. Consider a Lie symmetry 
operator Q = rdt + S,9x + V^u + O^d^s of system (j50|) . The coefficients of Q are functions 
of t, X, u and v'^. The infinitesimal invariance criterion applied to system ()50p implies, in 
particular, the following determining equations on the coefficients of Q: = Tx = Tt,s = 0, 
iu = — ^1 — 0- contrast to the case of single characteristics, the deduction of these 
simplest determining equations is much more involved. Certain tricks involving linear inde- 
pendence of a^, have to be used. Finding other simple determining equations which 
are typical for linear systems (e.g., r]uu = 0, r]uv^ = 0, 9la^<i = 0) demands still more calcu- 
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lations and tricks. A complete analysis of the whole system of determining equations seems 
impossible. 

Below by an iteration procedure we obtain another potential system associated with the 
characteristic tuple a = (a^, . . . ,a^), which is equivalent to system (fSOj) but is appropriate for 
the investigation of potential symmetries of the initial equation. Each step of the iteration 
procedure is similar to the consideration in the beginning of Section [71 For convenience an 
equation for an unknown function ip will be denoted by \ip\ . 

Step 1. Let us re-denote u ^ w'^ , B ^ B° , C ^ (7°, f\ a ^ (3^ and a' /JO-^ (We 

will use the old and new notations simultaneously.) Consider the (first-level) potential system 

fl = fl = P^'^Au^l - ((/3°'iA). - /3°'ii?°)u;° 

associated with the single characteristic = /3^'^. The tuple {w^, f^) is a solution of this system 
iff the modified potential = /P^'^ satisfies the equation = Awl.^ + B^w^. + C^w^, where 
B^ = B^ - and 

C' = C-B^ + A^^ + Afp + 2A[^^^y=C-B^ + A^^ + A^^^+2A (^^j^. 

Here and occasionally below the notation for the Wronskian W{a^, . . . ,a^) is used. The 
function it;^'^ = l//?"'^ is a solution of \w^\ . The Darboux transformation DT[i(;^'^] maps \w^\ to 
\w^\ = \u\. In view of Lemmall4lthe dual Darboux transformation DT[/3'^'^] maps \a\ = \f3^'^\ 
to the equation (3} = {A(3^)xx + {B'^l3^)x + adjoint to \w^\. Therefore, the functions 



(5^^' = DT[/3°'i](/3°'^) = al-^a 



satisfy the equation \I3^\ and G Chf([u;^J), i.e. they are characteristics of conservation laws 
of \w^\ . Note that (3^'^ = 0. 

Step 2. Using the conservation law of [w"*^], having the characteristic we introduce the 
potential and obtain the potential system 

fl = = (5^^^ Awl - ((/?^''^). - P'^'B')w\ 

(Its union with the constructed first-level potential system results in a second-level potential 
system.) The tuple {w^,f'^) satisfies this system iff the modified potential w'^ = /^//3^'^ is a 
solution of the equation Wt = Aw^^ + B'^wI + C^w^, where B"^ = B^ - A^ = B - 2Ax and 

C^ = C^- Bl + Axx + Ax^ + 2a($^] =C- 2Bx + SA^. + aM^ + 2A 



since = W"^. The function u)^'^ = is a solution of \w'^\. DT[ui^'^] maps \w'^\ 

in \w^\. Then the dual Darboux transformation DT[/3^'^] maps \(3^\ to the equation /3| = 
{Al3^)xx + {B^l3^)x + C2/32 adjoint to \w^\. Therefore, the functions P^'' = BT[P^'^]{P^'') 
satisfy the equation [/S^J and f3'^'^ G Chf(['u;2j), i.e. they are characteristics of conservation laws 
of \'w'^\. Since P^''^ = DT[a"^](a'^) then in view of the Crum theorem [201 

This formula implies, in particular, that /J^'^ = /32>2 _ q_ 
The next iteration is obvious. 
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Step s. Using the conservation law with the characteristic /?^~^'* of the modified potential 
equation [ti;^~^J from the previous step, we introduce the potential /^^ and obtain the potential 
system 

= p'-^''w'-\ // = p'-^''Aw'^-^ - ((r-i'M)^ - /3'-^''B'-^)w'-\ (51) 

(Its union with the (s — l)-level potential system constructed during the previous iterations 
results in an s-level potential system of[TJ) The tuple {w'^~^,f^) satisfies system (f5T]) iff the 
modified potential w'^ = f^/f3^~^''^ is a solution of the equation = Aw^^ + B'^w^ + C^w^, 
where = B"'^ - = B - sA^ and 

nS—l,S / nS — l,S 



Sjs-l) A OA f^x 

C - sB^^ - + A^— + 2A y^^^ 

since /3°'^ . . . = W^. The function w^''^ = 1//?*"^''* is a solution of \w'^\. DT[u;*'*] 
maps [iu^J in [u;''~"'^J . Then the dual Darboux transformation DT[/3'^~"'^'*] maps [/3^~^J to 
the equation = {Ap^)xx + {B^l3'^)x + adjoint to \w^\. Therefore, the functions 

^ -p:pj^s-i,s](^s-i,a) gj^^jgfy gquatiou \I3'\ and (5'^" G Chf([u;^J), i.e. they are charac- 
teristics of conservation laws of \w^\ . Since are constructed by iteration of the Darboux 
transformation from the characteristics a^, . . . , a^, in view of the Crum theorem we obtain 

^'-^ = DT[r-i'^](r-i''^) = - = ^/"/•••'"''""^ 

This formula implies, in particular, that (3^'" = if a ^ s and (3^'" 7^ if cr > s. 

The iteration procedure is stopped on step p after the construction of the equation \w'p\ since 
there are no nonzero functions /J^'*^. 

For any s < p the second equation in system (|5ip is a differential consequence of system (]5ip 
with s + 1 replacing s. Therefore, the 'minimal' combined potential system consists of the first 
equations of the potential systems from all steps and the second equation of the potential system 
constructed on the last, p-th, step. 

Let = = u, = = 1 and /?~^'** = 1 by definition. Excluding w** due to the 

formula tf* = /^//?*~^'^, we obtain the combined potential system in terms of only 

n = H'r-\ ff = HPAfs-'-GPr-\ (52) 



where 



0S-I,s ^rs^rs-2 

1^ " " / TJS A \ ITS TJS — l I o ZJS A ^x 



H' = lBzT-T = 77777^^ = {H' A), - B'-' + 2H^ A 



It can be proved that Hf + G| = for any s. 

System (j52p is the level form of the potential system of equation ([T|), associated with the 
characteristic tuple a = (a^, . . . , a^). The equations = H^Af^~^ — G^f^~^, s = 1, . . . ,p — I, 
are differential consequences of (f52]) . In the case p > 1 the derivatives with respect to x can be 
excluded from these equations with s > 1 as well from the last equation of (j52p . The resulting 
equations are = H'H'-^Af'-^ - s = 2,...,p. 

In view of Corollary [^Uj system (|52p should be equivalent with respect to point transforma- 
tions to a potential system of the first level. Below we explicitly construct a point transformation 
from system ([50]l to system (f52]l . 
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Let us define the functions g^'"^, cr ^ s, by the recursive formula 
For convenience it can be assumed that g^'"^ = if cr < s. 

Lemma 16. For any fixed s and a the function g^'^ is the potential of the equation \w^~^\, 
associated with the characteristic /3*~^''^. In particular, g^'^ = up to a trivial constant sum- 
mand. 

Proof. The lemma is proved by induction with respect to s. The statement of the lemma 
for s = 1 is obvious in view of the definitions of g'^''^ and Suppose that the statement is true 
for a fixed s and a = s, . . . ,p. Let us prove it for s + 1 and a = s + 1, . . . ,p. By assumption, 
the functions 5*''^, a = s, . . . ,p, satisfy the conditions 

and g^'^ = Then the first derivatives of g^'^^''^, a = s + 1, . . . ,p, are 

= (^p^'-{B^-'-2A,) - p^A - /?r''^^) + /3^''^^(< + friiJ^') 

Therefore, the function ^*+^''^ is a potential of the equation \w'^\ , associated with the charac- 
teristic P^^". Since c/x"^^'^+^ = f'+^ and 5'+^'*+^ = //+^ then c,^+i.^+i = up to a trivial 
constant summand which can be neglected. □ 

For our further considerations we need a simple but useful property of matrix minors. Let 
'^il 'i'' denote the submatrix of the square matrix M G M„^„, obtained by deletion of the rows 
with (different) numbers ii, . . . , iq and the columns with (different) numbers ji, . . . , jg, q ^ n. 
^il iq ~ ^^^-^ii corresponding minor. In particular, M^^"f^ = 1 by definition as 

the determinant of the empty matrix, detA^ = M {q = 0, i.e., there are no deleted rows and 
columns). Here and in the next lemma all indices run from 1 to n. 

Lemma 17. For any n ^ 2, M E M„^„, i^ j, k ^ I: 

MiMi - MlMl = sign(j - i) sign(i - I) M^^M. 

Proof. The lemma is also proved by induction over n. For n = 2 the statement reduces to 
the formula for detAl. Suppose that the statement is true for a fixed n. Let us prove it 
for n + 1. Consider M. G M^j+i^^+i. After permuting rows and columns, without loss of 
generality we can assume that i = k = n, j = l = n + l. Thus, it should be proved that 
K := M^^M^'+l - M^^^MJi+^ - M'^'^Hm = 0. We expand M^, M^^^ and M with respect to 
the elements of the (n + l)-st column of M: 

n—l 
1=1 
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(The terms with i = n and i = n + 1 cancel.) Then for each i £ {1, . . . , n — 1} we expand iW^^j^, 
M^"*"^ and M^'^^-,^ in the corresponding term of the sum with respect to the elements of the i-th 
row: 



n—l n—1 



i=l j=l 

(The terms with j = n cancel.) The coefficients in the latter sum vanish in view of the induction 
hypothesis applied to the matrices □ 



Lemma 18. g^'" = (-1)' 



_^ W{a^, . . . ,0" ^a'^)^i^a,_l 
PF(ai,...,a^-i) 



Here the notations introduced in Section [6] are used. Namely, the subscript s — 1 denotes 
the {s — l)-st order derivative with respect to x, and a function is considered as its zero-order 
derivative. The notation ~^ Ss-i" means that the derivatives are replaced by the 

function v'' for the range of subscripts in the corresponding Wronskian, i.e., = l,...,s — 1,0" 
in our case. Note that Lemma [18] gives significant values of g^'^ only for a ^ s and g'^''^ = if 
o" < s, in agreement with the definition of 5*''^. 

Proof. We again use induction on s. The statement of the lemma for s = 1 is obvious since 
gi,cT _ ycr ^ W{a'^)v^a = and the Wronskian of the empty tuple equals 1 by definition. 
Suppose that the statement is true for a fixed s and a = s, . . . ,p. Let us prove it for s + 1 and 
o" = s + 1, . . . ,p. In view of the assumption and the recursive formula for g^~^^'^ , the statement 
for s + 1 is equivalent to the formula 

W{a^, . . . , a'*) W{a^, a''"\ a'^)i}^a,_i - W{a^, a''"\ a"") W{a^, a-'')j)^a,^i 

The latter is a pure matrix equality and does not depend on the specific structure of Wronski 
matrices. It follows from Lemma [T71 □ 

Corollary 24. Systems (jSOp and (|52p are equivalent with respect to a point transformation of 
only the potential dependent variables, being linear in these variables. In other words, the p- 
level potential frame over class ([T|) is equivalent to the first-level potential frame of order p over 
the same class. They can be simultaneously considered in the framework of the general p- order 
potential frame. 

Corollary 25. On any level s the functions g^"^ , cr ^ s, can be expressed via the functions g''" , 
a' ^ <;^, of any lower level q: 



= (-l)'^-^ 



T¥(/?«-i'S...,/3^-i'^-i) 

Proof. We apply Lemma [TS] assuming that the iteration procedure is started from level q. □ 

Since tu* = f^/P'^~^'^ then in view of the formula for (3^'^'^^ we have one more corollary from 
Lemma 1181 



Corollary 26. For any s the final result of the s-th iteration (i.e., the expression of the potential 
via the potentials , . . . , and the form of the equation \w^\) is invariant with respect 
to nondegenerate linear transformations cf = ^^=1 0('' c^a , cr = l,...,s, of the characteristics 
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from the tuple {a^,...,a^). Here c^^j, a,<; = l,...,s, are constants such that det(c<jo-) 7^ 0. 
In particular, 



a 



Note 19. A nondegenerate linear transformation al^ = X]f=i o^''c<jo-, o" = l,...,s, in a 
characteristic tuple (a^,...,a'') necessarily implies the simultaneous linear transformation 
v'^ = X]f=i '^'"Cfcr in the associated potential tuple {v^ , . . . , v^) with the same coefficients. 

Note 20. We can say that the modified s-level potential w'^ and the equation \w'^\ corre- 
spond to the s-dimensional characteristics subspace (a^,...,a'') instead of the characteristic 
tuple (a^, . . . , a^) since the choice of a subspace basis is inessential in view of Corollary | 



Note 21. The Crum theorem can be proved in a way similar to Lemma [T8] by using Lemma [T71 

Below the notation "(a^, . . . . . . , a^y means that a!^ is absent in the corresponding tuple 

of a's. 

Lemma 19. For any fixed s the functions 

W{a^, . . . , a*) 

are linearly independent solutions of \w^\ . Moreover, W{w^'^ , . . . , w*'*) W{a^, . . . , a**) = 1. 

Proof. Since for any a v'^ = const is a solution of the corresponding potential equation, in view 
of Corollary [26l the tuples {v^, . . . , v^) = {dia, ■ ■ ■ , ^sa), cr = 1, . . . , s, where 5^0- is the Kronecker 
delta, give s solutions w'^''' of \w^\. By the formula on determinant expansion, w'^''^a^_i = if 
q = 1, . . . , s — 1 and w'^''^a'^_i = (—1)''^^ if ? = s. Combining differential consequences of these 
formulas, we derive that 



s,<T a 



0, ^ = 



Therefore, the product of the Wronski matrix of (a^ , . . . ,a^) and the transposed Wronski matrix 
of {w^'^, . . . is a matrix with zeros above the anti-diagonal (the diagonal going from the 

lower left corner to the upper right corner). Its anti-diagonal entries equal (—1)''"'^, ^ = 1, . . . , s, 
starting from the lower left corner. The determinant of such a matrix equals 1. Therefore, 
W{w^'', w'^') W{a^, ...,a') = l. □ 

Corollary 27. The s -level potential vu'^ is a linear combination of the potentials , . . . , with 
functional coefficients which are fixed solutions of the equation \w^\ : = w'^^v'^ . 

Lemma 20. DT[u>^'*](tt;'''?) = w^^-^'S ? = 1, . . . , s - 1. 

Proof. We substitute the expressions for w*'* and tf*'^ given by Lemma [19] into T}T[w^'^]{w^''^), 
simplify the obtain result and then apply Lemma [T71 □ 

The multiple Darboux transformation constructed with a tuple of linearly independent func- 
tions (-0"'^, . . . , ■0'^) is denoted by DT[?/;-'^, . . . ,^/''^], i.e.. 

As a result of the above considerations we obtain the following theorem: 
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Theorem 8. Let , he linearly independent solutions of the equation \a\ = \u\* 

adjoint to the equation \u\ and DT[q!^, . . . , a^] [aj = [/3^J • Then the adjoint equation 
\j3P\* = \wP\ to [/3PJ is the p- lev el potential equation of \u\, constructed from the characteristic 
tuple {a^, . . . , aP), the functions 

^_^ W{a^,...,h^,...,aP) 

=(-1 TTTT^ 7a ' ? = l'---'f'' 

W (a'-, . . . , aPj 

are its linearly independent solutions and DT[wP'^, . . . ^w^'^] \wP\ = \u\, i.e., 

ut = Au^^ + Bu^ + Cu i — — wl = Aw'P^ + BPwP + CPwP 

t 

at + (Aa),, - (Ba), + Ca = °^^°''-'"'') /3f + (ApP),, - (BPpP), + CPpP = 0. 

In view of reflexiveness in the duahty of hnear equations, Theorem [8] may also be reformulated, 
analogously to Lemma [T^ in terms of characteristics of conservation laws. 

Corollary 28. // , . . . , ijjP are linearly independent solutions of the equation C from class ([T]) 
and T)T[il}^, . . . , V'^](£ ) = C then C belongs to the class ([1]), 

and DT[a"^, . . . , a^] : Chf(£) ^ Chf(£ ). The tuples {ip^ , . . . , tpP) and (a^,...,aP) will be called 
dual to each other. 

Note 22. A statement on a connection of the same kind between the equations from two 
arbitrary steps of the iteration procedure can be formulated similarly to Theorem [51 It is 
sufficient to assume that the step of lower number is the start of the iteration and the step of 
greater number is the end of the iteration. In particular, 

DTK''^, . . . , w'^'] KJ = Iw"-^] ^ BTif]"-^^" , P"-^''] IP"-^] = \p'\ . 

Note 23. Analogously to simple Darboux transformations (see Note [THj) . for any linearly inde- 
pendent solutions ujP'^, . . . , vjP'P of \wP\ the multiple Darboux transformation DT[t(;P'^, . . . , wP'P] 
is a linear mapping from the solution space of \wP\ into the solution space of \u\ . The kernel 
of this mapping coincides with the linear span {wP'^, . . . ,wP'P). Its image is the whole solution 
space of \u\ since it is the composition of the simple Darboux transformations DT[tt;*'*]. For 
any s DT[w^''^] is a linear mapping from the solution space of \w'^\ onto the solution space 
of [tt;*~^J, where := u. Therefore, DT[wP'^, . . . ,wP'P] generates a one-to-one linear mapping 
between the solution space of \wP\ , factorized by the subspace {wP'^, . . . , wP'P), and the solution 
space of \u\ . 

The equations [nj and \wP\ are differential consequences of the p-level potential system ([52]) . 
There is a one-to-one correspondence between solutions of the potential system and the equa- 
tion \wP\ due to the projection (f^, . . . , fP) fP and the transformation ujP = fP/pP~^'P in 
one direction and due to the inverse transformation fP = pP~^'PwP and the backward recur- 
sive formula /*~^ = fx/H^ in the other. The correspondence between solutions of the initial 
equation \u\ and the potential system is one-to-one only up to arbitrary linear combinations of 
p fixed solutions of system ()52p . This follows, e.g., from the fact that every v'^ is determined 
via u up to an arbitrary constant summand and (/^, . . . , fP) is the product of {v^, . . . , yP) and 
a matrix-function with coefficients depending on t and x. Taking into account the above argu- 
ments and Note [20l we will call the equation \wP\ the modified potential equation, associated 
with the equation \u\ and the characteristics subspace {a^, . . . ,aP). 
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A system of the general form ()52p with . . . H^A 7^ is a p-level potential system of an 
equation from class ([T]) only under special restrictions on the coefficients. Namely, the following 
conditions are necessary and sufficient: 

H! + Gl = 0, HI = {AH%, - (^ ^ ^^^^ - H^j , s<p. 

Here the coefficients G^, s < p, are calculated from G^ and , o" ^ s, by the recursive formula 
G' = G"^' -yj^'^')- Hs _ ^jjs 

Therefore, for any fixed s < p the function should satisfy the Fokker-Planck equations with 
diffusion coefficient A and drift coefficient expressed via A, G^ and , a > s. 

10 General potential symmetries 

Consider a system of the general form (j52p with . . . H^A ^ 0, which is a p-level potential 
system of an equation from class ([T]). Such a system possesses s — 1 algebraically independent 
nontrivial differential consequences = H^Af^~^ — G^ f^~^, s = 1, . . . ,p — 1 having, as differ- 
ential equations, order one. System ()52p implies the second-order partial differential equation 

f! = AfS. - ^^^^fS (53) 



with respect to only the function /P. 



Definition 10. The Lie invariance algebra of a p-order (or p-level) potential system of an 
equation C from class ([1]) is called a p-order potential symmetry algebra of £. Any operator 
from this algebra is called a p-order potential symmetry operator of C 

Lemma 21. Let system ()52p be a p-level potential system of an equation from class ([T]). Then 
the maximal Lie invariance algebras of system (j52]) and equation ([53]) are isomorphic. Namely, 
for any Lie invariance operator Q = rdt + ^dx + r]du + O'^dfs of system (I52p its projection 
Q' = rdt + + O^dfp to the variables {t,x,fP) is a Lie invariance operator of equation (j53p . 
The coefficient 9^~^ of Q, where 6^ := r], is expressed via coefficients of the {p — s + l)-th 
prolongation of the operator Q' with respect to x in a backward recursive way in accordance with 
the equations f"'^ = f^ /II" . 

Proof. Similarly to Lemma [T3l the one-to-one correspondence between the solutions of the sys- 
tem and the equation gives an empiric argument in favour of the lemma. Since this argument is 
not sufficient, the best way again is to make direct calculations. An application of the infinitesi- 
mal invariance criterion [49] IST] to system ()52p has specific features due to the arbitrariness of p 
and the recursive form of the equations f^ = H" f'^~^ and needs a usage of tricks. That is why 
we provide some explanations on the derivation of the determining equations of the coefficients 
of the Lie invariance operators. 

The first trick is to consider also the determining equations obtained via the application of the 
infinitesimal invariance criterion to the nontrivial first-order differential consequences of ()52p . 
This allows us to simplify calculations at an early stage. The trick is admissible since the system 
extended by differential consequences possesses the same Lie invariance algebra as the initial 
system. (This is why usually differential consequences are used only under confining to the 
system manifold.) 

Let Q = rdt + ^dx + rjdu + O^djs be a Lie invariance operator of system ([52]) . Here the 
coefficients r, ^, r] and 9^ are smooth functions of the variables t, x, u and f^. Splitting the 
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infinitesimal invariance conditions of the extended system with respect to the derivatives Uf and 
Ux (there are no other possibilities for the initial split), we obtain, in particular, the equations 
Tu = = du — 0) "Huu = 0, Tx + TftrH'^ f'^~^ = 0. They allow us to further split with respect to u. 
Thus, splitting the already derived condition Tx+TftrH^ f'^~^ = 0, we have Tji = and, therefore, 
allows to also split with respect to in this condition. Iterating the splitting of this condition 
with respect to the dependent variables, we derive Tx = tjct = 0. Other infinitesimal invariance 
conditions imply, under splitting with respect to u, the equations ^fi = and 6*^1^1 = 0. 

Since Tu = = du = and u appears only in the equations corresponding to the value s = 1, 
these two equations with s = 1 have no influence on the further splitting in the infinitesimal 
invariance conditions for the other equations. We can consider the subsystem formed by the 
equations with s ^ 2 (this is the second trick). Here and play the same role as = u 
and in the whole system. Therefore, analogously to the whole system, for this subsystem we 
have the conditions 9p = 0, a ^ 2, =0 and ^^2^-2 =0. As a result of iterating the procedure 
of system contraction, we obtain 9ja = 0, a < s, ^/s = and Ofsp = 0. 

The third trick is that we can neglect the infinitesimal invariance conditions of the differential 
consequences from this moment on. Under the derived restrictions the infinitesimal invariance 
conditions for the equations f'^~^ = fx/H^ imply the backward recursive formula 

us ITS \ iTcr as 

^ I7>S 

with the start in s = p. Therefore, 0faj<; = for any s, a and in view of = 0, a < p and 
= 0. Finally, the infinitesimal invariance conditions for the last equation of system (|52|) 
result in the equations 

{2ix-rt)A = TAt + iAx, n = Aei, ^W^^""' 

6 + C^ePp - Cxx)A = T ^ — j^+C ^ — + Jfp — ' 

The equations Tu = = du = 0, tjs = S,f^ =0 and Oj,^ = 0, a < s, guarantee that the 
operator Q is projectable to the variables {t,x,f^, . . . , /^) for any s. The obtained system of 
determining equations contains, as a subsystem, the complete system of determining equations 
for Lie symmetry operators of equation ([53]) . Therefore, the projection Q' = rdt + S,dx + O^dp 
of Q to the variables {t, x, f^) is a Lie invariance operator of equation (I53p . And vice versa, for 
any Lie invariance operator Q' = rdt + C^x + O^dp of equation ([53]) the operator 

Q = Q' + 7]du + ^e''dp, 

cr<'p 

where the additional coefficients rj = 6^ and 6'^ , a < p, are determined by the above backward 
recursive formula, is a Lie symmetry operator of system (I52p . The formula can be re-written as 



I QS,x _ ^f±L fs _ tf^x fs 
ffs \ ffs S J X 



Hsp 



where there is no summation over s (since it is assumed fixed) and 6^'^ is the coefficient of dfs 
in the first prolongation of the operator rdt + 6.dx + O^dfs. Therefore the coefficient 9^~^ is 
expressed via the coefficients of the standard {p — s + l)-st prolongation of the operator Q' with 
respect to x in accordance with the backward recursive equations f^'^ = fx/H". □ 

Corollary 29. // ^ = for a fixed value of s and any a > s — 1 then 9p^ = for any ? ^ s 
and any a > s — 1. 
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Proof. If = 0, CT > s — 1, then the coefficient 9'' ^ for any q < s is expressed via the 
coefficients of the standard (s — ? + l)-st prolongation of the operator rdt + (,dx + 9^~^dfs-i 
with respect to x {s is fixed!) and the coefficients of the prolongation do not depend on f"^, 
a > s-l. □ 

In other words, if Q = rdt +(,dx + rjdu + O'^df; is a Lie invariance operator of system ([52]) and 
= for a fixed s and any a > s then the truncated operator Q = Tdt+S,9x+ri9u + Yl'^c;=i ^"^^f^ 
is a Lie invariance operator of the s-level potential system 

which is a subsystem of system (|52p extended by its differential consequences. This means that 
the corresponding potential symmetry of the initial equation is also induced by the truncated 
operator Q and, therefore, can be assumed to have order less than p. 

We also need to prove the statement on generalized potential symmetries of equations from 
class ([T]), which is similar to Lemma [211 This finally justifies the choice of conserved vectors in 
the canonical form (jlSp for the construction of potential systems. 

Lemma 22. Let system (j52p he a p-level potential system of an equation from class ([T|). Up to 
the equivalence of generalized symmetries, every generalized symmetry operator of system (|52p 
is obtained via the p-th order prolongation, with respect to only the variable x, of a generalized 
symmetry operator of the corresponding p-level modified potential equation (j53p and expressing 
the derivatives of f^ by f^ and derivatives of u according to system (j52p . In particular, up to the 
equivalence of generalized symmetries, the coefficients of every generalized symmetry operator of 
system (152]) depends at most on t, x, f^ and derivatives of u with respect to x and are linear 
with respect to the dependent variables and the derivatives of u. 

Proof. Suppose that Q = rjdu+9^dfs is a generalized symmetry operator of system (I52p . Here the 
coefficients 9^ := rj and 9^ are functions of t, x and derivatives of u and f^. Due to system ([52]) . 
its differential consequences and the equivalence relation of generalized symmetries, we can 
exclude the derivatives of of nonzero orders and derivatives of u containing differentiations 
with respect to t from the coefficients of Q. Therefore, they are assumed to depend, at most, 
on t, x, and derivatives of u with respect to x. 

We temporarily introduce the notations ■= d^u/dx^, k ^ 1, uq := f^ = u, Us '■= f^i 
ord9^ = maxjfc | d9^/duk ^ 0, k ^ —p}, f ■= oid9^, r ^ —p. The infinitesimal invariance 
condition |49i I51j applied to the equation /J = H^f^~^ implies the formula 

l-(9^^ + 9}.H'^f'^-' + Y,0luk+i). (54) 



k>0 



Iterating formula (|54p backward starting from s = p, we obtain the expression of 9^ ^ via 
the total derivatives of 9^ with respect to x up to order p — s + 1 according to the equation 
/^-i = {H')-^dx{{H'+^)-^dx{. . . {{HP)~^dxf'P) ...)). In particular, ord^-^ = r + p- s + l. 
In view of the other equations of system (|52p . the last equation = H^Afi ^ — f^~^ of 



the system is equivalent to equation (j53p . The equations on 9^ derived from the infinitesimal 
invariance condition for equation (|53p as a differential consequence of system (|52p coincide with 
the determining equations for the generalized symmetry operator Q' = 9^dfp of the single 
equation (|53p under the above relations between the ^-derivatives of and the functions "u^, 
k ^ —p. Therefore, there is a one-to-one correspondence between generalized symmetries of 
system (|52]) and equation ([53]) . established via the projection to the prolongation space over 
(t, X, fP) in the forward direction and the prolongation by formula (|54p in the reverse direction. 
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It is implicitly assumed that the x-derivatives of can be expressed in terms of the functions u^, 
k ^ —p, and vice versa. 

We can split the infinitesimal invariance condition for equation (j53p with respect to n^+i. 
(This is equivalent to splitting with respect to /p+^+i when (f53l) is considered as a single equa- 
tion.) Collecting the coefficient of u^+i, we derive the equation OZ^ur = 0. In view of the last 
equation, we can collect the coefficient of Ur+iUr, yielding OZrUr-i = 0. Iterating the procedure, 
in the l-th step {I ^ r + p + 1) we can collect the coefficient of Ur+iUr+2-i and obtain the 
equation 9?i^u,.+i-i = 0. The above procedure can be repeated for the terms containing Ur and 
so on. Finally we derive the equations 9uf,uy = 0, —p ^ k,k' ^ r. Then in view of the iterative 
formula (f54j) we have Ot^uy = 0, —p ^ k,k' ^ r + p — s + 1. □ 

Corollary 30. To exhaustively investigate potential symmetries of equations from class 

it suffices to only consider tuples of conserved vectors of the form (llSp . which are canonical 

representatives of the corresponding conservation laws. 

Proof. Consider p conserved vectors of an equation of the form ([T]), corresponding to p linearly 
independent conservation laws. They necessarily possess the representation 



where a** = a^(t, x), s = 1, . . . ,p, are linearly independent solutions of the adjoint equation ([18]) 
and <I>* are functions of t, x and derivatives of u. Due to the initial equations we can assume that 
the derivatives are only with respect to x. The corresponding potential systems are obtained 
by the substitution = ij^ — from the system (f50|) associated with the tuple of equivalent 
conserved vectors in the canonical form. Here are the potentials generated by the conserved 
vectors (I55p . In view of Lemmas 1161 and [TSl the induced substitution in terms of the potentials 
has the form = — For each value of s the function is a linear combination of 
a = 1, . . . , s, with coefficients depending on t and x, and the coefficient of does not vanish. 
Denote the maximal order of derivatives in by p and a value of s with ^'^^ 7^ by sq- We 
assume p ^ 1 since otherwise the functions and, therefore, the functions <I>* can be neglected 
due to the point transformation, which has, up to similarity, no influence on Lie symmetries of 
the involved potential systems. 

Every Lie symmetry of the system in and u corresponds to a generalized symmetry of 
system (I52p in and u. The question is when a generalized symmetry operator Q = r]du + 9^dfs 
of system ()52p induces a Lie symmetry operator Q = fjdu+O^dfs (in evolution form) of the system 
in and u. 

We use the notations of Lemma [22l Suppose that r + p ^ 1, where r = ord^^. Since 



the term rju^^^^l^^^Ur+p+p cannot be canceled with other terms of 0**". Therefore, ord^*" = r + 
p + p ^ 2. At the same time, if Q would be a Lie symmetry operator, ord 6^ ^ max(2 — s, 0) ^ 1. 

In the case r + p = the coefficient 6^ has the form = eP^{t,x)fP + 6lP°(t,x). The 
determining equations for the generalized symmetry of (j52p imply that 9P^{t,x) = C = const. 
Then 9'-'^ = Cf'-"" + 9'-^^^{t,x) in view of formula ([Ml)- Since fj = = Cu + 9^^{t,x), the 
operator Q gives a trivial potential symmetry of the initial equation, which corresponds to a 
trivial Lie symmetry of this equation. □ 

In view of Lemma [211 Lie symmetry analysis of any p-level potential system associated with an 
equation from class ([T]) is reduced to similar investigation for the corresponding p-level potential 



( 



•x 
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equation. In fact, we will investigate the modified level potential equation \wP\ instead of \f^\ ■ 
Such way has two advantages. Firstly, the potential does not depend on nondegenerate linear 
combining (including a rearrangement, cf. Corollary [26]) in the characteristic tuple. Secondly, 
the modified p-level potential equation is connected with the initial one via the multiple Darboux 
transformation expressed in terms of the corresponding characteristics. In terms of w's the p- 
level potential system has the form 



which will be called the modified p-level potential system associated with an equation from 
class ([I]) and the characteristic tuple (a^, . . . ,a^). 

Analogously to the case of simplest potential symmetries, the only possibility of obtaining 
pure potential symmetries is connected with the coefficient t]. Namely, the condition rjfs ^ 
(or 77^3 7^ in terms of w^s) should be satisfied for some s. In fact, we are interested in 
the investigation p-order potential symmetries only in the case when it is not reduced to the 
consideration of potential symmetries of a smaller order. The irreducibility is naturally defined 
in terms of the w^s. One of the reasons for this again is the independence of of nondegenerate 
linear combining in the characteristic tuple. Another reason is the following. Let Q = rdt + 
^dx + rjdu + 0''df^ and Q = Tdt+£,dx + r]du + C'dw'i be the representations of the same operator Q 
in terms of /'s and w^s. In view of Corollary [27] the conditions '^ja ^ = for a fixed value 
of s and any a > s — V and '^^^^ = for any <^ ^ s and any a > s — V are equivalent to the 
similar conditions in terms of the w^s, i.e., = for a fixed value of s and any a > s — V 

and 'C^^^ = for any ? ^ s and any a > s — 1\ Therefore, Corollary [29] can be completely 
reformulated in terms of the w's. 

Corollary 30'. // Cw°"^ ~ for a fixed value of s and any a > s — 1 then ^^^"^ = for any 
^ s and any a > s — 1. 

Definition 11. Let Q' = rdt + £,dx + C^d^p be a Lie invariance operator of the modified p-level 
potential equation associated with an equation from class ([1]) and the characteristic subspace 
(a^, . . . , a^). We will say that Q' generates a strictly p-th order potential symmetry of the initial 
equation if for any basis (d^, . . . , dP) in (a^, . . . , a^) the prolongation Q = Q' + C,^~^d^s-i of Q' 
to the corresponding potentials where := u, satisfies the following condition. For any s 

there exists a > s — 1 such that C,^~^^ 7^ 0. 

Definition 12. Let 5 be a Lie invariance algebra of the modified p-level potential equation 
associated with an equation from class ([1]) and the characteristic subspace {a^, . . . , a^). We will 
say that g generates a strictly p-th order potential symmetry algebra of the initial equation if 
for any basis (a^, . . . , qP) in {a^, . . . , a^) and for any s there exists Q' € Q whose prolongation 
Q = Q' + C*~^c?u,s-i to the corresponding potentials w'*"^, where := n, satisfies the following 
condition. There exists a > s — 1 such that C,^~^ ^ 0. 

Roughly speaking, a potential symmetry operator (resp. algebra) is strictly of order p if it 
cannot be obtained from a smaller number of conservation laws and potentials. 



Let us recall that = u and C,^ = rjhy definition. It follows from Lemma [21] in view of the 
formula = f^/l3^~^'^ that 



w = w 




(56) 



p 



p 



r = EC"(i,xV'^ + x) 
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and a more precise form of the coefficients ^ is calculated by an above backward recursive 
formula involving r, ^ and C^, in accordance with equations of (j56p . In particular, 

1 W(wP'P O'lujP'P]) 

1 Wiw'P'''^^ ujP'P O'lwP'P]) 

w'P'P ' wP'PW{'wP'''+^,. . . ,wP'P) 

by the Crum theorem. Here Q'[wP'P] = QPPvjP'P — tw^'^ — ^Wx'^. 

Consider the prolongation of transformations from the equivalence group G~ to the general 
potential frame. Any point equivalence transformation T in class ([T]) acts on the variables 
{t,x,u) by the formulas i = T{t), x = X{t,x), u = U^{t,x)u, where TtXxU^ / 0. In view 
of Proposition [7] it is prolonged to characteristics a* of conservation laws of an equation from 
this class: d'' = /{XxU^). The corresponding potential v'^ can be assumed to be transformed 
identically under prolongation of T. (See Section [71) Taking into account the constructed 
representations for Z?*'"", g^'^ , a ^ s, and w'^ via and f we obtain the following statements. 

Lemma 23. For any p £ N any transformation T G G~ is prolonged to the p-order potential 
frame over class The prolonged transformation TP 

a" 

t = T{t), X = X{t, x), u = U^{t, x)u, a"" = , = ^;^ 

A = -2iA, B = { B -2--^A] ^— — —, C = -—L—-, 

where TtXxU^ 7^ 0, realizes a simultaneous transformation between the tuples consisting of the 
initial, potential and modified potential equations and the potential systems in terms ofv^ and f^ . 
The transformations from , prolonged to the p-order potential frame over the class form 
the group Gj^j called the equivalence group of this frame. 

Note 24. In general, under the transformation T: t = T[t), x = X{t,x), ip = ^{t,x)'ip, 
ip'^ = ^{t,x)'ip^ we have 

Note 25. The prolonged transformations from C do not exhaust all possible equivalence 
transformations of the p-order potential frame over the class ([T]). They can be extended, e.g., 
with linear combining of characteristics. If the variables (t, x, u) (and, therefore, the arbitrary 
elements {A,B,C)) are not transformed and = a'^Cas, where Cas = const, det(co-s) / and 
Ccts = 0, (t > s, then the corresponding transformation of the other functions appearing in the 
potential frame is easily constructed: 

V'=V''C^S, r'"=r''c,., r"=/''c,., f = fCss, W'=W', W^^'^ =W'^'C,,. 

Here (c^o-) is the inverse matrix to (c^o-)- 

The transformations from GjJ^j preserve the determining equations derived in the proof of 
Lemma [211 Let Q = rdt + ^dx + rjdu + O^dfs be a Lie invariance operator of system (j52p . The 
coefficients of Q are transformed under the operator mapping generated by TP £ GJ^j by the 
formulas 



f = rTt, i = rXt + ^Xx, r? = tU^u + ^U^u + U'rj, 
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Therefore, both the conditions rjfcr ^ and ijfcr = as well as the conditions 9ja. ^ and 
9jrT = are preserved by the transformations from Gj^] for any s and a. This means that 
pure p-order potential symmetries of any equation from class ([1]) are not mixed under the 
transformations from GjJJ] with either potential symmetries of lesser orders or Lie symmetries. 
The dimension of the factor-space of p-order potential symmetry operators corresponding to a 
p-tuple of characteristics with respect to the subspaces of potential symmetry operators of lesser 
orders are also not changed. Therefore, p-order potential symmetries of equations from class ([T]) 
can be studied up to the equivalence relation generated by transformations from G'^y 

There are different ways of employing this equivalence relation. One of them is to emphasize 
the simplification of the form of the equations under consideration. In particular, we can put 
A = 1 and B = (and re-denote C by —V) in ([T|), which implies = B — sA^ = 0. As a result, 
the symmetry analysis of the p-order potential frame over class ([1]) is reduced to the symmetry 
analysis of the p-order potential frame over class ([7|). The equivalence group of the reduced 
class ([7]) is canonically isomorphic to a subgroup of the equivalence group G~ with special 
restrictions of the parameter- functions X and U^. In view of Lemma [23] the transformations 
from Gi are prolonged to equivalence transformations of the whole reduced potential frame. The 
prolonged transformations form the group Gj^j which is canonically isomorphic via projection 
to the equivalence group of the class of modified p-level potential equations for the equations of 
the form ([Tj). Hence, the classification of potential symmetries of the class ([7]) follows from the 
group classification of the same class in terms of {w^, V^) instead of (u, V). 

Another way to proceed is to simplify the form of the operators. A Lie invariance operator 
Q' of an equation from class ^ with a nonvanishing coefficient of dt (or a vanishing coefficient 
of dt and a nonvanishing coefficient of dx) can be reduced by transformations from G~ to the 
form Q' = df [Q' = dx)- Note that an appropriate simplification of the form of the equations 
leads to a simplification of the form of their symmetry operators and vice versa. The choice of 
how to employ the above equivalence relation depends on the problems under consideration. 

At first, we formulate a criterion when a Lie invariance operator of the modified p-level 
potential equation generates a strictly p-th order potential symmetry of the initial equation. A 
simplification of the form of symmetry operators will be very helpful in this problem. 

Theorem 9. Suppose that \wP\ is an equation from class ip'^ , ip^ are its linearly 
independent solutions and DT['(/'-'^, . . . , -0*'] \wP\ = \u\ . 

1) Trivial Lie invariance operators of \w'^\ generate only trivial Lie invariance operators 
of \u\. Namely, w^d^p generates udu and ip{t, x)dwP generates DT[tlj^ ,iljP]{ip)du- Here the 
function tp = ip{t,x) runs through the solution set of \wP\ and, therefore, DT[ip^ , . . . , ipP]{tp) 
runs through the solution set of \u\ . 

2) Let Q' he an essential Lie invariance operator of [w^J . Then Q' generates a strictly p-th 
order potential symmetry of \u\ iff any subspace of {tp^, . . . , tp^) is not invariant under the action 
of the associated first-order differential operator Q' . 

Proof. Item 1 is obvious in view of the properties of the Darboux transformation. It is significant 
only complement to item 2. 

To prove item 2 we show the equivalence of the respective negations. Namely, Q' does not 
generate a strictly p-th order potential symmetry of [nj iff there exists a subspace of {^p^, . . . , ■0^) 
which is invariant under the action of the associated first-order differential operator Q' . 

Let us recall that any (nontrivial) essential Lie invariance operator Q' of \w^\ has the form 

Q' = rdt + ^dx + e'wPd^r>, 

where r = T{t), ^ = S,{t,x) and = C^^{t,x) are smooth functions of their arguments and 
(tjC) / (0)0). The first-order differential operator Q' associated with Q' acts on functions 
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of t and X according to the formula Q'tp = Q'[ip] = C^^ — Tipt — CV'x- Due to the equivalence 
relation generated by transformations from Gjp] we can assume without loss of generality that 
Q' G {dt,dx}- Then the standard prolongation Q'^^^ of Q' of any order q formally coincides 
with Q'. 

Suppose that a g-dimensional subspace T in the linear space (■0^, . . . , V'^) is invariant under 
action of Q' . Let s = p — q. Without loss of generality we can choose a basis {w^'^, . . . ,it;^'^} 
in {ip^ , . . . , ipP) in such a way that the functions wP'^~^^, . . . , w^'^ form a basis of Z. Therefore, 

p 

Q'K"] = Yl '"-^'^"''^ a = s + l,...,p, 

<;=s+l 



where k^c;, o", ? = s + are constants. Since = DT[wP'^~^^ , . . . ,wP''^]{wP) then the 



coefficient of d^a in the corresponding Lie invariance operator Q of system ([56|) is = Q', .w^ ^ 



I.e., 



{k^o- - K„a)w'' = 0. 



■W 



a=s+l 
P 



a=s+l 



Corollary [30j implies = for any ? ^ s and any a > s. This means that Q' does not generate 
a strictly p-th. order potential symmetry of \u\ . 

Conversely, suppose that Q' does not generate a strictly p-th order potential symmetry of [nj . 
Therefore, there exists a basis {ujP'^, . . . ,wP'P} in {ip^, . . . jipP) such that the corresponding Lie 
invariance operator Q = Q' + C'^~^d^^-i of system ([56|) satisfies, for some fixed value s, the 
condition ^^^^ = for any ? ^ s and any a > s — 1. In particular, 



1 ^ W{wP''+\...,wP'P,Q'lwP'P]) 
^"^^ ~ wP'PW{wP'''+'^,...,wP'P) ~ ' 



i.e., Q'lujP'P] E {wP'^^^, . . . ,wP'P). In view of Corollary [27] the condition = for any ^ ^ s 
and any o" > s — 1 is equivalent to C^7^ = for any <; ^ s and any a > s — 1. Thus the 
order in the tuple (u'^'*"'"^, . . . ,wP'P) is inessential for the condition on the C's and, therefore, 
Q'lwP'''] G (wP'^+i, . . . , wP'P), a = s + 1,... ,p. It follows that the subspace {wP''+'^, wP'P) of 
the space {ip'^, . . . , ipP) is invariant under the action of Q'. □ 

Corollary 31. An essential Lie invariance operator Q' of \wP\ generates a strictly p-th or- 
der potential symmetry of \u\ = DTlip'^, . . . ,■0^] \wP\ iff any one- or two-dimensional subspace 
of {tp^, . . . jtpP) is not invariant under the action of the associated first-order differential opera- 
tor Q'. 

Proof. In view of Theorem [9l it is enough to prove that a finite-dimensional space I invariant 
under the action of Q' has a one- or two-dimensional invariant subspace. Let q = dimX and 
if^, ip'^ form a basis of X. Then = K^^ip'^, where k^,; are constants. Hereafter the 

indices a and <j run from 1 to If (ci, . . . , Cq) is an eigenvector of the matrix (ko-^ ) then {ca^'^) is 
a one-dimensional invariant subspace of T. In the real case, if the matrix {ku^) has no real eigen- 
values, we can take an eigenvector (ci, . . . , Cg) of the complexification. Then {ip'^ Re c^, f'^ Im c^) 
is a two-dimensional invariant subspace of I without proper invariant subspaces. It is obvious 
that in the complex case it is enough to consider only one-dimensional subspaces. □ 
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Note 26. Similarly to Propositionll21 Theorem [9] can be re-formulated in different terms. Thus, 
a subspace {w^'^ , . . . ,wP''^) of the solution space of \wP\ is invariant under the action of the 
operator Q' and Q'[wP''^] = k^^w^''' iff (w^'^, . . . is an invariant solution of the uncoupled 

system of q copies of \wP\ with respect to the operator 

Q' = Tdt + id^ + C^wP^'^d^P,. - K^,wP''d^P,.. 

Here constants. The indices a and ? again run from 1 to q. 

The characterization of linear second-order parabolic equations possessing potential symme- 
try algebras of a fixed order is given by the following theorem. 

Theorem 10. A linear (1 + 1)- dimensional second-order parabolic equation admits a strictly 
p-th order potential symmetry algebra iff it is equivalent with respect to point equivalence trans- 
formations to an equation from class ([7]) in which 

where ip^ = ip''{t,x) are linearly independent solutions of the equation tpt — ipxx-'r P{x)'4^ = and 
either P = ^ = const, or no subspaces of {ip'^, . . . , ip^) are invariant under the action of dt 

if P is inequivalent to fj,x~'^ with respect to point equivalence transformations. The associated 
characteristic subspace is {a^, . . . , a^) , where 

In the case P = = 0) the potential symmetry algebra contains at least one (two) inde- 

pendent operators which are linearly independent up to Lie symmetries and essentially involve 
the p-th order potential. For the general value P = P{x) a sufficient condition for the potential 
symmetry algebra corresponding to {ip^, . . . .,ipP) having strictly p-th order is the following. For 
an arbitrary choice of basis in {ip^ , . . . ,-0^) all the modified potential equations \w^~^\ including 
\w^\ = \u\ have non- stationary functions as values of the arbitrary elements V^'"^. 

Proof. The equivalence relation generated by transformations from GJJ^j reduces the p-th order 
potential frame over class ([T|) to the p-th order potential frame over class ([7]) in which there is an 
analogous equivalence relation generated by transformations from Gj^] . In view of Lemma [2H 
the investigation of p-th order potential symmetries can be replaced by the investigation of Lie 
symmetries of p-th order potential equations. The projection of the group GjJ] to {t,x,wP , P), 
where P := V^, coincides with the equivalence group of the class of modified p-level potential 
equations for the equations of the form ([7]). This class is a copy of the class d?]) written in 
terms of {t, XjW^ , P), i.e., both these classes have the same (up to the notation of variables, 
arbitrary elements and group parameters) equivalence group. Therefore, the Lie-Ovsiannikov 
classification described in Theorem [2] can be used here. Since the initial equations of form ([7]) 
and the corresponding modified p-level potential equations are connected via the multiple Dar- 
boux transformation, the Crum theorem [20^ [46] provides a connection between their arbitrary 
elements. To complete the proof, we have to examine all the inequivalent cases of extension of 
Lie symmetries from Theorem [21 re- written in terms of {t,x,wP,P). The general case without 
extension is eliminated in view of item 1 of Theorem [9l 

Let P = Ijlx~'^, i.e., the modified p-level potential equation has the form w^ — Wxx+^J'X~'^w'P = 0. 
We prove by contradiction that there are no finite dimensional subspaces in the solution space 
of [u;^J which are simultaneously invariant under the action of the differential operators dt, D 
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and n corresponding to dt, D and 11. Suppose that such a subspace I exists. Let q = dimX 
and , . . . , if'^ form a basis of X. Then 

n[v9-] = -At\l - Atxipl - {x^ + t)ip^ = 

where /t^^, and are constants. In what follows the indices a and <; run from 1 to q. The 
above equality implies that {K'^'' — x'^5a<;)'p'" = 0, where (Jo-? is the Kronecker delta and 

For any fixed t, the determinant of the matrix {K'^'' — x'^d^q) vanishes only for a finite number 
of values of x. Therefore, for any fixed t ip''{t,x) does not vanish for at most a finite number of 
values of x, i.e., by continuity all ip'^ are equal to identically. This implies a contradiction. 

In the case /i = we additionally have the operators dx and G. They also possesses no simul- 
taneously invariant finite dimensional subspaces in the solution space of \wP\ since otherwise by 
an analogous reasoning we have 

where and are constants. Then (k^^ — K^^t — x6o-<;)(p'^ = and identically, which 
implies the same contradiction. 

There is only one independent nontrivial Lie symmetry operator dt for the general value 
P = P{x). If it does not induce a strictly p-th order potential symmetry operator of ^ then it 
has an invariant subspace I in {ip^, . . . , tp^). Let q = dimX and (^^, . . . , ip'^ form a basis of I. 
Then ifl = K^^if'' for some constants Ka<; and, therefore, for the modified {p — q)-level potential 
equation obtained from \w^\ with DT[(/9"'^, . . . , ip'^] we have 

i.e., VP~'^ = V^~''{x). As a result, the problem for the p-order potential frame is completely 
reduced to the same problem for the {p — g)-order potential frame. □ 

Based on Theorem 1101 we can formulate pure symmetry criteria on the existence of potential 
symmetries of arbitrary order without involving equivalence transformations. 

Corollary 32. A linear second-order parabolic equation admits the strictly p-th order potential 
symmetry algebra associated with its p-dimensional characteristic subspace (q^, . . . only if 
the corresponding p-level potential equation possesses nontrivial Lie symmetry operators. If the 
p-level potential equation has more than one independent nontrivial Lie symmetry operators then 
the potential symmetry algebra is of strictly p-th order. More precisely, if the p-level potential 
equation has more than one (three) independent nontrivial Lie symmetry operators then for 
any choice of basis in {a^,... ,af) the potential symmetry algebra contains at least one (two) 
independent operators which essentially involve the p-th order potential. 
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11 On number and order of potential symmetries 



The title of the section is shghtly vague. Our aim here is to investigate the following basic 
questions about potential symmetries of linear parabolic equations: Can a fixed equation from 
class (dl) possess an infinite series of potential symmetry algebras? How can equations having 
potential symmetry algebras of all orders be constructed? What orders of potential symmetries 
are possible for a fixed equation? When are the orders of potential symmetries bounded? We 
study only certain examples that, nevertheless, allow us to formulate quite general statements 
answering some of these questions. 

It is easy to construct equations having potential symmetry algebras of all orders. Indeed, we 
take the linear heat equation wt = Wxx as a potential equation (i.e., P = 0) and choose, for any 
p G N, the p-tuple "0 = {Pq, . . . , Pp-i) of its solutions. Here is the canonical heat polynomial 
of degree A;, 

P2m{t,x) = + 1! (2m - 2)! ^ 2! (2m - 4)! ^ ■ ■ ■ ^ (m - 1)! ¥ ^ ^' 

™2m+l J. „2m— 1 j.2 ™2m— 3 j.m—1 „3 j-m „ 

2m+i{t,x) = (2m + 1)! ^ 1! (2m - 1)! ^ 2! (2m - 3)! ^ ■ ■ ■ ^ (m - 1)! ¥ ^ ^ 1! ' 

A;,m € NU {0}. Note that dP^+i/dx = Pk hence d^Pk/dx^ = 1, d^Pk+i/dx^ = x. A direct 
calculation implies that W{Pq, . . . , Pp-i) = 1 and that the corresponding value of the arbitrary 
element V equals for any p G N. Therefore, DT[Po7 • • • , Pp-i] \w\ = \u\ , where \u\ also is the 
linear heat equation ut = u^x- In fact the multiple Darboux transformation DT[Po) • • • > -fp-i] is 
nothing but p-order differentiation with respect to x: u = DT[Pq, . . . , Pp-i] w = d^w/dx^. Let 
(a^^, . . . , a^P) be the characteristic tuple of \u\ , dual to the solution tuple [Pq, . . . , Pp-i) of \w\ , 
i.e., 

aP^ = {-iy'^W{Pi,...,Pp.s), s = l,...,p-l, aPP = {-ir-\ 

(See Corollary [281) The Wronskians W = W{Pi, . . . , Pg), g E N, and := 1 are solutions of 
the backward heat equation and additionally satisfy the conditions 

dWydx = W'^-^, H^'?(0,0)=0. 

Therefore, = Pq{—t,x) is the backward heat polynomial of order q and 

aP^ = (-ir-iPp_,(-t,x), s=l,...,p. 

In view of Theorems [8l and \T0\ for any p G N the potential symmetry algebra Qp of \u\ , associated 
with the p-dimensional characteristic subspace {a^^, . . . , a^P), is of strictly j?-th order. For any 
choice of basis in (a^^, . . . , a^P) the potential symmetry algebra contains at least two independent 
operators which essentially involve the p-th order potential. Summarizing these results, we can 
formulate the following statement. 

Proposition 13. The linear heat equation admits an infinite series {Qp, p G N} o/ potential 
symmetry algebras. For any p G N the algebra Qp is of strictly p-th potential order and is 
associated with p-tuples of the linearly independent lowest order polynomial solutions of the 
backward heat equation. Moreover, it is the standard p-th prolongation, with respect to only x, 
of the Lie invariance algebra Qq, re-written in terms of{t,x,w) and, hence, is isomorphic to go- 
Other examples can be constructed in a similar way. Thus, for the same potential equa- 
tion and the p-tuple ip = (Pq, Pi, . . . , Pp_2, Pp) we have W{Po, Pi, ... , Pp-2, Pp) = x and the 
corresponding value of the arbitrary element V is equal to 2x~^ for any p G N. Therefore, 
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T)T[Pq, Pi, ... , Pp_2, -Pp] \w\ = \u\ , where \u\ denotes the equation ut — u^x + 2x~^n = 0. Let 
(a^^, . . . , aPP) be the characteristic tuple of \u\ dual to the solution tuple {Pq, Pi, ... , Pp-2, Pp) 
of \w\ , i.e., 

^ps ^ . . . , s = 1, . . . ,p - 1, oPP = {-lY-^X-\ 

The Wronskians = W{Pi, . . . ,Pq-2,Pq), q ^ 2, satisfy the conditions dW^dx = xW'i-'^, 
W'^{0,0) = 0. The quotients W'^/x are solutions of the adjoint equation vt + Vxx — 2x~'^v = 0. 
Therefore, W'^ = xWi - W and 

^ps ^ ^_^y-l^^p-s _ ^-l^p-s+l^^ s = 1, . . . ,p - 1. 

Here denotes the backward heat polynomial Pq{—t,x) of order q. Analogously to the pre- 
vious example, for any p G N the potential symmetry algebra 0p of \u\ associated with the 
p-dimensional characteristic subspace (a^^, . . . , oP^), is of strictly p-th. order. For any choice of 
basis in {a^^ , . . . , a^^) the potential symmetry algebra contains at least two independent op- 
erators which essentially involve the p-th order potential. As a result, the equation [uj has 
an infinite series {flp, p G N} of potential symmetry algebras isomorphic to the Lie invariance 
algebra go of the linear heat equation. 

The described construction can be generalized. Let [uj be the equation ut — Uxx + Vu = 0, 
where the function V has the form (j57p with P = 0, i.e., tp^ = ip^{t,x) are linearly independent 
solutions of the linear heat equation -0^ = ^l^xx■ This means that \u\ is the image of the linear 
heat equation \w\ under the Darboux transformation DT['0-'^, . . . , ip^]. Then for any q £ N 

\u\=BT[Po,...,Pq-i,r\...,rP] \w\, 

where ip'^^ for any fixed s is a solution of the linear heat equation, being a preimage of under 
the Darboux transformation DT[Po, . . . , Pq-i], i.e., DT[Po, . . . , Pq-l]^p'^^ = tp^ . The function 'p'^^ 
is found by g-fold integration of the function 'p^ with respect to x with a special choice of an 
'integration constant' depending on t. The best way is to employ of the recursive formula 
'p'x = ip'^~^'^, 'Pi^ = ip^^ := ip^. In view of Proposition 1131 this implies the following 

statement. 

Corollary 33. Suppose that a linear (1 -|- 1)- dimensional second-order parabolic equation is 
equivalent with respect to point equivalence transformations to an equation from class ([7]) in 
which 

y^_, ( {w{^p\...,r))A 

[ W{i;\...,i^P) jj 

where 'p^ = tp'^{t, x) are linearly independent solutions of the linear heat equation -pt = 'pxx- Then 
this equation possesses an infinite series {ip+fe, /c G N U {0}} of potential symmetry algebras. 
Each algebra flp+fc from the series is of strictly {p+k)-th potential order, is isomorphic to the Lie 
invariance algebra Qq of the linear heat equation and contains at least two operators essentially 
involving potentials of the {p + k)-th level. 

It is easy to see that a basis of the above consideration is formed by the construction of the 
infinite series of the auto-Darboux transformations {DT[Po, . . . , Pp-i],p £ N} for the linear heat 
equation. Similar series of auto-Darboux transformations exist for all the equations of the form 

Ut — Uxx + fJ'X~'^u = 0. (58) 

To show this, it is sufficient to prove that equation (j58p has an infinite series of linearly inde- 
pendent solutions such that the Wronskian of p first solutions from the series is constant for an 
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infinite set of p's. We restrict ourself to solutions which are polynomials in t. Let tp^^, i = 1,2, 
be linearly independent stationary solutions of ([58]) . i.e., 

The following values of c^''* can be taken: 

^02 ^ |^|l/2|j^|^| •£ 1 + 4^ = 0, 





\x\ 


V— 


VpOi = 


\x\ 


1/2 


^'' = 


\x\ 


1/2 



cosKln|x|, (^''^ = sin/tin |a;|, k = y^— 1/4 — ^, if 1 + 4^ < 0. 

Consider the functions c/?^* = 11^ (/j*^*, where 11 = —4t'^dt — Atxdx — x"^ — 2t. They are linearly 
independent and polynomial in t. They are solutions of (I58p as a result of the action of the 
symmetry operator 11 on solutions. The Wronskian = W {^p^^ , Lp^"^ , . . . , Lp^^ , Lp^"^) does not 
depend on t since (^^* necessarily is a linear combination of the functions (p^^ , k' < k. Therefore, 
it is enough to evaluate W'^ for a single value of t. Putting t = 0, we prove by induction 
that = 0, i.e., W'' is a nonzero constant. As a result, BT[(p^^ , ip^"^ , . . . , ip''^ , ip'^'^] is an 
auto-Darboux transformation of equation (j58p for any G N U {0}. 

Proposition 14. Equation (jSSp admits an infinite series {Q2q, 9 G N} of potential symmetry 
algebras. For any g G N the algebra Q2q is of strictly 2q-th potential order and is associated with 
2q-tuples of the linearly independent solutions of (j58p which are lowest order polynomials in t. 
Moreover, it is isomorphic to the Lie invariance algebra go of equation (jSSp and contains at least 
one operator essentially involving potentials of the 2q-th level. 

Combining Darboux, auto-Darboux and equivalence transformations similarly to the proof 
of Corollary [33l we derive the following corollary of Proposition [TH 

Corollary 34. Suppose that a linear (1 + 1)- dimensional second-order parabolic equation is 
equivalent with respect to point equivalence transformations to an equation from class ([7]) in 
which 



[ T^(v^...,v^) 



where ip^ = tp^{t,x) are linearly independent solutions of the equation tpt — ipxx + fJ'X~'^tp = 0. 
Then this equation possesses an infinite series {0p+2fc5 k €^ N U {0}} of potential symmetry 
algebras. Each algebra Qp-\-2k from the series is of strictly {p+2k)-th potential order, is isomorphic 
to the Lie invariance algebra Qq of equation (j58p and contains at least one operator essentially 
involving potentials of the {p + 2k) -th level. 

Note 27. An equation from class ([7]) with a stationary value of the arbitrary element V has 
auto-Darboux transformations constructed with solutions polynomial in t. In contrast to the 
special case V = /ix'^, this does not imply conclusions in the general case V = V{x) since, in 
particular, Theorem 1 1 1 does not give a sufficiently powerful criterion on the existence of potential 
symmetries for such a situation. 

12 Discussion 

In the present paper we investigate symmetries and conservation laws of linear (l-l-l)-dimensional 
second-order parabolic equations. In our opinion, the most important results of the paper are 
the following. 
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• It is proved that any potential conserved vector of a linear parabolic equation is equivalent 
to a local one (Theorem [5]). The local conservation laws are described in Theorem [H 
Namely, the space of characteristics associated with the local conservation laws of a linear 
parabolic equation can be identified with the space of the functions depending only on t 
and X which are solutions of the adjoint equation. 

• Different criteria on the existence of potential symmetries of arbitrary order are formulated. 
Theorem [9] gives a criterion on operators of potential equations, generating potential sym- 
metries with the same orders as the level numbers of the potential equations. The shape 
of the equations having potential symmetries of a strongly fixed order is characterized in 
Theorem [101 Corollary [32] supplies a pure symmetry criterion in terms of the number of 
independent nontrivial Lie symmetries of potential equations. 

Extensive preparatory considerations were needed for obtaining these results. Thus, the group 
classification of class ([T]) has a twofold application for the investigation of potential symmetries 
in this class. Firstly, an exhaustive knowledge on Lie symmetries is necessary for focusing our 
attention on pure potential symmetries. At the same time, the modified potential equations of 
any fixed level for the equations from class ([1]) form a copy of this class. Therefore, results on the 
group classification in class ([1]) also are directly used in the description of potential symmetries. 
This entails the need for a careful revision of these results. The study of normalization properties 
of class ([1]), their subclasses and the associated classes of inhomogeneous equations justifies the 
choice of the gauge A = 1, B = under both the group classification of single equations 
and the analysis of the whole potential frame over the class ([1]). It also gives a well-founded 
explanation of the difficulties arising in the classification of the Kolmogorov and Fokker-Plank 
equations. 

The investigation of local conservation laws of equations from class ([T]) naturally leads to 
the consideration of systems of two mutually adjoint equations from class ([TJ. Such systems 
form the basis of the potential frame over class ([T]) and are directly connected with the so-called 
adjoint variational principle. The framework of the adjoint variational principle is extended 
to admissible transformations. A number of auxiliary statements on a hierarchy of normalized 
classes of second-order evolution systems are proved. The group classification of Fokker-Plank 
equations is obtained from the group classification of the Kolmogorov equations in a simple way 
based on the application of the adjoint variational principle. 

A feeling for the general problem, understanding possible ways of solving it and a rough 
shape of formulating the final results arise during the consideration of simplest potential sym- 
metries. The techniques resulting from this approach involve different ideas. Thus, the po- 
tential systems should be simultaneously studied with the associated initial, adjoint, potential, 
modified potential and adjoint modified potential equations which together form the potential 
frame over the class ([T]) (of the first order in the case of simplest potential symmetries). For 
applications of this idea to general potential symmetries, an analogue of the potential equa- 
tion, corresponding to a tuple of characteristics, should be proposed at first. The equivalence 
group of the initial class ([1]) is prolonged to the whole potential frame including characteris- 
tics of the initial equations, potentials, modified potentials and characteristics of the modified 
potential equations. It is shown that the problem can be investigated up to the equivalence 
relation generated by the prolonged equivalence group. Moreover, the classification of Lie sym- 
metries of potential systems with respect to the above equivalence relation is reduced to the 
same classification for modified potential equations. This allow us to use the Lie-Ovsiannikov 
classification. A connection between different objects of the potential frame is provided via the 
dual Darboux transformation which, for this reason, is an important component of the proposed 
technique. 
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The study of general potential symmetries requires a development of the above ideas in re- 
lation to simplest potential symmetries and the creation of specific tools. It appears that the 
p-order and p-level potential systems associated with the same tuple of characteristics should be 
investigated simultaneously. (This conclusion cannot be drawn in the consideration of simplest 
potential symmetries since the corresponding 1-order and 1-level potential systems coincide.) 
The potential equations associated with characteristic tuples are introduced in a natural way 
due to the special iterative procedure for the construction of potential systems with levels higher 
than 1. Moreover, the modified potential equations are in fact associated with characteristic 
subspaces spanned by the corresponding characteristic tuple. This allows us to study potential 
symmetries up to forming linear combinations of elements of a characteristic tuple. Enhanc- 
ing results on the multiple dual Darboux transformation between equations from the class ([T]), 
presented in Theorem [U is of fundamental importance for both the consolidation of the whole 
potential frame and for deriving definitive statements on potential symmetries. The proof of 
Lemma [2T] on Lie symmetries of potential systems is quite intricate and involves a number of 
tricks. It is essentially based on the higher-level representation of potential systems. The next 
important step is the prolongation of the equivalence group of the class ([I]) to the whole poten- 
tial frame in Lemma [23j The confluence of these three components (the known shape of the 
Lie symmetries of potential systems, the multiple dual Darboux transformation and the pro- 
longed equivalence transformations) results in criteria on the existence of potential symmetries, 
formulated in Theorems [9] and [10] and Corollary [32l as well as in subsequent estimations of the 
number of potential symmetries for subclasses of the class ([T]). 

Another new problem on potential symmetries, which was first posed in general and solved for 
the (1 + l)-dimensional linear parabolic equations in the present paper, is to justify the choice 
of natural representatives among the equivalent conserved vectors for introducing potentials. 
Different choices of representatives are equivalent only with respect to generalized potential 
symmetries of arbitrarily high order. The construction of the natural potential frame for the 
class ([1]) was justified by Corollary [30l To study (usual) potential symmetries of equations 
from the class ([T]), it is in fact sufficient to only consider conserved vectors of the canonical 
form ()18p . which have the lowest orders in both the flux and the density. Any tuple of conserved 
vectors of any equation from the class ([T]) , containing a conserved vector of a higher than lowest 
order, gives only symmetries trivial in all senses. In contrast to the case of simplest potential 
symmetries, the proof of the above statements for an arbitrary number of potentials needed a 
description of the generalized potential symmetries the equation under consideration. As shown 
in Lemma [22l each generalized symmetry of every potential system associated with conserved 
vectors in canonical form is linear in the dependent variables and their derivatives up to standard 
equivalence of generalized symmetries. Moreover, there is a one-to-one correspondence between 
generalized symmetries of a potential system and those of the associated potential equation. 

In spite of the wide range of the performed investigations, a number of interesting and 
difficult problems concerning potential symmetries of linear (1 -|- l)-dimensional second-order 
parabolic equations remains unsolved. We list some of them, without making any claim on 
completeness. 

Problem 1. What is the maximal number of inequivalent characteristics leading to simplest 
pure potential symmetries for a fixed equation from class ([I])? Classify such characteristics. 
A similar question can be asked about p-tuples of linearly independent characteristics giving 
strictly p-order potential symmetries. 

In Section [8] Problem [1] is solved only for simplest potential symmetries of the linear heat 
equation. The answer to the question about the number of characteristics in this case is two. 
More precisely, any appropriate characteristic is equivalent to 1 or x with respect to the essential 
part of the point symmetry group of the linear heat equation. 
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Problem 2. Given a fixed equation from class ([T]), are strict orders of its potential symmetries 
bounded and if so, what is their maximal value, or does it possess potential symmetries of 
arbitrarily large orders? In other words, does there exist an integer p such that the equation 
has strictly p-order potential symmetries and all potential symmetries of orders higher than p 
are reduced to potential symmetries of lesser orders? 

Wide classes of linear parabolic equations possessing infinite series of potential symmetry 
algebras, whose sequences of potential orders are unbounded, are constructed in Section [11] 
making use of auto-Darboux transformations. The examples given are connected with the cases 
V = and V = fj,x~^ of the Lie-Ovsiannikov classification. The problem is to study the reduced 
potential equations with V = V{x). An obstacle which should be surmounted is the absence of 
a powerful criterion for the existence of potential symmetries for this case. 

Problem 3. Propose a necessary and sufficient criterion on the existence of potential symmetries 
for the cases when a potential equation possesses a single linearly independent nontrivial Lie 
symmetry operator. 

Theorem 1101 gives only a sufficient criterion on the existence of potential symmetries in such 
cases, in contrast to the cases with a higher number of nontrivial symmetry operators. 

The criteria proposed in Section [9] establish, in fact, a connection between the existence of 
potential symmetries of an equation from class ^ and the reducibility of the equation to a 
special form. The investigation of point equivalence of linear parabolic equation was stimulated 
by the celebrated paper of Kolmogorov [42]. He posed the problem of describing Kolmogorov 
equations (C = 0) which are reduced to the heat equation by point transformations of a special 
form. This problem was completely solved in [18] . A symmetry criterion on the reducibility 
naturally arises in the framework of the Lie-Ovsiannikov group classification [45 i [5T] . It was 
additionally discussed in a number of papers. See, e.g., [Qj [Ml [TT] [TSJ [73]. Nevertheless, the 
symmetry criterion in its present forms is not as constructive as Cherkasov's. 

Constructive criteria on the reducibility can be obtained via the calculation of differential 
invariants and semi-invariants of the equivalence group used. Second-order differential semi- 
invariants of the linear transformations of the dependent variable in class ([T|) were calculated in 
[33j by the infinitesimal method. The same method was used in [39] for finding a necessary and 
sufficient invariant condition on the coefficients of the equations from class ([T|) , which reduced to 
the linear heat equation by point transformations. The more effective approach to calculations 
concerning invariants and equivalence problems in classes of differential equations is given by 
Cartan's method of moving frames [50] in its Fels-Olver version [241 125] . Within the framework 
of the method of moving frames, the equivalence problem for equations from class ([T]) is neatly 
investigated in [47] . 

In [Ij] the generalized Kolmogorov problem on reducibility by combining equivalence and 
Darboux transformations was posed. Note that in fact the term 'Darboux transformation' 
was not used in this paper, as well as the Crum representation [201 HB] of multiple Darboux 
transformations between linear parabolic equations. Only the iterative procedure in terms of a 
sequence of potential systems was presented. A simple application of the Crum theorem does 
not give an exhaustive solution of the generalized Kolmogorov problem. Additional tools should 
be used to create a constructive criterion for the generalized Kolmogorov problem, similar to the 
Cherkasov criterion for the classical Kolmogorov problem and the criterion in terms of differential 
invariants of the equivalence group. 

Problem 4. Does there exist a constructive criterion on equations from class ^ to be connected 
via compositions of point equivalence and Darboux transformations? In particular, is it possible 
to formulate explicit conditions on the arbitrary elements A, B and C under which equation ([T|) 
is reduced by a composition of point equivalence and Darboux transformations to the linear heat 
equation (or to the equation of the form ([7|) with V = ^x^^ or general V = V{x) ). 
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Potential symmetries can be used to construct new exact solutions of equations from class ([T]), 
especially equations possessing no nontrivial Lie symmetries. An obvious way for this is Lie 
reduction of associated potential systems. An alternative but equivalent possibility is based on 
Theorem [9] and other statements of Section [TOj Namely, we can at first find exact solutions 
of the corresponding potential equations and then map them to exact solutions of the initial 
equations by appropriate Darboux transformations. In fact, wide families of exact solutions are 
already known for the linear parabolic equations admitting nontrivial Lie symmetries. That is 
why the second way is preferable. 

The results presented in this paper can be extended to other subjects. For example, they can 
be applied to the investigation of potential symmetries of nonlinear equations which are linearized 
to equations from class ^ (the Burgers equations, u~^-diffusion equation, etc.). Since the field 
(real or complex numbers) in which the dependent and independent variables take values does 
not have an appreciable influence on our investigations, most of the obtained results are easily 
extended to (1 + l)-dimensional linear Schrodinger equations. The nonclassical symmetries of 
the equations from class ^ were described in [53ii55j- We hope that a simultaneous application 
of tools from [29l [53l [SU [67] and this paper will allow us to investigate the potential nonclassical 
symmetries in the class ([T|. 
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